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ABSTRACT 


The  increase  in  terminal  currents  and  noise  observed  at 
high  drain  bias  in  MOS  field-afreet  transistors  is  explained 
by  quantitative  models  that  give  predictions  in  good  agreement 
with  experiment.  We  describe  the  fabrication  and  resulting 
properties  of  a  new  type  of  Schottky  barrier  phototiode.  An 
explanation  is  proposed  for  the  fall-off  in  quantum  yield 
observed  in  metal-semiconductor  photodiodes  for  excitation 
wavelengths  in  the  near  UV  region  of  the  spectrum.  New  results 
are  given  for  the  electron  and  hole  capture  rates  associated 
with  the  defect  centers  produced  by  nickel  in  germanium.  The 
results  of  basic  optical  and  electrical  measurements  are  reported 
for  amorphous  semiconductor  As2Se^Sb2Se^  films  fabricated  by 
evaporation.  Monochromated  Guinier-deWolf f  and  Guinier-Lenne 
x-ray  powder  cameras  and  a  scanning  microdensitometer  are 
combined  into  a  data  system  for  studying  solid  state  reaction 
kinetics.  The  changes  in  the  short  range  ordered  structure  of 
an  Ni^Fe  alloy  are  determined  for  a  series  of  anneals  between 
0  and  40  hours  at  480°C. 
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SUMMARY 


This  report,  for  the  ninth  semi-annual  period  of  contract  support, 
describes  technical  findings  in  two  main  subject  areas:  semiconductors 
and  semiconductor  devices;  and  measurement  techniques. 

SEMICONDUCTORS  AND  SEMICONDUCTOR  DEVICES: 

The  increase  in  terminal  currents  and  noise  observed  at  high  drain 
bias  suggests  the  presence  of  avalanche  multiplication  near  the  drain  of 
MOS  field-effect  transistors.  Quantitative  models  are  developed  that 
predict  values  of  noise  and  terminal  currents  in  good  agreement  with 
experiment.  The  results  demonstrate  that  surface  scattering  strongly 
affects  the  ionization  rates  and  that  one  type  of  carrier  predominates  in 
the  avalanche  process. 

We  describe  the  fabrication  and  properties  of  a  new  type  of  Schottky 
barrier  photodiode.  The  novelty  in  structure  consists  in  the  mask  for 
metallization  which  forms  a  grating  of  gold  film  contacting  the  n-type 
silicon  substrate.  In  spectral  response,  quantum  yield,  speed  and  noise, 
the  device  shows  better  performance  than  that  reported  for  previous  silicon 
Schottky  barrier  diodes. 

An  explanation  is  proposed  for  the  fall-off  in  quantum  yield  observed 
in  metal-semiconductor  photodiodes  for  excitation  wavelengths  in  the  near  UV 
region  of  the  spectrum.  Calculations  based  on  supposing  the  existence  of 
an  inversion  region  bounding  the  metal-semiconductor  interface  agree  well 
with  experimental  observations. 

New  results,  believed  to  be  more  accurate  than  values  reported 
previously,  are  given  for  the  electron  and  hole  capture  rates  associated 
with  the  defect  centers  produced  by  nickel  in  germanium.  In  contrast  with 
previous  studies,  we  deduce  the  rates  using  the  photomagnetoelectric  and 
photoconductive  effects,  taking  careful  account  of  the  effects  of  excess 
carrier  trapping.  Moreover,  we  ensure  sensitive  control  of  the  compensation 
ratio  of  the  nickel  and  the  shallow  impurity  densities. 

The  results  of  basic  optical  and  electrical  measurements  are  reported 
for  amorphous  semiconductor  As2SegSb2Se2  films  fabricated  by  evaporation. 
Transmittance  and  absorption  studies  made  in  the  0.7  to  16  micron  range 
provide  information  about  the  refractive  index,  the  static  relative  dielectric 
constant,  the  relaxation  time,  the  mobility  and  the  thermal  coefficient  of 
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energy.  Measurement  of  the  electrical  conductivity  as  a  function  of 
temperature,  applied  field,  and  illumination  provides  additional  basic 
information  about  transport  processes  and  trapping  levels.  From  the 
understanding  lent  by  the  combination  of  these  measurements,  we  propose 
a  plausible  energy  band  configuration. 

MEASUREMENT  TECHNIQUES: 

The  combined  use  of  monochromated  Guinier-deWolf f  and  Guinier-Lenne 
x-ray  powder  cameras  and  a  scanning  microdensitometer  provides  a  powerful 
tool  for  studying  solid  state  reaction  kinetics.  Application  of  the  room- 
temperature  and  high-temperature  cameras  to  both  interface-controlled  and 
diffusion-controlled  kinetics  is  discussed.  Additionally,  we  describe 
requirements  for  quantitative  sample  preparation,  a  graphical  analysis 
method  and  high  temperature  data  refinement.  An  experimental  comparison 
of  detection  limits  with  other  x-ray  methods  in  qualitative  and  quantitative 
uses  of  the  Guinier  cameras  demonstrates  approximately  a  factor  of  30 
advantage  for  the  Guinier  cameras  in  certain  applications. 

The  changes  in  the  short  range  ordered  structure  of  an  Ni3Fe  alloy 
containing  74  a/o  Ni,  26  a/o  Fe  are  determined  for  a  series  of  anneals 
between  0  and  40  hours  at  480°C.  Three-dimensional  short  range  order 
parameters  are  determined  by  diffuse  x-ray  scattering  analysis  from  a  single 
crystal  sample.  Two  main  conclusions  result.  The  ordered  arrangement 
changes  from  one  dominated  by  nearest  neighbor  interactions  to  one  with 
longer  range  interactions  in  the  period  between  8  and  12  hours  at  480°C. 

The  time  rates  of  change  of  the  magnetic  coercive  force,  hardness,  and 
electrical  resistivity  all  are  approximately  equal. 
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I.  Introduction 

The  original  general  objective  of  this  research  program  was  to 

establish  at  the  University  of  Florida  a  "Center  of  Competence  in 

Solid  State  Materials  and  Devices"  capable  of  doing  research  of  benefit 

to  the  Department  of  Defense.  At  this  stage,  four  and  one  half  years 

following  initiation  of  contract  support,  ample  evidence  exists  that 

the  University  of  Florida  has  indeed  become  such  a  center  of  competence. 

Among  the  evidence  supporting  this  conclusion  are  the  technical  findings 

.ostered  by  this  contract.  These  are  documented  in  the  eight  previous 

1—8 

reports  in  this  series  and  in  109  related  papers  published  in  research 
journals. 

During  the  past  year  the  objectives  of  the  research  program  have 
shifted  decisively.  In  accordance  with  an  agreement  made  with  the 
contract  monitor  in  January  1971,  we  have  focussed  the  objectives  to  give 
greatest  emphasis  to  the  following  problem  areas: 

1.  The  development  of  new  types  of  detectors  with  stress 
given  to  utilizing  the  Dember  and  photomagnetoelectric 
effects  in  semiconductors  and  to  utilizing  new  mask  geometries 
for  Schottky  barrier  detectors  compatible  with  conventional 
silicon  technology. 

2.  The  incorporation  of  the  effects  of  irradiation  into 
the  characterization  of  device  behavior;  junction  field- 
effect  devices  are  to  receive  special  attention.  Emphasis 
is  placed  on  the  development  of  characterizations  enabling 
computer  simulation  and  design  of  circuits, 

3.  Full  static  and  dynamic  characterization  of  bipolar  and 
field— effect  devices  aiming  both  toward  improved  device  design 
and  toward  improved  circuit  implementation.  The  characterization 
is  to  deal  with  operation  over  a  wide  range  of  temperatures, 
with  emphasis  on  the  range  between  ordinary  room  temperature 

and  cryogenic  temperatures. 

4.  Noise  studies  proceeding  collaterally  with  each  of  the 
efforts  listed  above, 

5.  Radiation  studies  of  amorphous  semiconductors. 

For  the  most  part,  the  content  of  the  present  report  reflects  this 
change  in  emphasis.  Section  II,  concerning  semiconductors  and  semiconductor 
devices,  describes  studies  in  detectors,  field-effect  devices,  noise,  and 
amorphous  semiconductors.  Section  III,  however,  which  deals  with  measure¬ 
ment  techniques,  describes  the  results  of  research  begun  before  the 
implementation  of  our  change  in  emphasis. 
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H*  Semiconductors  and  Semiconductor  Devices  (E.  R.  Chenette, 
S.  S.  Li,  F.  A.  Lindholm  and  A.  van  der  Ziol) 


A.  NOISE  PERFORMANCE  OF  METAL-OXIDE-SILICON  FIELD  EFFECT 
TRANSISTORS  (C.  S.  Kim,  E,  R.  Chenette,  A.  van  der  Ziel) 


1 .  INTRODUCTION 

Since  the  Metal-Oxide-Semiconductor-Field-Effect-Transistor  (MOSFET) 

was  first  introduced,  much  work  has  been  done  on  its  basic  device  charac¬ 
teristics  and  noise -1’ 2  With  the  advances  in  fabrication  technology, 
large  scale  MOS  integrated  circuits  have  become  popular  and  a  new 
memory  device,  called  the  charge-coupled  MOS  device,  has  been  intro¬ 
duced.  MOS  devices  are  widely  favored  over  bipolar  devices  in  cryogenic 
applications  and  in  many  high  frequency  communication  circuits.  It  is 
probable  that  research  on  MOS  devices  will  continue  to  provide  for  the 
development  of  new  and  important  FET  applications. 

Avalanche  ionization  occuring  in  the  channel-to-drain-depletion 
(CDD) J  region  of  an  FET  has  recently  received  considerable  attention. 

This  phenomenon  was  first  noticed  in  the  form  of  the  excess  gate  cur¬ 
rent  in  junction  gate  FET’S  operating  at  low  temperatures.*^  The  anom¬ 
alous  noise  enhancement  observed  when  operating  at  low  temperatures  of 
far  beyond  saturation  has  been  credited  to  this  avalanche  ionization. 
Klaassen  drew  from  his  experiments  the  qualitative  conclusion  that  the 
noise  associated  with  the  excess  substrate  current  in  a  MOSFET  resulted 
from  avalanche  ionization  was  a  superlinear  type  of  shot  noise.7  The 
expression  he  presented  for  this  noise  is 


*sub^  ~  ^  q  ^sub 


(1.1) 
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where  M(lgub)  is  the  current-dependent  shot  noise  multiplication 
and  1  k  is  the  average  substrate  current. 

The  purpose  of  the  present  work  is  to  give  a  clear  physical  pic¬ 
ture  of  how  this  avalanche  ionization  in  the  CDD  region  of  a  FET 
occurs.  A  quantitative  model  for  the  noise  associated  with  the  ava¬ 
lanche  current  was  developed;  it  compared  favorably  with  our  experimental 
results. 

The  study  of  avalanche  mechanism  in  FET*s  led  to  a  convenient 
method  of  measuring  avalanche  ionizarion  rates.  The  previous  exper¬ 
imental  method  on  a  two-terminal  device  involved  the  necessity  of  dis¬ 
tinguishing  between  two  currents  of  nearly  identical  magnitudes.  This 
limited  the  accuracy  of  the  measurement.  Furthermore  the  evaluation 
of  the  ionization  rate  of  majority  carriers  from  the  data  required  a 
proper  correction  for  the  contribution  to  the  avalanche  current  from 
minority  carriers. 

These  difficulties  are  eliminated  in  experiments  with  the  avalanche 
current  in  four- terminal  devices  such  as  FET's,  due  to  the  uniqueness 
of  its  avalanche  mechanism.  The  reason  is  that  the  minority  carriers 
are  collected  by  the  gate  (in  a  JFET)  or  by  the  substrate  (in  a  MOSFET) ; 
this  considerably  simplifies  the  avalanche  picture. 

The  contribution  of  this  work  is  establishing  a  clear  physical 
picture  of  the  avalanche  ionization  mechanism  in  a  MOSFET  and  sub¬ 
stantiating  it  with  noise  study. 
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II.  AVALANCHE  IONIZATION  AND  RELATED  NOISE  IN  A  MOSFET 


The  drain  terminal  current  of  a  MOSFET  at  the  onset  of  saturation 
can  be  expressed  in  terms  of  terminal  voltages  and  other  device  para- 

g 

meters  in  the  following  way. 

=  CHoffCoxZ/lO  {VVDS-vDS2-(4/3)fFvB  ^(VDS)3/2-  S(0)3/2]  J  (2.1) 

Here,  B(V)  =  1  +  (V  -  VgUB)/  2<Pp  and 

v  '  =  V  -lb  -  4  '  +  V 
G  G  ’so  ^MS  ss 

Other  symbols  are  defined  as  follows: 

y  ^  =  effective  mobility 

C  =  oxide  capacitance 
ox 

L  =  channel  length 

Z  =  channel  width 

c|>  =  bulk  Fermi  potential 

r 

ibsQ  =  surface  potential  at  source  junction 

<J>  ,=  work  function  difference  between  metal  and  semiconductor 

MS 

V„  =  gate  terminal  voltage 

U 

V  =  substrate  terminal  voltage 

SUJd 

V  =  drain  terminal  voltage  at  the  onset  of  saturation 

Uu 

V  =  Q  /C  =  surface  and  oxide  charge  in  equivalent  voltage, 
ss  ss  ox 
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In  deriving  the  above  equation  the  mobile  charge  density  in  the 
conducting  channel  is  usually  assumed  to  depend  only  on  the  component 
of  electric  field  which  is  normal  to  the  surface.  (Gradual  Channel 
Approximation)  This  assumption  deteriorates  in  the  vicinity  of  drain 
junction  as  the  MOSFET  is  biased  beyond  saturation,  and  mobile  carriers 
are  depleted  from  this  region. 

The  FET  is  usually  divided  into  two  regions  along  the  channel 
in  order  to  describe  the  saturation  characteristics;  one  extends  from 
the  source  to  the  virtual  drain  where  the  gradual  channel  approximation 
is  still  valid  and  the  other  from  the  virtual  drain  to  the  actual  drain 
where  the  depletion  approximation  can  be  used.  Then  the  channel  length 
and  the  saturation  drain  voltage  in  Eq.  (2.1)  must  be  replaced  by  the 
length  from  the  source  to  the  virtual  drain  and  the  actual  channel 
potential  at  the  virtual  drain,  respectively,  to  obtain  the  saturation 
drain  current  . 

When  the  drain  is  biased  well  beyond  saturation  the  electric 
field  along  the  direction  of  electron  flow  in  the  channel-to-drain- 
depletion  (CDD)  region  becomes  sufficiently  high  to  cause  avalanche 
ionization.  Avalanche  ionization  can  occur  before  the  depletion  region 
reaches  through  to  the  source  junction  or  before  the  drain-substrate 
diode  breaks  down  in  the  reverse  direction. 

Figure  (2.1)  shows  a  sketch  of  the  CDD  region  of  a  FET.  It  is 
assumed  that  the  drain  current  is  carried  in  a  thin  slab  near  the  inter¬ 
face  of  the  oxide  and  the  semiconductor.  It  is  not  necessary  to  con¬ 
sider  the  effect  of  the  transverse  electric  field  in  detail.  The 
transverse  electric  field  does  not  affect  the  avalanche  process  directly; 
it  only  serves  to  remove  minority  carriers  generated  when  ionizing 
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gate 


i 


id=  electrons  collected  at  drain  terminal 
ic, .wholes  collected  at  substrate 

vlt  I J 


ir=  electrons  entering  channel  from  source 
terminal. 


Figure  2.1--Sketch  of  CDD  Region 
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collisions  occur.  This  is  probably  the  most  important  difference 
between  the  avalanche  process  as  it  occurs  in  the  CDD  region  of  a  FET 
and  in  the  avalanche  diode.  With  this  device  it  is  possible  to  study 
avalanche  multiplication  of  electrons  alone  or  of  holes  alone.  The 
simple  one-dim emsional  model  of  Fig.  (2.1)  provides  a  useful  description 
of  the  processes  involved. 

Consider  an  n-channel  MOSFET .  Here  the  current  is  carried  by 
electrons.  The  channel  current:  from  the  source  region,  Ig,  enters 
from  the  left  hand  side  and  passes  through  the  CDD  region,  w,  is  smaller  than  a 
Debye  length  for  most  of  the  current  densities  of  interest  the  secondary 
holes  are  scattered  out  of  the  CDD  region.  These  holes  are  collected 
at  the  substrate  terminal  as  Igub  because  of  the  transverse  component 
of  the  electric  field.  Both  primary  and  secondary  electrons  are  col¬ 
lected  at  the  drain  terminal  and  cause  the  drain  current  1^. 

The  sign  of  the  current  is  taken  positive  in  the  direction  of 
electron  flow.  Obviously  we  have 


Id  =  ^s  ■*"  ^sub 


(2.2) 


The  conventional  multiplication  factor  is  defined  as  the  ratio  of  the 
average  total  electron  current  emerging  out  of  the  CDD  region  to  the 
electron  current  entering  it. 

H  »  ij  /  il  »  Id/  Is  (2‘3) 

The  secondary  holes  do  not  participate  in  the  avalanche  multi¬ 
plication  because  the  transverse  component  of  the  electric  field  is  not 


8 


sufficiently  high  to  cause  ionizing  collisions.  If  we  neglect  the 
thermal  generation  and  recombination  of  carriers  in  the  CDD  region,  the 
current  increase  in  an  infinitesimally  small  region  between  y  and  y+dy 
can  be  described  by  the  following  differential  equation. 


di(y)  =  Ky)a(y)  dy  (2.4) 

a 

where  (y)  is  the  ionization  rate,  defined  as  the  number  of  ionizations 
prt  unit  length.  Solving  the  above  equation  with  a  proper  set  of 
boundary  conditions  we.  have 


I(y)  =  Is  exp  (/*  ct(y')  dy’) 
Id  =  Is  exp  (/  a(y)  dy) 


(2.5) 


The  multiplication  factor  can  thus  be  written  as 

M  =  exp  (/£  Oi(y)  dy)  (2.6) 

A  firsv  order  approximation  of  the  electric  field  distribution 
in  the  CDD  region  can  be  obtained  from  the  one-dimensional  solution 
of  Poisson's  equation. 

dE(y)/dy  -  q  1^/K,  E()  (2.7) 

where,  q  =  electron  charge 

=  acceptor  impurity  density 

K  =  dielectric  constant  of  silicon 
s 

£q  =  permittivity  of  free  space 

The  CDD  region  starts  from  the  virtual  drain  located  at  y  =  i. 

The  electric  potential  at  the  virtual  drain  can  be  approximated  by  the 
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drain  saturation  voltage  .  The  solution  of  Eq.  (2.7)  gives 
E(y)  =  Es  +  a  (y  -  A) 


Ed  =  Es  +  a  d  =  a  d 


where  Ed  >>  Eg  is  assumed  (2.8) 


d  =  L  -  l  =  *-U2  ^  .  Vds)1/2 
where  the  boundary  values  and  constant  a  are  defined 


as 


E<L)  e  Ed,  E(£)  =  ESJ  V (L)  =  Vd,  V(l)  =  Vds 
a  =  lN a/  Ks  eQ 


(2.9) 


The  field  dependence  of  the  ionization  rate  for  the  case  of  a  small 
electric  field  is  well  approximated  by 

a(E)  =  A  exp  (-B/E)  ^ 

where  A  and  B  are  constants  determined  by  the  scattering  parameters 
and  the  threshold  energy  for  the  avalanche  ionization.  In  the  MOSFET 
case  surface  scattering  plays  a  dominant  role  in  decreasing  the  kinetic 
energy  of  the  electron  gained  from  the  electric  field.  How  this  addi¬ 
tional  mechanism  of  surface  scattering  will  affect  the  field  dependence 
of  the  ionization  rate  in  Eq.  (2.10)  is  not  available  in  the  literature. 
For  a  value  of  M  close  to  unity,  Eq.  (2.6)  can  be  approximated  by 

-  d  d 

M  =  exp  (fQ  ot(E)  dy)  *  1  +  /  a(E)  dy  (2>U) 

Combining  Eq.  (2.2),  (2.3),  (2.8),  (2.10)  and  (2.11),  the  substrate 
current  can  be  written  as 
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^sub  (W  -  1) 


h  fc  ME)  dy 


MM 

*  Is  A  exp  (-B/E)  dy 
Ed 

rs  (MM  -(g  exp  (-B/E)  dE,  where  dy  =  a”1  dE 


r 

Is  (A/a)  /  j  exp  (-Bz)/z2,  where  E  =  ]./z, 


dE  =  -z  dz 


Is  CAM)  (Ed  S*p  (-B/E,)  -  Eg  exp  (-B/KJ  +  B  E±  (B/Ed) 


"  B  Ei  (B/Ed) ] 


(2.12) 


where  the  symbol  (x)  refers  to  the  exponential  integral  defined 


E-^  (x)  =  /”  exp  (-t)/t  dt  with  x  >  0 


(2.13) 


If  we  assume  E,  »E  ,  we  obtain 
a  s 


Isub  =  Is  (A/a)  Ed  exP  (-B/Ed> 

=  Is  A  a"1/2  (Vd  -  Vds)M2  exp  t-2B/a1/2  (Vd  -  Vdg)1/2]  (2 

The  calculation  of  the  integral  in  Eq.  (2.12)  enables  us  to  rewrite 
Eq.  (2.6)  as 

In  M  -  (A/a)  ^  exp  (-B/Ed)  -  Es  exp  (-B/Eg)  +  B  E-  (B/Ed) 


"  B  Ei  (B/Eg)] 

"  (A/a)  fed  e*P  (-B/Ed)  +  B  E£  (B/Ed) ] 
where  E^  »  Eg  is  again  assumed. 


(2.15) 
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Summarizing  the  result  we  have  the  following  expressions  for  the 
terminal  currents. 

rs  =  ^effCox^1)  {vG'vds  "  Vds2  "  (4/3)<f>pVB  [3(Vds)3/2  -  6(0)3/2]} 
rsub  =  ISA  &  ^  -  V^)1/2  exp  [ -2B/  a1^2^  -  V^)1/2] 

Id  =  Is  ^  1/2^d  “  Vds^^2  exp  ^-2B/  fll/,2^vd  "  V^)1^2]}  (2.16) 

The  multiplication  noise  for  the  case  of  avalanche  diodes  has  been 
well  formulated  by  many  authors.  ’  The  avalanche  multiplication 
noise  in  MOSFET's  has  to  be  treated  differently  from  that  in  p-n 
junctions  for  two  main  reasons.  First,  only  one  type  of  current  carrier 
(electrons  in  n-channel,  holes  in  p-channel  devices)  is  active  in  the 
multiplication  process.  Secondly,  the  current  entering  the  multipli¬ 
cation  region  does  not  exhibit  fi’1!  shot  noise.  This  is  obvious  in 
the  sense  that  the  channel  noise  in  FET’s  in  its  limiting  form  is 
thermal  noise  proportional  to  the  channel  conductance. 

The  multiplication  process  is  a  random  process.  Therefore,  when 
one  electron  entering  the  multiplication  region  (in  this  case,  the  CDD 
region)  causes  M  electrons  to  emerge  out  of  it,  M  is  a  random  integer. 
Since  each  electron  goes  through  an  independent  random  process  we  can 
write  for  the  total  number  of  electrons  collected, 

ns 

n.  =  I  M.  (2-17) 

d  i=l  i 

where  ng  and  M.^  are  both  fluctuating  quantities  in  the  time  domain. 

In  calculating  nd>  the  average  number  of  electrons  collected  at 
the  drain  terminal,  we  first  set  up  a  subensemble  of  unit  time  elements 


during  which  the  same  number  of  n  primaries  enter  the  multiplication 

s 

region.  We  take  the  average  over  this  subensemble,  i.e.. 


nd  "  ns  M 


(2.18) 


Then  we  set  up  the  entire  ensemble  consisting  of  all  these  subensembles, 
and  take  the  average,  i.e., 


— s 


n ,  =  n .  =  n  M  =  n  M 

d  d  s  s 


(2.19) 


Applying  the  same  concept  involved  in  Eq.  (2.18)  and  (2.19),  we 

can  calculate  the  variance  of  n,  as  follows: 

d 


ns 


ns 


Hg  ns 


"0  ■  ( A  v  (j-i  V  -  i=i  A  Mi  "i 


(2.20) 


In  Eq.  (2.20)  the  double  summation  results  in  n  (n  -  1)  terms  with 

s  s 

different  i,  j  and  in  ng  terms  with  same  i,  j.  Hence, 


nd  =  ns  (ns  "  1)  MiMj  +  ns 


,  2 


“  ns(ns  ”  1)  (Mi)2  +  ng  M ^ 


n  Var  M  +  n  2  (M)2 


(2.21) 


Therefore, 


nd2  =  ns  Var  M  +  ns2  (M)2 


and. 


Var  nd  =  ng  Var  M  +  ns2  (M)2  -  (ng)2  (M)2 
=  Var  ns  (H)2  +  Ng  Var  M 


(2.22) 
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The  derivation  of  Eq.  (2.22)  from  Eq.  (2.17)  actually  furnishes 

proof  of  the  so-called  Burgess'  Variance  Theorem.  According  to  another 

theorem  (van  der  Ziel,  Noise;  Source  ■,  Characterization,  Measurement, 

Prentice-Hall,  Inc.,  Englewood  Cliffs,  N.  J.,  1970,  p.  15)  the  low 

frequency  spectrum  of  the  random  variable  n,  is 

a 


Snd  (0)  =  2  Var  nd 

Combining  Eq*.  (2.22)  and  (2.23),  the  noise  spectrum  of 
rent  i,  will  be 

<3 


(2.23) 


the  drain  cur- 


sid  (0)  =  q  2  Var  nd 


2  —  2 

=  2  q  (M)  Var  n  +  N  Var  M 

o  S 


(2.24) 


We  can  represent  these  noises  by  their  corresponding  equivalent 

noise  diode  currents,  1  =  S(0)/  2q. 

eq  n 

Xeqs  =  Sis  (0)/  2q  =  q  Var  ns  (2.25) 

Ieqd  =  Sid  (°)/  2q  =  hqs  ^  +  h  Var  M  (2.26) 

where  ^  =  q  Ng.  Note  that  Ieqg  £  2  q  Ig.  In  other  words,  the  electron 
current  entering  the  multiplication  region  does  not  show  full  shot 
noise  as  in  the  case  of  the  avalanche  diode;  instead  it  shows  the  low- 
field  channel  noise  of  the  FET. 

lae  next  step  is  to  calculate  the  variance  of  M.  Let  us  assume 
that  the  multiplication  process  is  a  random  process  with  its  probability 
sensity  obeying  Poisson's  law.  This  assumption  should  be  valid  since 
the  individual  ionizing  events  are  independent.  We  can  write  the  prob¬ 
ability  that  an  electron  entering  the  multiplication  region  will  produce 


“  lonization  events  as 
pm  c  u  exp  C-u)/  mi 

(2.27) 

where  u  -  J  =  J  d>  ls  . 

“  *•  «.e  averaee  lonlzatlo:;;:r0rr;fi-i2-»”  «* 

If  we  define  P(M)  as  th  e^r°n  per  un*t  distance. 

^  the  drain  ter®inai  due  to  onTeLcT"  ^ 

region,  We  haVe  ectro"  entering  the  multipliCation 

*  PG  *  exp  C-u) 

(2.28) 

in  other  words  pm  i c 

pass  through  the  ““  a"  Neutron  win 

ionizing  collision.  8100  Wlthout  going  through  an 

When  one  elec 

P«h  through  the  multipart Z°Tell7  t°nl2ln8  COlUSl°"  dur1"*  its 
0f  d™iu  terminal  and  one  hole  „  U  eleCtr°"S  "U1  —»•  - 
terminal .  Now  ,(2)  has  t„  he  ca!c„ l  f  a  at  the  substrate 

^  «“  -ting  collision  a  I  ^ 

the  probability  of  both  primary  and  s  7‘  the  <>“» 

the  rest  of  the  path  without  having  any  i  Tt  eleCtr0"S  — "g 

summing  this  product  over  the  probabili  COUlsl°"  a"d  ««■  by 

initial  ionization  at  different  p„i„ts  7  dlStrIbuti»”  »f  having  the 

“ne  of  reasoning  we  can  get  the  genii  e7  ^  Path'  fcU»—  tbla 

expression  of  p(m)  as 


P(*°  =  C«  ex?  ("Mu)  fexp(u)  -xjM  -  1 
where  is  a  constant  coefficient. 


(2.29) 


Obviously, 


and, 


(2.30) 


FT  fe'W1  “h!^  p((0  *  0  (2-31) 

Combining  Eq.  (2.29),  (2.30)  and  (2 . 31) .yields 
M  =  exp  (u) 


=  2  (M)2  -  M 


Var  M  =  M2  -  (M)2  =  M  (M  -  1) 


(2.32) 


The  detailed  derivations  of  Eqs.  (2.29)  and  (2.32)  are  presented  in  the 
appendix.  Eq.  (2.32)  was  already  given  by  Tager. 

Substituting  Eq.  (2.32)  into  Eq.  (2.26)  and  using  Igub  -  Ig  (M  -  1) , 
we  get. 


W  ‘  Ieqs  M2  +  ^ub  M 


(2.33) 


In  the  appendix  the  detailed  calculation  of  the  noise  in  the  substrate 
terminal  current  and  the  correlated  part  in  the  drain  noise  and  the 
substrate  noise  is  carried  out.  The  result  is 

Wub  ’  W  W  -  1)2  +  rsub  5  <2'34> 

leqd.sub  *  Ieqs  ««  -  1)  +  ISub  S  <2-35> 

The  above  derivation  was  based  on  a  corpuscular  picture  of  the 
avalanche  noise  process.  The  same  result  can  be  obtained  by  a  more 
collective  approach.  Rewriting  Eq.  (2.4),  including  the  time  dependence 
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and  solving  it  we  have 


i(y,t)  =  1(0, t)  exp  (fVQ  ot(y ' , t)  dy.)  (2<36) 

Xd  =  is  exP  C/c  ct(y,t)  dy)  =  is  M(t)  (2  37) 

Since  there  is  no  physical  reason  for  spatial  correlation  in  the 
fluctuation  of  the  avalanche  ionization. 


di(y>t),  di(y',t)*  =  0 


(2.38) 


If  we  assume  that  the  avalanche  ionization  in  the  infinitesimally 
small  subregion  is  a  pure  random  process,  we  can  expect  that  each 
incremental  current,  di,  will  show  full  shot  noise  with  a  white  spectrum. 

di(y,t)2  =  dS-  df  =  2  q  dl(y)  df  (2.39) 

Substituting  Eq.  (2.4)  and  (2.5)  into  the  above  equation  we  have 

dS^y)  =  2  q  Is  a(y)  exp  [/  a(y')  dy*  ]  dy  (2.40) 

This  shot  noise  current  is  amplified  through  successive  avalanche  multi¬ 
plications  while  flowing  the.  rest  of  the  distance,  d-y,  toward  the  drain 
terminal.  It  appears  at  the  drain  terminal  as, 

dSi(d)  =  dSj,(y)  exp  [2/^  a(y')  dy'] 

=  2qlg  a(y)  exp  [2 a(y')  dy«]  exp  [/£  a(y")dy"]  dy  (2.41) 
Summing  over  all  these  infinitesimal  shot  noise  generators,  we  have 
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si  -  ft  dvd) 

"  2qIs  {exP  [/o  a&">  dy'3>2  /o  exP  ^  “Cy")  dy"]  a(y)  dy 
“  2qls  (M)“  JQ  M  exp  (-z)  dz 
"  2qlg  (M)2  (1  -  1/M) 


"  2qlg  M  (K  -  1) 


(2.42) 


where  z  =  fy  a(y)  dy,  and  dz  =  ot(y)  dy. 

This  only  gives  the  additional  noise  due  to  the  fluctuation  in 
the  avalanche  multiplication.  In  other  words  we  have  only  considered 
the  fluctuation  of  the  exponential  part  of  Eq.  (2.37).  The  fluctuation 
of  i  was  not  considered.  If  we  adH  the  noise  of  i  we  obtain 

S  Q 


-  2  _  _ 

00  +  2qls  K  (M  -  1) 


(2.43) 


In  a  similar  way  we  can  derive  Eqs.  (2.34)  and  (2.35). 

It  is  not  surprising  that  both  approaches  give  identicalresults, 
for  both  are  based  on  the  assumption  that  the  multiplication  process 
is  a  random  process  obeying  Poisson's  law  which,  in  turn,  implies 
shot  noise. 

Eqs.  (2.33),  (2.34)  and  (2.35)  summarize  the  noise  behavior  of 
the  multiplication  process  in  FET's  in  a  simple  form.  We  can  make  some 
interesting  observations  while  studying  these  equations. 

First,  the  noise  in  the  primary  current  is  amplified  during  the 
multiplication  process,  whether  the  noise  is  flicker  noise  or  thermal 
noise.  At  low  frequencies  the  flicker  noise  predominates  and  the 
noise  at  the  drain  is  an  amplified  thermal  noise. 


i 


l 


| 
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The  noise  resulting  only  from  the  fluctuation  in  the  multiplica¬ 
tion  is  a  shot  noise  with  a  white  spectrum.  The  magnitude  of  this 
noise  is  Igub  M,  where  both  Igub  and  M  are  sensitive  functions  of  drain 
/oltage.  At  high  frequencies  this  component  of  noise  predominates  and 
the  amplified  thermal  noise  is  usually  buried  below  this  noise  level 
when  there  are  appreciable  multiplication  currents. 

Finally,  the  multiplication  noise  in  FET’s  is  in  general  smaller 
in  magnitude  than  that  of  an  avalanche  diode.  The  reason  is  due  to 
the  difference  pointed  out  earlier  i.e.,  only  one  type  of  carrier  is 
active  in  the  multiplication  process  and  the  primary  current  does  not 
display  full  shot  noise,, 

When  both  types  of  carriers  go  through  multiplication  the  fluctu¬ 
ation  of  M  is  much  greater.  If  the  ionization  rates  are  assumed  to  be 

equal  for  both  electrons  and  holes,  the  variance  of  M  in  the  avalanche 
diode  is 

Var  M  =  (if)2  (M  _  i) 

(2.44) 


If  the  channel  current  entering  the  multiplication  region  shows  full 
shot  noise. 


(2.45) 


Then  substituting  Eq.  (2.44)  and  (2.45)  into  (2.33)  gives 


^eqd  “  Is  (M)2  +  Is  (M)2  (fl  -  1) 

-  Is  W)3 


(2.46) 


Eq.  (2.46)  is  identical  to  the  formula  for  multiplication  noise  in 

avalanche  diodes  vased  on  equal  ionization  rates  for  electrons  and  holes 

obtained  by  others.  For  example,  one  can  obtain  Eq.  (2.46)  by  sub- 

1 2 

stituting  k=l  into  the  formula  presented  by  McIntyre.  Here  k  is  the 
ratio  of  hole  ionization  rate  to  electron  ionization  rate. 
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III.  EXPERIMENTAL  METHODS 


A.  Device  Fabrication 

Both  n-channel  and  p-channel  si-MOSFET's  were  fabricated  on  a 
1  ohm-cm  substrate  with  (1,0,0,)  orientation.  Fig.  (3.1)  shows  the 
geometry  of  the  devices.  Devices  with  channel  lengths  of  1,  2  and 
^  ® ^ ,c  were  made  on  the  same  wafer.  The  gate  was  allowed  to  overlap 
about  one  tenth  of  a  mil  on  each  side  toward  the  source  and  the  drain. 

The  wafer  was  stripped  of  the  oxide  after  diffusion  of  the  source 
and  the  drain.  The  gate  oxide  was  regrown  by  placing  the  wafer  first 
in  dry  oxygen  for  30  minutes,  then  in  wet  oxygen  for  15  minutes;  and 
finally  returning  it  to  a  dry  oxygen  environment  for  one  hour.  It 
was  then  placed  in  a  nitrogen  furnace  at  900°  C  for  annealing.  After 
etching  away  the  first  few  hundred  Angstroms  of  oxide  from  the  top, 
the  wafer  was  boiled  in  hot  hydrochloric  acid.  It  was  thoroughly 
rinsed  and  placed  in  a  phosphorus  diffusion  furnace  for  about  ten  min¬ 
utes  to  grow  a  phosphor-silicate  glass  layer. 

The  back  side  of  the  substrate  was  also  doped  heavily  to  form  a 
degenerate  layer  which  provides  a  good  ohmic  contact.  Gold  was  used 
to  mount  individual  chips  on  the  header. 


B.  Device  Characterization 

The  capacitance  versus  bias  voltage  mtasurements  on  the  MOS 
capacitors  which  were  fabricated  close  to  each  device  site  were  per¬ 
formed  to  determine  the  oxide  thickness,  flat  band  voltage  and  bulk 
13 

doping.  Fig  (3.2)  shows  a  typical  C-V  (Capacitance  vs.  Bias  Voltage) 
curve  measured  at  1  MHz.  The  high  frequency  capacitance  minimum  gives 
information  on  the  density  of  the  bulk  impurity  which  contributes  to 
the  depletion  layer  capacitance,  the  minimum  capacitance  being  the 
series  combination  of  the  oxide  capacitance  and  the  depletion  layer 
capacitance. 

The  channel  conductance  was  plotted  on  an  x-y  recorder  as  a 
function  of  the  drain  bias  with  the  Wayne  Derr  Model  B-641  Autobalance 
Bridge  operating  at  1,592  Hz.  According  to  the  elementary  gradual 
channel  solution  these  quantities  show  a  linear  relationship  described 
by  the  following  equation. 

Sd  -  c*ffco*z/L>  <V  -  vd>  ai) 

where  the  symbols  are  as  defined  in  Chapter  II. 

The  channel  conductance  ploted  as  a  function  of  drain  bias  exhib¬ 
its  a  good  linear  relationship  as  shown  in  Fig.  (3.3).  The  effective 
gate  voltage  V^'  lumps  the  effect  of  bulk  charge,  oxide  charge  and  the 
interface  surface  states..  The  effect  of  the  substrate  bias  voltage 
is  also  taken  care  of  by  the  effective  gate  voltage.  The  effective 
gate  voltage  was  determined  by  fitting  a  straight  line  to  the  experi¬ 
mental  plot  and  finding  the  intersection  points  on  the  horizontal  scale. 
Vg’  thus  obtained  is  equivalent  to  the  drain  saturation  voltage  . 
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MOS  CAmCITANCE  2S  6IAS 
DEVICE  No.  36 


Figure  3.2— Typical  C-V  Curve 


DRAIN  CONDUCTANCE  3Z£  DRAIN  VOLTAGE 
DEVICE  No.  38- L 
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Figure  3.3— Typical  Drain  Conductance  vs.  Drain  Voltage  Curve 


The  slope  of  this  line  which  corresponds  to  the  factor  u  C  Z/L  can 
,  ,  ,  eff  ox 

be  used  to  determine  the  effective  mobility  by  substituting  the  other 

parameters. 

Finally  the  terminal  currents  (source,  drain  and  substrate  cur¬ 
rents)  were  simultaneously  recorded  as  a  function  of  the  drain  bias  on 
an  x-y  recorder.  The  bias  was  varied  from  zero  to  where  the 
effect  of  avalanche  ionization  in  the  CDD  region  appears  obvious.  The 
substrate  current  was  separately  monitored  with  a  nanoampere  meter.  The 
average  value  of  the  multiplication  factor  M  was  determined  from  this 
measurement  using  the  following  relation. 

M  =  Id  /  Is  =  (I,  +  Isub)/Is  =  1  +  Isub/Is  (3,2) 


C.  Noise  Measurement 

The  low  frequency  portions  of  the  noise  current  spectra  were  mea¬ 
sured  with  the  G.  R.  1921  Real  Time  Spectrum  Analyzer  and  the  310  AH.  P. 
spectrum  analyzer.  The  calibration  was  achieved  by  a  G.  R.  Random  Noise 
Generator  whose  spectral  intensity  was  determined  by  graphically  inte¬ 
grating  the  output  voltage  spectrum  over  the  entire  frequency  range. 

The  high  frequency  port5.cn  of  the  noise  spectrum  was  measured  by 
the  Hanbury  Brown  Twiss  type  correlation  circuit  recommended  by  Chen 
and  van  der  Ziel.  The  system  consists  of  two  modified  Collins 
Communication  Receivers  and  Cross  Correlator  circuit  built  by  Kalter.1^ 
The  system  set-up  is  described  in  Fig.  (3.4)  The  effective  noise 
band-width  is  determined  by  the  intermediate  frequency  filter  of  the 
receivers. 

The  noise  current  source  was  transformed  into  a  noise  voltage  source 
through  the  source  follower.  The  multiplier  vd  integrator  circuit 
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Figure  3.4— Cross  Correlation  System  Set-up 


performs  the  following  calculation. 

I^77Wvs':  Vb  +_v? +  %v  + 

=  G  G,  2 
a  b  e 

n 

Here  enea*  =  e^*  =  e^*  =  0,  because  the  amplifier  noise  0“  one 
channel  is  not  correlated  with  that  of  the  other  channel  and  the  device 
noise  is  not  correlated  with  the  amplifier  noise.  The  output  of  the 
integration  circuit  therefore  responds  only  to  the  device  noise.  It 
was  possible  to  detect  noise  currents  down  to  I  =  2  jj  amp  expressed 
in  terms  of  equivalent  saturated  diode  current,  when  a  high-Q  tuned 
circuit  was  used  as  a  load  for  the  device# 

The  device  was  connected  in  common  substrate  (in  analogy  to  common 
gate)  configuration  with  the  gate  and  substrate  a.c.  short-circuited 
to  ground.  Fig.  (3.5)  shows  the  noise  equivalent  circuit. 

If  we  neglect  the  thermal  leakage  ciirrent  in  the  reverse  biased 
source  to  substrate  diode  and  in  the  drain  to  substrate  diode,  the 
substrate  current  results  only  from  the  multiplication  process.  There¬ 
fore  the  substrate  current  noise  is  fully  correlated  with  the  drain 
noise  current.  The  short-circuited  output  noise  current  is  given  by 


i  2  _  ^  _ 

out  ~  |Ys  +  Ysub  +  Ymp  (  4kT  gs  A  f 


x  — iw  ,  „  *  sub  Yin  „ 

|Ym|2  I  s  Ysub  +  ^  |2]  (3.3) 

Since  the  transconductance  associated  with  the  substrate  terminal 

is  much  greater  than  the  admittance  of  the  high-Q  tuned  circuit  at  the 

source  terminal  or  the  admittance  of  the  source-substrate  diode  biased 
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Figure  3.5  Equivalent  Noise  Circuit  of  Common  Substrate  Configuration 


(3.4) 


in  the  reverse  direction,  we  can  impose  a  condition, 

*in  >>  Ys,  ^  >>  Yguk 

This  condition  reduces  Eq.  (3.3)  to 

W2  =  4kT  Ss  4  £  +  isub2 

=  4kT  gs  A  f  +  ig2  (M  -  l)2  +  2q  Isub  M  A  f  (3.5) 

When  the  source  terminal  is  a.c.  short-circuited,  the  magnitude  of 
the  output  noise  is  simply 

7  2  _  7T 

*-out  1d 

=  is2  (M)2  +  2q  Isub  M  A  f  (3.6) 

The  drain  current  noise  was  also  measured  with  all  the  other 
terminals  a.c.  short  circuited  to  ground.  The  same  measurement  was 
done  for  substrate  current  noise  with  all  the  other  terminals  a.c. 
short-circuited  to  ground. 
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IV  RESULTS  AND  DISCUSSIONS 


The  basic  device  parameters  were  determined  by  use  of  the  expert 
mental  methods  described  in  the  previous  chapter  and  are  presented  in 
Table  I.  The  parameter  »M,  which  lumps  the  surface  states  and  the 

fixed  charges  in  the  oxide,  is  higher  than  is  common  for  commercial 
devices. 


A  typical  set  of  the  terminal  currents  measured  as  a  function  of 
drain  bias  is  shown  in  Fig.  (4.1).  The  increase  of  drain  current  and 
substrate  current  which  occurs  well  below  the  drain-substrate  diode 
reverse  breakdown  voltage  is  evidence  of  the  presence  of  avalanche 
multiplication.  The  sum  of  the  source  and  the  substrate  terminal  cur¬ 
rents  equals  the  drain  terminal  current  throughout  the  drain  bias 
range  investigated.  Furthermore,  the  source  terminal  current  remains 
insensitive  to  drain  bias  compared  to  drain  and  substrate  terminal 
currents.  These  facts  provide  firm  support  for  the  idea  that  only 
electrons  are  active  in  the  multiplication  process  and  that  the  secon¬ 
dary  holes  are  not  collected  at  the  source.  The  possibility  of  reach- 
through  can  be  eliminated  because  of  the  difference  in  the  dependence 
of  source  and  drain  terminal  currents  on  drain  voltage. 


F±g.  (4.2)  shows  a  plot  of  the  log  of  the  substrate  current  versus 
(Vd  ~  Vds)  '  Tt  shows  a  good  linear  relationship  in  agreement  with 
Eq.  (2.16).  This  supports  the  assumption  that  the  first  order  field 
approximation  obtained  from  the  one-dimensional  solution  of  Poisson's 
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equation  is  adequate  to  estimate  the  magnitude  of  electric  field  in  the 
CCD  region  for  the  study  of  avalanche  phenomenon.  It  Also  supports  ifefe 
idea  that  the  field  dependence  of  the  ionization  rate  in  Eq.  (2.10)  is 
a  good  approximation  even  when  surface  scattering  is  involved.  Eq.  (2.10) 
was  derived  for  the  case  where  optical  phonon  scattering  is  mainly 
responsible  for  reducing  the  carrier  energy  obtained  from  the  electric 
field.  However,  it  seems  to  hold  equally  well  when  surface  scattering 
is  present. 

The  value  of  the  average  multiplication  factor  M  was  determined 

from  the  ratio  of  the  experimental  values  of  the  terminal  currents  as 

given  in  Eq.(2.3).  The  electric  field  at  the  drain  junction,  E,,  was 

d 

determined  with  the  help  of  Eq.  (2.8)  using  appropriate  device  parameters. 

These  data  are  plotted  and  presented  in  Figs.  (4.3),  (4.4)  and  (4.5).  A 

semi-logarithmic  scale  is  used  so  that  these  plots  are  really  log1  (In  M) 
-1  10 

versus  E, 
d 

The  values  of  M  given  by  Eq.  (2.15)  were  compared  with  these  data. 

A  computer  program  was  used  to  determine  the  parameters  A  and  B  of 
Eq.  (2.15)  such  thht  it  will  give  least  sum  of  square  error  between  the 
two  sets  of  values  of  M.  The  points  marked  as  "theoretical"  are  the 
results  of  that  computer  program. 

Table  II  shows  the  values  of  the  parameters  A  and  B  as  determined 
from  the  procedure  described  above  for  both  n-  and  p-channel  devices. 
These  values  of  A  and  B  are  used  in  Eq.  (2.10)  to  calculate  the  ion¬ 
ization  rates  a(E)  which  are  shown  in  Fig.  (4.6). 

It  should  be  noted  that  the  values  of  A  and  B  obtained  here  differ 

considerably  from  those  presented  in  the  literature  for  avalanche  ion- 

1 6  17 

izarion  in  bulk  silicon.  ’  To  make  clear  the  nature  of  this 
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TABLE  II 
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u  m  ys; 


DEVICE  No.  18-L 


(l06cm/voIJ) 

t-o 

Figure  4.5 —  In  M  vs.  Ej  *  for  Device  18-L 
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difference,  Fig.  (4.6)  also  shows  a  plot  of  the  expression 

o(E)  =  1.6  x  106  exp  (-1.66  x  106/  E)  (4.1) 

which  is  accepted  to  be  valid  for  electrons  in  bulk  silicon. 

That  these  differences  exist  is  not  surprising  when  one  remembers 
that  surface  scattering  is  absent  in  the  case  of  bulk  silicon. 

Fig.  (4.7)  shows  a  plot  of  log  B  versus  log  for  the  devices  of 

Table  II.  This  plot  suggests  a  strong  correlation  between  the  two 
quantities.  The  slope  corresponding  to 

B  *  ^eff”1  (4.2) 

is  also  shown  for  comparison. 

This  result  is  reasonable  since  the  lower  the  effective  mobility, 

the  higher  the  electric  field  required  for  the  kinetic  energy  of 

electrons  to  exceed  the  threshold  energy  of  avalanche  ionization. 

Drain  current  noise  was  measured  at  three  different  drain  bias 

values.  Typical  spectra  are  shown  in  Figs.  (4.8)  and  (4.9).  The 

theoretical  limiting  value  of  the  thermal  noise,  given  by  (2/3)  4kT  gdQ, 

is  shown  in  the  lower  right  corner.  At  a  drain  voltage  of  5  volts 

there  was  negligible  substrate  current,  which  implied  the  absence  of 

the  avalanche  multiplication.  As  is  common  in  MOSFET’s  the  flat  portion 

of  the  spectrum  at  V  =  5  volts  is  above  the  limiting  value  of  I  = 

d  °  eqo 

(2/3)  4kT  g  .  The  excess  noise  may  be  due  to  one  or  more  of  the 
do 

following  causes: 

18  19 

a)  the  effect  of  ionized  bulk  impurities,  * 

b)  white  noise  of  non-thermal  origin  operating  in  the  channel, ^  or 

21 

c)  hot  electron  effect. 
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_  2 

In  the  upper  left  hand  corner  of  these  figures  the  factor  of  (M) 
is  shown.  M  was  determined  from  terminal  currents  as  described  earlier. 
It  is  seen  that  the  1/f  spectrum  at  different  drain  voltages  varies  by 

—  2 

the  factor  (M)  .  However  the  high  frequency  spectrum  tends  to  vary  as 
*sub  ^is  is  in  agreement  with  Eq.  (2.33)  which  is  presented  again 
below. 

Je,d  <f>  "  ^eqs  <f>  W2  +  JSub  S  (4.3) 


At  the  lower  frequency  end  of  the  spectrum  the  thermal  noise  com¬ 
ponent  of  I  ^  is  buried  underneath  the  flicker  noise  component,  and 
the  first  term  of  Eq.  (2.33)  predominates.  Therefore  the  drain  current 


noise  in  the  flicker  region  of  the  spectrum  can  be  approximated  as 

I  (M)2. 

eqs 


At  the  high  frequency  end  of  the  spectrum  I  g  consists 


mainly  of  the  thermal  noise  component.  If  high  drain  bias  makes  the 
substrate  current  larger  than  I  g  in  the  thermal  noise  region,  the 
second  term  predominates. 

These  points  appear  to  be  more  obvious  when  the  drain  current 

noise  at  different  frequencies  is  measured  as  a  function  of  drain  bias. 

Fig.  (A. 10)  shows  the  results  of  the  measurement  on  device  #  36-R.  In 

—  2  - 

this  figure  the  variation  of  (M)  and  of  Igub  M  with  drain  bias  are 

also  plotted  for  comparison.  The  low  frequency  noise  follows  the 
—  2 

shape  of  the  (M)  variation  whereas  the  high  frequency  noise  varies 
as  the  Isub  M  curve.  This  is  an  evidence  that  Eq.  (4.3)  properly 
describes  the  drain  current  noise  when  avalanche  ionization  is  present 
in  the  CDD  region  of  a  MOSFET. 

Further  evidence  of  the  validity  of  the  noise  model  is  presented 
by  the  date  of  Figs.  (4.11)  and  (4.12)  which  show  both  drain  noise  and 
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substrate  noise  currents  as  a  function  of  drain  bias.  At  high  drain  bias 

/o  teges  these  two  noise  currents_approach  the  same  value  of  I  M. 

Theoretical  value*?  j  2  .  sub 

d  sub  '  obtained  from  Eqs.  (2.33)  and (2. 34) 

are  shown  on  all  these  figures  to  emphasise  the  excellent  agreement 
with  experiment. 

The  noise  measurements  were  repeated  with  the  FET  in  common  sub¬ 
strate  connection  with  the  gate  a.c.  grounded  and  a  high-impedance  in 
the  source  lead.  Fig. (4.13)  shows  the  data  in  the  same  way  as  Figs. (4.11) 
and  (4.12).  Again  there  is  good  agreement  with  them.  As  discussed 
in  Chapter  III  this  is  experimental  proof  that  the  substrate  current 

noise  and  the  drain  current  noise  are  fully  correlated  in  avalanche 
multiplication. 

In  general  the  quantitative  model  derived  in  Chapter  II  describes 
the  avalanche  multiplication  noise  in  FET's  with  good  accuracy.  However 

there  are  certain  possibilities  which  might  cause  discrepancies.  These  ’ 
possibilities  are  listed  beiow. 

a)  At  high  current  densities  it  is  possible  that  a  certain  portion 
Of  the  secondary  holes  are  not  scattered  out  of  the  CDD  region  but  are 
collected  back  at  the  source  terminal.  The  rise  of  the  source  terminal 
current  at  the  far  end  of  the  drain  bias  in  Device  II  35-L  might  indicate 

this  possibility.  See  Fig.  (4.14).  Other  units  did  not  show  this. 

b>  As  the  CDD  region  is  driven  into  higher  fields  it  becomes 
possible  that  the  field  induced  localised  breakdown  around  chemical 
impurity  sites  and  local  defect  sites  will  contribute  to  extra  low 
frequency  noise.  We  have  no  conclusive  evidence  for  this  mechanism. 

The  accuracy  of  Eqs.  (2.33)  and  (2.34)  verified  experimentally  does, 
together  with  the  facts  established  from  the  d.c.  terminal  current 
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measurements,  prove  that  only  one  type  of  carrier  is  active  in  the 
avalanche  multiplication  process.  If  both  types  of  carriers  had  gone 
through  the  avalanche  multiplication,  the  noise  resulting  from  the 
random  fluctuation  of  avalanche  multiplication  should  have  been  much 
larger  than  the  magnitudes  predicted  by  these  equations. 
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V  CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  work,  a  detailed  study  of  the  avalanche  multiplication 
phenomena  in  the  CDD  region  of  an  MOSFET  has  been  carried  out.  The 
present  chapter  summarizes  the  important  results  obtained  and  discusses 
several  areas  which  deserve  further  study. 

Experiments  were  performed  on  both  n-channel  and  p-channel  silicon 
MOSFET* s  fabricated  for  this  study.  The  experimental  results  were  com¬ 
pared  with  the  theoretical  study  on  the  avalanche  current  and  asso¬ 
ciated  noise.  The  following  conclusions  are  drawn. 

1)  Experimental  evidence  is  established  for  the  presence  of 
avalanche  ionization  in  the  CDD  region  of  a  MOSFET,  when  the  drain  is 
biased  far  beyond  saturation. 

2)  Only  one  type  of  current  carrier  (electrons  in  n-channel 
devices  and  holes  in  p-channel  devices)  is  active  in  the  avalanche 
multiplication  process  in  FET's.  The  minority  carriers  generated  from 
the  avalanche  multiplication  are  removed  from  the  CDD  region  and  col¬ 
lected  at  the  substrate  terminal.  Both  the  terminal  current  measurement 
and  noise  measurement  substantiate  this  fact. 

3)  The  noise  magnitude  associated  with  the  avalanche  multiplication 
in  MOSFET  s  can  be  described  with  accuracy  by  the  following  two  equations 
derived. 
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(5.1) 


Icqd  "  Ieqs  ^)2  +  I8ub  M 
•*-eqsub  “  ^eqs  (^  “  1)  +  lsub  ^  (5.2) 


4)  The  first  order  approximation  for  the  electric  field  intensity 
in  the  CDD  region  obtained  from  the  one-dimensional  solution  of  Poisson's 
equation  provides  sufficiently  accurate  estimation  of  field  intensity 
for  the  study  of  avalanche  phenomena - 

5)  The  field  dependence  of  ionization  rates  in  MOSFET's  is 
accurately  described  by 

a(E)  =  A  exp  (— B/E) 

where  A  and  B  are  constants  for  the  following  range  of  electric  field 
investigated. 

electron:  1.0  x  10"*  volt/cm  to  3.0  x  10"*  volt/cm 

hole:  2.0  x  10"*  volt/cm  to  5.0  x  10"*  volt/cm 

The  value  of  parameters  A  and  B  determined  differs  somewhat  from  the 
bulk  silicon  case  and  shows  a  strong  correlation  with  surface  scattering. 

1)  Eqs.  (5.1)  and  (5.2)  should  also  be  valid  expressions  for 
the  avalanche  current  noise  in  Junction  FET's.  The  substrate  current 

and  noise  have  to  be  replaced  by  the  gate  current  and  noise,  respectively. 
These  equations  deserve  an  experimental  verification  for  Junction  FET 
case. 

2)  At  low  temperatures  the  avalanche  multiplication  phenomenon  is 
expected  to  be  more  appreciable.  The  extension  of  present  work  to  low 
temperatures  would  add  to  present  knowledge. 


3)  In  Chapter  IV  the  correlation  between  the  parameter  B  and 
the  effective  mobility  was  discussed.  A  further  study  leading  to  a 
quantitative  relationship  between  the  parameters  A,  B  and  other  device 
parameters  is  recommended.  As  a  matter  of  fact  if  this  quantitative 
relation  can  be  obtained,  the  avalanche  current  measurement  will  pro¬ 
vide  a  new  method  of  measuring  the  effective  mobility.  Previous 
measurements  of  effective  mobility  were  mostly  done  by  measuring  the 
channel  conductance. 

4)  Expecting  that  at  low  temperatures  or  in  a  short  channel 
device  the  avalanche  multiplication  noise  will  be  significant,  a 
study  of  a  circuit  which  will  minimize  this  noise  deserves  immediate 
attention.  As  discussed  in  the  chapter  on  experimental  method  the 
common  substrate  configuration  can  be  less  noisy  than  common  source 
configuration. 


APPENDIX  I 


CALCULATION  OF  THE  VARIANCE  OF  THE  MULTIPLICATION  FACTOR  M 

Assume  that  the  ionization  events  are  independent  events  whose 
occurance  is  equally  probable  at  any  point  along  the  path  of  length  d. 
Then  the  probability  of  having  m  ionization  evenps  along  the  path  d  is 
given  by  Poisson's  distribution; 

pm  =  u»  exp  (-u)/n!  (A-l.l) 

where  u  =  ad  is  the  average  number  of  ionization  events  and  ~ a  is  the 
average  ionization  rate  per  unit  length.  Obviously  the  probability  of 
having  no  ionization  events  along  any  fixed  length  Z  is  given  by 

pQ  (£)  -  exp  (-«£)  (A-1.2) 


The  probability  density  of  path  length  between  two  successive 
ionization  events  has  an  exponential  distribution  as  a  corollary  of 
Poisson's  law.  Hence,  the  probability  of  having  an  ionization  event 
at  a  point  between  y  and  y  +  dy  (assuming  the  last  ionization  event 
at  y=0)  can  be  written  as 


f(y)  dy  =  a  exp  (-ay)  dy 


(A-1.3) 


drain 

P(l)  = 


Defining  P(M)  as  the  probability  of  M  particles  emerging  at  the 
terminal,  we  have 


Po(d)  = 


exp  (-y) 


(A-1.4) 


Suppose  only  one  ionization  event  occurs  at  y  =  y'.  Refer  to 
Fig.  (A-l) .  Then  from  Eq.  (A-1.3)  the  probability  of  the  particular 
ionization  event  being  located  at  y'  is 


f(y')  =  a  exp  (-  ay') 


(A-l. 5) 


Therefore  the  probability  of  two  particles  emerging  at  the  drain  terminal 

due  to  one  ionizing  collision  event  located  at  y1  is 

dP(2)yt  =  f(y«)  PQ (d-y ' )  p0(d-y')  dy'  (A-1.6) 

(a)  (b)  (C) 

where  (a),  (b) ,  etc.  indicate  the  particular  paths  involved.  P(d-y')2 
accounts  for  the  fact,  that  two  particles  cover  a  distance  of  d-y'  without 
8°ieg  through  another  ionization.  The  ionization  event  which  we  assumed 
to  occur  at  y  can  happen  anywhere  along  the  path.  Thus  integrating 
over  the  total  distance  of  d,  we  get 
P(2)  =  /0  dP(2)y, 


=  /  a  exp  (-ay')  exp  [-  a(d-y')]2  dy' 


-  a  exp  (-2  a  d)  fQ  exp  (a  y')  dy’ 


=  exp  (-2u)  (exp  (u)  -l]  (A-l. 7) 

Now  in  order  to  have  three  particles  emerging  at  the  drain  terminal, 
two  ionization  event  should  occur  somewhere  along  the  path.  There  are 
two  possible  ways  in  which  this  can  happen,  as  indicated  in  Fig.  (A-2) . 
Case  I  is  when  the  primary  particle  goes  through  two  successive  ionizing 
collisions ;  Case  II  is  when  the  primary  and  the  first  secondary  particle 
each  go  through  one  ionizing  collision.  In  a  similar  way  as  in  the 
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Figure  A-l--Illustration  for  P(l)  and  P(2) 
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Figure  A-2--Illustration  for  P(3),  Case  I  and  Case  II 
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P(2)  calculation,  we  have  for  Case  I 


dP(3) 


y',  y"  -  f(y')  f(y"  -  y')  pQ(d  -  y*>  Po(d  -  y»)  Po(d  _  ytl) 

(a)  (b)  (c)  (d) 

For  Case  II,  we  have  the  same  form  of 


(A- 1.8) 


(e) 


dP(3) 


y',  y»  ~  f(y')  f(y"  -  y')  P0(d  -  y 1  )  Po(d  -  y")  Po(d  -  y»)  (A-1.9) 

(a)  Cb)  (c)  (d)  (e) 

Since  both  the  first  and  second  ionizing  collision  can  occur  anywhere 
along  the  path  as  far  as  they  occur  in  succession,  we  have  to  integrate 
by  dy'  and  dy"  over  the  entire  length  of  d. 

P(3)  =  dP (3)  dy'  dy" 

*  ^  f(y,)  P0(d  -  *">  -  y)  P0(d  -  y")2  dy")  dy’ 

J~Q  [a  exp  (-ay  )]  exp  [~  a(d  -y')]  /y,  a  exp  [-  a(y"  -  y«)] 

exp  [-  a^d  -  y")]2  dyu  dy. 

d 

o  6XP  (~2cd)  exP  (a  y')  [exp  (a  -  d)  -  exp  (a  -  y')j  dys 

=  1/2  exp  C-3y)  [exp  (u)  -  l]2  (A-1.10) 

Following  the  same  reasoning,  we  can  go  on  calculating  P(4) ,  P(5) 
and  so  on.  Summarizing  the  results  we  can  write  in  general 


P(M)  =  CM  exp  (-Mu)  [exp  (u)  -  1  ]M  “  1  (A-l.ll) 

where  C^  is  a  constant  number. 

The  sum  of  P(M)  for  all  the  possible  value  of  M  gives  the  normali¬ 
zation  of  the  probability  density  P(M) . 
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E  P(M)  =  1 
M=1 


(A-1.12) 


Differentiating  both  sides  of  Eq  (A-1.12)  with  respect  to  u  gives 


(A-1.13) 


Substituting  Eq.  (A-l.ll)  into  (A-1.13)  we  have 

xJt  1  =  PC0  -  1)  exp  (u)  [exp  (u)  -  l]"1  -  M) 

1  oil  JJ=:  1 

*=  E  P(M)  (-M)  +  exp  (u)  [exp  (u)  -  l]“* 

M=1 

E  P(M)  M  -  exp  (u)  [exp  (u)  -  l]“^  E  P(M)  =  0 
—  M=1  M=1 


Now  M  and  h  can  be  defined  as 


(A-1.14) 


M  =  E  P  (M)  M 
M=1 

OD  A 

M2  =  E  P  (M)  MZ 
M=1 


(A-1.15) 


Combining  Eq.  (A-1.14)  and  (A-1.15)  we  have 


M  =  exp  (u)  [exp  (u)  -  l]”1  M  -  exp  (u)  [exp  (u)  -  1] 


.-1 


=  exp  (u)  (A-1.16) 

Differentiating  both  sides  of  Eq.  (A-1.16)  again  with  respect  to  u,  we 
have 


—  M  =  ■—  [,  LP(M)  M  ] 

3u  3u  H=1 


=  E  M  —  P  (M) 
M=1 


=  exp  (u)  [exp  (u)  -  l]"1  E  P(M)  M  (M  -  1)  -  f  M2  P(M) 

M=1  M-l 


=  exp  (u) 


(A-1.17) 


Substituting  Eq.  (A-1.16)  into  the  above,  we  get 
M2  =  2  exp  (2u)  -  exp  (u) 

*  2  (M)2  -  M 
Hence, 

Var  M  =  M2  -  (M)2 

»  M  (M  -  1) 


(A-1.18) 


(A-1.19) 
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APPENDIX  II 


CALCULATION  OF  i  1  AND  i  i  * 
sub  d  sub 


terminal,  ngub,  can  be  written  as 


The  number  of  secondary  particles  emerging  out  of  the  substrate 
written  as 
nsub  =  nd'  “  ns 

ns 

”  ns 


nc 


=  Z  (Ms  —  1) 

i-1  '  (A**2. 1) 

where  nd  and  ng  are  the  number  of  particles  emerging  out  of  the  drain 

terminal  and  the  number  of  particles  entering  from  the  source  terminal 
respectively.  Applying  the  variance  theorem,  we  have 


^ar  nsub  =  “  l)2  Var  n„  +  n_  Var  M 


(A-2.2) 


Hence, 


Sisub  <0)  ’  q2  Var  "sub 


=  2q2  [ (M  -  l)2  Var  ns  +  ns  Var  M] 


(A-2.3) 


Expressing  the  above  equation  in  terms  of  equivalent  noise  diode 
current,  we  have 


Icqsub  =  IeqS  <W  -  l)2  +  Ifiub  M 

where  the  result  of  Eq.  (A-1.19)  is  also  used, 


(A-2.4) 


63 


For  the  correlated  part,  we  can  write 


Var  (n^  ngujj)  ~  nj  Hgu5,f  ~  nj  nsuj, 


(EM;)  [  E  (M;  -  1)]*  -  (n'  M)  [n~  (M  -  1)] 
i=l  j=l  J  s  s 


Hg  ng 

»  ill  jIxMi  CMj  "  1)*“  (ns)2  M  (M  -  1) 


»  ns  (ns  -  1)  M(M  -  1)  +  ns  (M2  -  M)  -  (ns)2  M(M  -  1) 
=  [ns2  -  (ns)2]  M(M  -  1)  +  ns  [M2  -  (M)2  ] 


=  Var  ns  M(M  -  1)  +  ns  Var  M  (A-2.5) 

Hence , 

Sidisub  (0)  =  q2  2  Var  (ndnsub) 

=  2q2  [M(M  -  1)  Var  ns  +  ns  Var  M  ]  (A-2.6) 

Again  expressing  the  above  in  terms  of  equivalent  saturated  noise  diode 
current,  we  have 


^eqdsub  *eqs  ^  +  *sub  ^ 


(A-2.7) 
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B.  A  NEW  GRATING  TYPE  Au-NSi  SCHOTTKY-BARRIER  PHOTODIODE  FOR  0.4-1.06 
MICRON  PHOTODETECTION  (S.S.  Li  and  C.T.  Wang) 


I.  Introduction 

The  qualities  desired  in  a  solid-state  photo-detector  are  high 


speed,  high  sensitivity  and  low  noise  operation  for  the  photon  ener¬ 
gies  of  interest.  For  years  it  has  been  known  that  a  Schottky-barrier 
(S.  B.)  photodiode  could  yield  these  qualities  if  incident  optical 
radiation  could  be  efficiently  coupled  into  the  depletion  layer  of  the 
diode  to  create  electron-hole  pairs  there.  Because  of  losses  associated 
with  coupling  light  through  the  metal  contact  into  the  semiconductor 
below,  however,  the  potential  of  an  S.  B.  photodiode  has  yet  to  be 
achieved . 

Schneider1  has  attempted  a  solution  to  the  problem^of  coupling 
optical  radiation  by  fabricating  contacts  of  thin  (~!00A),  semi¬ 
transparent  gold  film  coated  with  zinc  sulfide  as  an  antireflecting 
agent.  Photodiodes  made  in  this  way  suffer  two  main  disadvantages: 

(1)  they  require  complicated  processing  compared  with  that  used  in 


the  conventional  technology  of  silicon  devices;  and  (2)  the  detec¬ 
tor  shows  sensitivity  only  in  a  very  narrow  range  of  the  spectrum. 

In  this  paper  we  present  a  grating- type  Au-nSi  S.  B.  photo¬ 
diode  fabricated  with  a  new  structure  of  contact  mask.  The  mask 
forms  a  grating  of  gold  film  over  the  contacting  silicon  substrate. 
The  new  structure  can  be  readily  constructed  by  employing  the 
standard  silicon  planar  technology.  The  resulting  device  shows 
properties  (such  as  quantum  yield  and  spectral  response  et  al) 
much  superior  to  those  of  previous  S.  B.  photodiodes. 

In  the  following  sections,  a  brief  report  is  given  on  the 
I-V  characterisitics ,  spectral  response,  noise  performance  and  the 
technology  used  to  fabricate  such  grating  type  Au-nSi  S.  B.  photo¬ 


diode. 

II.  Device  Fabrication 

The  structure  and  geometry  of  the  grating  type  Au-nSi  S.  B. 
photodiode  is  displayed  in  Fig.  1.  The  device  was  constructed  by 
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(a)SIDE  VIEW  (c)  ENLARGED  PORTION 

OF  THE  WINDOW 
AREA  IN  (a). 


Au-nSi  SCHOTTKY  DIODE. 


Fig.  1 


Structure  and  geometry  of  a  grating  type 
Au-nSi  Schottky-barrier  photodiode. 


68 


first  growing  an  Si02  mask  on  a  1300  ohm-cm  n-type  silicon  wafer 
approximately  10  mil  thick.  A  window  for  metalization  was  opened 
on  the  oxide  mask  by  using  standard  photoelectric  and  etching  processes. 
A  thin  gold  film  was  deposited  onto  the  window  area  (which  is  the  bare 
silicon  surface)  and  extended  to  part  of  the  oxide  layer  (see  Fig. 

1(a))  by  vacuum  evaporation  technique.  The  thickness  of  the  gold 
film  is  not  crucial  in  this  case  and  may  be  of  a  few  microns.  The 
grating  pattern  of  the  gold  film  shown  in  Fig.  1(b)  was  made  by 
using  the  Shipley  process.  The  selection  of  grating  spacing  is 
determined  by  the  resistivity  of  silicon  substrate  and  the  device 
operating  bias  condition  (since  the  depletion  layer  width  depends  on 
both  the  doping  density  and  the  applied  reverse  bias  voltage) .  For 
example,  in  the  present  case,  with  a  1300-ohm-cm  silicon  wafer  being 
used,  the  grating  spacing  is  chosen  equal  12.5  ym,  which  is  suffi¬ 
cient  to  form  a  complete  depletion  layer  between  gold  gratings  and 

the  silicon  surface  at  zero  bias  condition  (see  Fig.  1(c)).  The 

-4  2 

total  area  of  the  device  is  25  X  10  cm  and  the  effective  area  for 

photodetection  is  9.3  X  10  ^  cm2.  However,  the  area  of  the  device 

could  be  either  enlarged  or  reduced,  depending  on  specific  applications. 

Note  that  extension  of  the  gold  film  into  part  of  the  oxide  layer  is 

to  serve  as  a  field  plate  so  that  soft  breakdown  could  be  avoided  un- 

3 

der  large  reverse  bias  condition. 

Ill,  Measurements  and  Discussions 

In  this  section,  we  report  the  results  of  our  measurements  on  the 
I-V  characteristics  of  the  device  under  reverse  and  forward  bias  condi¬ 
tions,  the  spectral  dependence  of  the  quantum  yield,  determination  of 
the  barrier  height  and  the  noise  characteristics. 

A.  I-V  Characteristic  curves: 

The  I-V  plot  for  the  device  G-6  and  G-9  under  forward  bi 3S  condi¬ 
tion  is  shown  in  Fig.  2,  for  temperatures  between  263°K  and  317°K. 

These  curves  can  be  described  by  the  Schockley  diode  equation  that 

qVF 

tf  "  Ts  exp  to’  (1) 

where  V  is  the  applied  voltage,  n  is  found  (from  Fig.  2)  to  vary  be- 

F 

tween  1.06  and  1.10  for  both  G-6  and  G-9,  and  I0  is  the  reverse  satura- 

4  b 

tion  current  which  is  given  by 

Ig  =  SA*  T2  exp(  -  -£~)  (2) 
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Fig.  2  The  I-V  characteristic  curves  for  devices 
G-6  and  G-9  under  forward  bias  condition 
and  for  different  temperatures.  I  is 
the  reverse  saturation  current.  Tfie 
curves  are  in  accord  with  the  prediction 

that  I  =  I  exp(— r^r),  with  n  varying 
from  1.06  tS  1.10. 


Here  A*  is  the  effective  Richardson  constant,  S  is  the  cross  section 

-3  2 

area  of  the  diode  (which  is  equal  to  1.57  X  10  cm  for  the  present 
case),  and  is  the  barrier  height. 

2 

The  barrier  height,  is  determined  from  the  slope  of  Ig/T 

versus  1/T  plot,  as  is  shown  in  Fig.  3.  The  result  yields  <£>  •=  0.8075eV, 

Bn 

which  is  in  good  agreement  with  the  reported  value  for  the  Au-nSi 
4 

system.  Note  that  values  of  Ig  used  in  Fig.  3  are  obtained  from  Fig. 

2.  The  effective  Richardson  constant  A*  is  calculated  from  Eq.  (2) 

and  is  found  to  be  84.2  for  T  =  300°K  and  for  e  <  5  X  103  V/cm. 

m 

Fig,.  4  shows  the  I-V  characteristic  cun/e  for  Device  G-6  under 
reverse  bias  condition  and  for  T  =  300°K. 

To  compare  the  result  with  theory,  we  use  the  following  expres¬ 
sion  for  the  reverse  leakage  current,  I  ,  which  is  given  by:3 

K 

IR  =  SA*TZexp(-  [1-exp  (-—)  ]  (3) 

where  V  is  the  applied  reverse  bias  voltage,  and  the  barrier  height, 
taking  into  account  the  image  force  lowering  effect  and  the  sta¬ 
tic  dipolar  charge  layer  effect,  can  be  written  as3 


V  ■  <*M  -  X)  -  (qem/4ireoes)’s  -  aeB  (4) 

where  <J>M  is  the  metal  work  function,  x  is  the  electron  affinity  of 
the  semiconductor,  a  is  a  parameter  to  be  determined  from  the  exper¬ 
iment,  and  em  is  the  maximum  electric  field  strength  in  the  depletion 
region  which  is  given  by^ 


€„  =[(2,ND/eoesXvD  +  vR)]** 


m 


(5) 


where  V_  is  the  diffusion  potential  and  N_.  is  the  donor  density. 

}  D  u 

/  The  second  term  on  the  right  hand  side  of  Eq.  (4)  represents  the 
^Inage  force  lowering  potential,^  and  the  third  term  denotes  the  lower- 
/  ing  potential  due  to  the  static  dipolar  charge  layer  presenting  at 

the  metal-semiconductor  interface.3  The  dipole  layer  is  thought  to  be 
a  fundamental  consequence  of  electronic  wave  function  penetrating 
from  the  metal  into  the  forbidden  gap  of  the  semiconductor.  This 
correction  to  the  barrier  height  is  assumed  to  be  a  linear  function 
of  electric  field. 

The  solid  line  in  Fig.  4  is  the  theoretical  curve  computed  by 
using  Eqs.  (3),  (4),  and  (5).  The  parameters  used  for  this  computation 
are  S  =  1.57  X  10-3  cm2,  A*  =  84.2,  T  =  300°K,  and  a  =1.74  X  10-5  cm. 
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3.3  3.4  35 

1000/T  (K  '  ) 

Fig.  3 

2 

The  reverse  saturation  current,  Is/T 
versus  1/T  for  devices  G-6  and  G-9. 
barrier  height,  <f>_  ,  deduced  from  the 
slope  of  this  plotnis  0.8075eV. 


THEORY 

EXPERIMENT 


The  I-V  characteristic  curve  for  Device  G-6  under  reverse  bias  condition  and  for  T  =  300°K. 
The  solid  line  is  computed  by  usin?  Eqs.  (2),  (3)  and  (4)  of  the  main  text.  The  parameters 
used  are  <J>  =  0.8075V,  A*  =  84.2  and  a  =  1.74  x  10-5  cm. 


Tt  is  noted  from  Fig.  4  that  the  results  are  in  excellent  agreement 
with  the  theory  for  VR  <  30V.  For  VR  >  30  V,  the  theoretical  pre¬ 
diction  exceeds  the  measured  values.  This  discrepancy  could  be 
corrected  if  a  different  value  of  a  is  used.  The  value  of  a  can  be 
determined  from  the  slope  of  £nIR  versus  (  VR  +  VR)  plot.  The 
barrier  height,  $Bn>  decreases  with  increasing  reverse  bias  voltage. 
For  example,  $Bn  =  0.8075v  at  VR  =  0  and  $Bn  =  0.7069V  at  VR  =  30V. 

B.  Spectral  Response  and  Quantum  Yield 

Fig.  5  shows  the  spectral  response  curves  for  devices  G-6  and 

G-9  for  V  =  0  and  10V,  and  for  T  =  300°K,  The  results  show  that 
R 

at  zero  bias  voltage,  the  net  quantum  yield  (i.e.  photo-current  per 
absorbed  photon)  exceeds  0.6  for  wavelength  between  0.5  ym  and 
0.85  ym  and  is  about  0.06  at  X  =  1.06  ym;  at  VR  =  10V,  the  quantum 
yield  varies  between  0.8  ~  0.85  for  X  between  0.55  ym  and  0.9  ym, 
and  is  about  0.15  at  X  =  1.06  ym.  This  result  can  be  explained  as 
follows.  The  present  device  is  ooerating  in  depletion-mode  photo¬ 
detection.  The  main  contribution  of  the  photocurrent  in  such  device 
comes  from  photocarriers  generated  in  the  depletion  region  and  the 
carriers  generated  within  one  diffusion  length  outside  the  depletion 
layer  region.  At  zero  bias  condition,  it  is  believed  the  depletion 
layer  of  the  present  device  may  not  well  cover  the  whole  silicon 
window  areas  between  the  gold  gratings ,  and  hence  the  quantum  yield 
may  not  be  very  high  due  to  the  recombination  losses  outside  the 
depletion  region.  However,  by  increasing  the  reverse  bias  voltage, 
the  depletion  layer  width  is  increased,  as  a  result  more  photo¬ 
carriers  are  generated  within  the  depletion  region,  and  hence  the 
quantum  yield  is  increased. 

The  responsivity  (i.e.  photocurrent/incident  photon  energy 
density/sec)  of  devices  G-6  and  G-9  is  also  calculated.  The  results 
yield  that  R  =  0.44  A/W  at  X  =  0.63  ym  and  VR  =  10V,  and  R  =  0.63 
A/W  at  X  =0.9  ym.  The  spectral  dependence  of  the  responsivity  for 
device  G-6  is  also  displayed  in  Fig.  5,  for  VR  =  10V. 

C.  Response  Time 

Since  the  present  device  is  operating  under  depletion  mode 
photodetection,  the  upper  limit  of  device  response  time  will  be  equal 
to  the  carrier  transit  time  across  the  depletion  layer  (  which  is  of 
the  order  of  lO-"^  sec.).  However,  since  the  RC  time  constant  of 
the  device  is  longer  that  the  carrier  transit  time,  the  response  time 
of  the  present  device  will  be  limited  by  the  RC  time  constant.  The 
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series  resistance  of  the  device  is  equal  to  1,5  KH  and  the  capacitance 
at  =  10V  is  found  to  be  3pf.  Therefore,  the  RC  time  constant  of 
the  present  device  at  V  =  10  V  is  4.5  nsec,  which  should  be  the  upper 
limit  of  the  device  response  time.  We  have  made  use  of  device  G-6 
to  detect  the  light  pulse  of  a  Strobotac  light  source  and  observed  the 
rise  time  of  light  pulse  to  be  0.15  ysec  and  pulse  width  of  0.4  ysec., 
in  good  agreement  with  the  value  specified  by  the  instructions  manual. 
This  indicates  that  our  detector  has  a  response  time  less  than  0.15 
ysec. 

D.  Noise  Measurements  and  Noise  Equivalent  Power  (NEP) 

The  uoise  measurements  on  our  detector  are  performed  by  employ¬ 
ing  a  low  noise,  high  input  impedance  (10MH)  preamplifier  incorporated 
with  a  noise  resistance  of  1.24  Kfl.  A  general  radio  (G.R.)  1381 
Random-Noise  generator,  a  G.R.  1925  multifilter  and  a  G.R.  1926 
multichannel  RMS  detector  have  been  used  for  noise  analysis. 

The  resulting  noise  spt'  .crum  is  plotted  in  Fig.  6,  for  device 
G-9,  with  the  reverse  bias  voltage  =  8.6V  and  the  primary  current, 

I  -  132  nA. 

P 

The  primary  current,  1^ ,  is  the  sum  of  the  dark  leakage  current, 

1^,  and  the  photocurrent,  I  The  broken  line  in  Fig.  6  represents 

the  full  shot  noise,  I„,  calculated  from  I„  =  /2ql  'Af . ^  The  measured 

N  N  p 

equivalent  noise  current  is  somewhat  lower  than  the  computed  full 
shot  noise.  This  is  expected  at  this  low  level  of  primary  current  due 
to  the  fact  that  effect  of  traps  in  the  interface  region  and  the 
recombination  centers  in  the  junction  region  may  not  be  completely 
negligible  .  If  this  is  the  case,  the  measured  equivalent  noise 

g 

current  is  usually  less  than  the  calculated  full  shot  noise. 

The  equivalent  excess  1/f  noise  current  for  G-9  (and  G-6) 

is  about  553  nA  at  *L0Hz  and  decreases  to  15nA  at  100  Hz.  It  has 

been  pointed  out  by  Hsu  that  the  1/f  noise  for  an  S.B,  diode  using 

3 

field  plate  structure  will  not  be  a  strong  function  of  the  reverse 
bias  voltage;  this  has  also  been  observed  for  the  present  case. 

The  noise  equivalent  power  (NEP)  is  defined  as  the  noise  equi¬ 
valent  current  per  responsivity  of  the  detector.  In  the  present 
device,  we  have  calculated  the  NEP  for  device  G-6  for  X  =  0.63  ym 
and  0.90  ym.  Table  1  summarizes  all  the  important  measured  and 
calculated  parameters  for  G-6  and  G-9  grating  type  Au-nSi  Schottky- 
Barrier  photodiodes. 
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TABLE  I 


Measured  and  Computer  Parameters  for  a  Grating-type 
Au-nSi  Schottky  Barrier  Photodiode, 

VR  =  10V  and  T  =  300°K 

/  0 

1.  Device  active  area:  9.3  X  10-  cin 

2.  Silicon  substrate  resistivity:  1300  f^-cm 

3.  Spacing  between  gold  gratings:  12,5  ym 

4.  Series  resistance:  1.5  Kf2  at  V_  =  IV 

F 

5.  Capacitance:  C  =  45  pf  at  V„  =  0,  C  =  3  pf  at  V„  =  15V 

R  R 

6.  full  shot  noise  current:  I„  =  2  X  10~14A  at  V„  =  IV 

N  R 

IN  =  4  X  10"UA  at  VR  =  10V 

7.  Photoresponse:  (at  V  =  10V  and  X  =  0.9  ym) 

K 

Quantum  yield:  0.85 
Responsitivity:  R  =  0.63  A/W 
NEP:  6.3  X  10~14  W 

Detectivity  (D*) :  0.48  X  1012  cm  -  Hz^/W 


REFERENCES 


1.  M.  V.  Schneider,  Bell  Sys.  Tech.  J..  45,  1611  (1966). 

2.  H.  Melchion,  M.  B.  Fisher  and  F.  R.  Arams,  Proc.  IEEE.  58 

1466  (1970).  - ’ 

3.  A.  Y.  C.  Yu  and  C.  A.  Mead,  Solid  State  Elec..  13,  97  (1970). 

4.  S.  M.  Sze,  Phys.  of  Semi-Devices,  p.  399,  Wiley  (1969). 

5.  J,  M.  Andre*7s  and  M.  P,  Lepselter,  Solid  State  Elec  13 

1011  (1970).  - - -  ’ 

6.  R.  J.  Minniti,  Jr.,  G.  W.  Neudeck  and  R.  M.  Anderson  J. 

Appl.  Phys..  42,  1886  (1971).  — 

7.  S,  T,  Hsu  ,  I.E.E.E.  Trans,  on  Electr.  Devices.  ED-17,  495 


79 


c.  QUANTUM  YIELD  OF  METAL- SEMICONDUCTOR  PHOTODIODES 
(S.  S.  Li,  F.  A.  Lindholm  and  C.  T.  Wang) 


I.  Introduction 

Gartner (1)  has  expressed  the  quantum  yield,  or  internal  quantum 
efficiency,  (defined  as  the  photon  current  produced  per  absorbed  photon) 
in  a  metal-n- type  semiconductor  photodiode  as 


0  =  1- 


(1+aL  ) 
P 


In  this  expression  a  is  the  absorption  coefficient  of  the  semiconductor, 
the  magnitude  of  which  depend,  strongly  on  the  photon  wavelength,  W  is  ’ 

the  depletion  layer  width,  and  Lp  is  the  diffusion  length  of  minority 
carriers . 

Eq.  (1)  predicts  that  the  quantum  yield  will  approach  unity  if 
<*W  »  1.  For  most  metal-semiconductor  photodiodes  this  inequality  is 
indeed  obeyed  in  the  short  wavelength  region  (near  UV  region)  of  the 
spectrum  where  the  absorption  coefficient  exceeds  104  cm"1. (2)  In  conflict 
with  the  prediction  of  Eq.  (1),  however,  experimental  observations (3)  show 
a  sharp  decrease  in  quantum  yield  with  decreasing  wavelength.  Although 
qualitative  explanations  have  been  suggested  to  account  for  this 
decrease,  no  quantitative  analysis  has  yet  been  given. 

In  this  paper  we  show  that  the  observed  fall-off  can  be  attributed 
to  the  presence  of  an  inversion  layer  near  the  metal-semiconductor  inter¬ 
face.  A  principal  result  of  the  paper  is  a  theoretical  expression 
including  the  effects  of  an  inversion  layer  and  displaying  the  functional 
dependence  of  quantum  yield  upon  doping  concentration  and  photon 
wavelength.  This  expression,  together  with  related  expressions  for  the 
width  of  the  inversion  region  and  for  the  electric  field  in  the  inversion 
region,  are  derived  in  the  Appendices.  In  the  presentation  immediately 
below,  we  apply  these  expressions  to  a  Au-nSi  Schottky  barrier  diode. 

The  theoretical  predictions  agree  well  with  observed  behavior,  including 
the  decrease  in  quantum  yield  seen  at  short  wavelengths. 

For  concreteness,  we  shall  assume  in  our  discussion  below  that  the  semi¬ 
conductor  is  n-type.  A  similar  treatment  will  describe  the  behavior  of  ^ 
p-type  Schottky  photodiodes. 
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11 •  Existence  of  an  Inversion  Layer  in  a  Metal-n-type  Semiconductor  Diode 

According  to  Shockley  and  Spenke, ^  an  inversion  layer  forms  in 
a  metal-n-type  semiconductor  diode  if  either: (1)  the  barrier  height  of 
the  diode  exceeds  one  half  of  the  forbidden  band  gap  of  the  constituent 
semiconductor,  or  (2)  a  large  concentration  of  acceptor  surface  states 
is  present  at  the  interface  between  metal  and  semiconductor.  Condition  (1) 
is  fulfilled  for  Schottky-barrier  photodiodes  such  as  Au-nSi,  PtSi-nSi, 
and  any  metal-n  GaAs.^"^  In  the  absence  of  donor  surface  states,  there¬ 
fore,  one  can  expect  a  p-type  inversion  region  to  exist  in  these  diodes. 

Fig.  1(a)  shows  the  energy  band  diagram  for  a  metal-n  semiconductor 
diode  when  an  inversion  layer  is  present  at  the  interface  between  metal 
and  semiconductor.  Fig.  1(b)  shows  a  schematic  representation  of  the 
carrier  profile  (not  to  scale)  in  the  inversion  region  and  in  the  bulk. 
Henish  finds  it  instructive  to  regard  the  inversion  layer  as  a  thin 
highly  conductive  film,  termed  a  channel  or  gutter,  bordering  the  surface. 
There  is  much  experimental  evidence  for  its  existence  in  germanium  and 
silicon. <7"9> 

111  •  Computation  of  the  Inversion  and  Depletion  Layer  Widths  for 
a  Au-nSi  Schottky  Diode 


As  is  derived  in  Appendix  A,  the  inversion  layer  width  x.  at  zero 
bias  voltage  is  given  by: 
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where  £q  is  the  free  space  permittivity,  £R  is  the  dielectric  constant 
of  the  semiconductor  (e.g.  £R  =  11.7  for  Si),  y  is  the  donor  density, 

VD  is  the  diffusion  potential,  ^  is  the  barrier  height  of  the  diode, 

Vg  '*'S  the  energy  gap  t*ie  semiconductor  expressed  in  electron  volts, 
vt  is  the  thermal  voltage  kT/q,  and  Ny  is  the  effective  density  of  states 
in  the  valence  band. 


For  a  Au  nSi  Schottky  photodiode,  the  logarithmic  term  in  the  bracket 


of  Eq. 


(2)  is  much  smaller  than  the  term 


Thus 
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(a)  BAND  DIAGRAM 


x 


(b)  CARRIER  PROFILE 


Fig.  1  (a)  Energy  band  diagram  for  a  metal-n-type  semiconductor  photo¬ 

diode  when  an  inversion  layer  is  present. 

(b)  Spatial  dependence  of  electrons  and  holes  in  the  bulk  of 

semiconductor  and  in  the  inversion  layer  region,  respectively. 
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Using  Eq.  (3),  we  compute  the  inversion  layer  width  for  a  Au-nSi  Schottky- 
barrier  photodiode  as  a  function  of  donor  density.  The  result  is  shown  in 

Fig.  2. 


The  total  depletion  layer  width,  W,  is  not  materially  affected  by 
the  presence  of  an  inversion  layer,  and  thus  the  conventional 
depletion  approximation  gives 


2£o£R(VD+Vf 

«ND  J 


(4) 


The  functional  dependence  of  W  for  a  Au-nSi  Schottky  barrier  photo¬ 
diode,  computed  by  using  Eq.  (4),  is  also  shown  in  Fig.  2;  VD  denotes  the 

reverse  bias  voltage. 

The  electric  field  strength  within  the  inversion  layer  is  given  by 
2qND  *  f  N  V (x)+V 

e(x)  -  [-(VV(X))  +  Vt(N“)exP(-  — . S)j  (5) 

The  derivation  appears  in  Appendix  A.  Eq.  (5)  predicts  that  the  electric 
field  holds  nearly  constant  in  the  inversion  layer,  as  shown  in  Table  AI. 
This  constancy  justifies  an  approximation  used  in  Appendix  B  to  derive 
the  expression  for  quantum  yield. 


IV.  Calculation  of  the  Quantum  Yield  for  a  Au-nSi  Schottky-Barrier 
Photodiode 

The  presence  of  an  inversion  layer  at  the  interface  between  metal 
and  semiconductor  has  an  important  effect  on  the  quantum  yield  of  a  metal- 
semiconductor  photodiode,  particularly  in  the  short  wavelength  region 
of  the  spectrum.  Derivation  of  the  quantum  yield  for  a  metal-n-type 
semiconductor  photodiode  in  which  an  inversion  layer  is  present  is  given 
in  Appendix  B,  and  the  result  is 
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where  e =  kT/qL^  is  called  the  critical  field,  e  is  the  mean  electric 

field  strength  in  the  inversion  region  deduced  from  Eq.  (5),  L  is  the 

n 

electron  diffusion  length  in  the  inversion  region,  and  L  is  the  hole 

P 

diffusion  length  in  the  n-region.  The  first  term  on  the  right  hand  side 
of  Eq.  (6)  is  due  to  the  contribution  from  the  inversion  region,  while 
the  second  term  represents  contributions  from  both  depletion  layer  region 
and  the  bulk  of  the  semiconductor.  Note  that  Eq.  (6)  reduces  to  that 
of  Gartner's  expression  (Eq.  (1))  if  x^  =  0. 
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Fig.  2  Dependence  of  the  inversion  layer  width  and  the  depletion 

layer  width  W  on  the  doping  densities,  Np,  for  a  Au-nSi  Schottky 
barrier  photodiode,  for  =  0  and  10V. 


We  now  apply  Eq.  (6)  to  calculate  the  quantum  yield  for  a  Au-nSi 

Schottky-barrier  photodiode.  For  a  Au-nSi  Schottky-barrier  photodiode, 

the^barrier  height  ^  is  0.81  eV(5>  and  the  band  gap  V  =  1.12  eV  at 

300%  is  expected  to  be  ^  >  Vg/2.  Thus,  because  <J>8>  V  /2,  one 

may  expect  an  inversion  layer  to  exist  in  the  absence  of  large 

concentrations  of  donor  surface  states.  In  computing  the  quantum  yield, 

0,  we  use  the  data  for  the  absorption  coefficient  of  silicon  given  by 

Dash  and  Newman  and  the  carrier  diffusion  length  deduced  from  the 

lifetime  data  given  by  Ross  and  Madigan. (12)  The  results  are  shown  1, 
Figures  3-5. 

Fig.  3  show,  the  quantum  yield  versus  photon  wavelength  for  a  Au-nSi 
photodiode,  with  doping  concentrations,  as  a  parameter  and  for 
VR  -  0.  The  result  indicates  that  the  quantum  yield  depends  strongly 
on  the  doping  densities  both  in  the  short  wavelength  and  the  long  wave¬ 
length  sides  of  the  spectrum.  In  the  long  wavelength  side,  the  quantum 
yield  increases  with  decreasing  doping  densities,  while  in  the  short 


wavelength  region  the  quantum  yield 
densities. 


increases  with  increasing  doping 


This  behavior  can  be  explained  by  using  the  results  shown  in  Fig.  2 
together  with  Eq.  (6)  for  the  quantum  yield.  From  Eq.  (6)  note  that 
the  depletion  layer  width,  W,  determines  the  quantum  yield  in  the  long 
wavelength  region;  the  quantum  yield  will  Increase  with  Increasing  w. 

From  Eq.  (4),  it  is  also  noted  that  the  depletion  layer  width  is  Inversely 
proportional  to  the  square  root  of  the  doping  density.  As  a  result, 
the  quantum  yield  is  expected  to  increase  with  decreasing  doping  density 
in  the  long  wavelength  region.  On  the  other  hand,  the  inversion  layer 
width,  x.,  controls  the  quantum  yield  in  the  short  wavelength  side  of 
the  spectrum;  the  quantum  yield  increases  with  decreasing  values  of  x  . 
Furthermore,  Eq.  (3)  shows  that  x.  is  inversely  proportional  to  the  1 
square  root  of  the  doping  concentration  tID.  Therefore  the  quantum  yield 
will  increase  with  increasing  doping  density  in  the  short  wavelength 

region  of  the  spectrum,  and  decrease  with  increasing  doping  densities  in 
the  long  wavelength  region. 
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QUANTUM  YIELD 


0.10 


0.01 :- 
0.2 


Fig,  3  The  quantum  yield,  n,  versus  photon  wavelength,  X,  for  a 
Au-nSi  photodiode  with  doping  density  as  a  parameter,  and 


for  V  =  OV.  The  result  shows  that  H  increases  with 

J  _  V _ I  -  1  -  J  -  1  •  .  I  .  1  4  4  .4 


increasing  doping  density  in  the  short  wavelength  region, 
and  decreases  with  increasing  doping  density  in  the  long 
waveleneth  reeion. 


"  \  l _ LLLZ 


To  compare  Che  cheoretical  curve  with  experience!  data,  we 
abricated  a  Au-nSi  Schottky-barrier  photodiode  having  a  gold  film 

— If*  thlCk  a"d  a  doplnS  density  of  the  n-Si  substrate 

yield*  d-CV  The  "eaSUred  “nd  the0retidal  «■!“«  for  the  quantum 

yreld  displayed  in  Fig.  4.  Note  that  the  experimental  data  in 

t  e  short  wavelength  regime  are  in  excellent  agreement  with  the 

theoretical  values,  while  in  the  long  wavelength  region  the  measured 

values  are  lower  than  the  calculated  values.  The  discrepancy  can  be 

a  tribute!  to  the  fact  that  the  reflection  coefficient  for  gold  film 

increases  in  the  long  wavelength  side  of  the  spectrum.1  To  lend 

credibility  to  this  explanation,  we  measured  the  quantum  yield  of  a 

grating-type  Au-nSi  Schottky  barrier  photodiode14  and  compared  with 

the  theoretical  curve  computed  from  Eq.  (6).  The  result  displayed  in 

ig.  5,  shows  that  the  theoretical  predictions  and  the  experimental 

data  are  in  excellent  agreement  over  the  whole  spectral  range. 

It  is  concluded  that  the  sharp  decrease  of  the  quantum  yield  in 

the  short  wavelength  region  for  a  metal-n-type  silicon  photodiode  can 

e  attributed  to  the  presence  of  an  inversion  layer.  By  increasing 

-  e  oping  density  of  the  semiconductor  one  can  reduce  this  effect. 
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QUANTUM  YIELD  (v) 


Fig.  4  Comparison  of  the  computed  quantum  yield  with  experimental 


data  for  a  thin  film  type  Au-nSi  photodiode  for  V  =  0  and 

14  -3  K 

Np  =  4x10  cm  .  The  theoretical  curve  has  been  normalized 


QUANTUM  YIELD  {v) 


Fig.  5  Comparison  of  the  computed  quantum  yield  with  experimental 
data  for  a  grating  type  Au-nSi  Schottky-barrier  photodiode, 
for  VR  =  10V  and  NQ  =  1013  cm  3.  The  grating  structure  is 
designed  to  overcome  the  reflection  loss  due  to  gold  film.14 
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Appendix  A:  Derivation  of  the  inversion  layer  width  for  a  metal-n-type 
semiconductor  diode 

If  we  neglect  the  majority  carriers  (electrons)  in  the  inversion 
layer  of  the  diode  shown  in  Fig.  1(b),  the  space  charge  carrier  density  is 
p(x)  =  q(ND  +  p (x) ), for  0  <  x  <  x±  (A.l) 


where 


-  qND,  xi  <  x  <  W 

<f>  (x)  V  +V(x) 

p(x)  =  Nv  exp  (-  -E - )=  Nv  exp  (-  - ) 

£  t 


(A.  2) 


In  writing  Eq.  (A. 2),  we  have  assumed  that  in  the  inversion  region  the 

quasi  Fermi  level  for  holes  is  nearly  constant.  This  assumption  is  valid 

for  light  intensities  small  enough  that  Ap  «  jp  where  AP  and  p^  denote 

o  o 

the  excess  density  and  the  equilibrium  density,  both  averaged  over  the 
width  of  the  inversion  layer.  Substituting  Eqs.  (A.l)  and  (A. 2)  into 
Poisson's  equation  gives 


^2~  =  ‘  (F^)LND  +  Nv  eXp(~  “V~)] 

dx  o  R  t 


(A. 3) 


The  electric  field  strength,  e(x),  within  the  inversion  layer  region  is 

obtained  by  integrating  (A. 3)  once  and  then  using  the  boundary  condition 
dVCx") 

that  — — —  =  0  at  x  =  W.  This  yields: 
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0  <  x  <  W 


(A.4) 


Eq.  (A.4)  is  used  to  estimate  the  electric  field  strength  within  the 
inversion  region  of  a  Au-nSi  Schottky-barrier  diode.  The  result  for 
different  doping  densities  is  computed  in  Table  A.l,  for  x  =  0  and  x  =  x.. 

For  0  <  x  <  x^,  note  that  the  electric  field  strength  within  the  inversion 
region  varies  by  at  most  25%  for  all  doping  densities. 

The  inversion  layer  width  is  obtained  by  integrating  Eq.  (A.4)  once 
and  using  the  boundary  conditions  that  V(x)  =  -  $  at  x  =  0,  and  V(x)  =  -V  /2 


for  x  =  x. 


The  result  is 
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(A. 5) 
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Eq.  (A. 5)  is  derived  for  the  ease  of  zero  bias  voltage,  and  can  be 
further  simplified  for  Au-nSl  Schottky-barrier  photodiodes.  The 
detailed  discussion  is  given  in  the  main  text. 

The  result  given  in  Eq.  (A. 5)  is  obtained  using  the  treatment 
employed  first  by  Bamburg(13)  and  Henisch. (16) 


ND(cm  3) 

io13 

io14 

io15 

io16 

io17 

£  /volt/cm. 
m^at  x=0 

1 . 20xl03 

3. 89xl03 

1. 29xl04 

4. 29xl04 

1.42xl05 

(at  x=x±) 

7.81xl02 

2. 79xl03 

9. 82xl03 

3. 37xl04 

1.15xl05 

Table  A. 1  Electric  field  strength  in  the  inversion 
region  of  Au-nSi  Schottky  barrier  photo¬ 
diodes  with  donor  density  as  a  parameter. 
Computed  from  Eq.  (A. 4). 
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Appendix  jB:  Derivation  of  the  onrnitiim  xHoia  e -  ...  , 

- - - — a  tor  a  metal-n-type  semiconductor 

photodiode  when  an  inversion  layer  is  prpspnt- 


The  current  density  Jph  produced  by  the  lncldant  (non.reflected) 

0  1.Q  fllO  eim  /vv«n  m>S  r  i  i 


|ju  -  -  - -  \ nun  l  i 

flux  «  is  the  superposition  of  hole  and  electron  components, 


*1  u  =  j  +  J  , 
ph  n  p  * 


which  applies  at  any  plane  x  in  the  diode.  Thus 


(B.l) 


where 


V  *  VV  +  W 

'  W  +  Jp<«>  +  Jdep 


Jdep=  W  '  yw> 


CB.2) 


(B.  3) 


IS  the  component  of  hole  current  density  arising  from  carrier  generation 
in  the  depletion  region.  To  find  the  functional  dependence  of  the  quantum 
y  eld  n,  therefore,  we  need  only  express  J  (x  ) ,  J  (W) ,  and  J  as 
functions  of  $.  nip  dep 

The  derivation  is  based  on  the  following  assumptions:  (1)  small 
njection  (An  or  Ap  «  nQ,  PQ),  (2)  small  reverse  bias  condition  (no 
avalanche  multiplication  occurs  within  the  depletion  region) . 

For  monochromatic  radiation,  the  rate  of  generation  of  photocarriers  is 


g(x)  =  $ae" 


(B.4) 


Jdep  =  ~q  f  g(x)dx  =  -q  /  $ae-ax  dx 

X.  x 

1  xi 

-ax. 

=  c<Ke  1  -  e'aW)  (B<5) 

where  x.^  is  the  inversion  layer  width  given  by  Eq.  (A. 5) 

In  the  bulk  of  the  n-type  semiconductor,  holes,  the  minority  carriers 
flow  mainly  by  diffusion.  Thus  the  continuity  equation  becomes 


Dp  T"!2  "  ^  +  s(x)  =o,  x  >  w 


P  9x2  TpTKW-w.  x>W  (B.6) 

where  Dp,  Tp  denote  the  hole  diffusion  constant  and  lifetime,  respectively. 
From  the  solution  '  of  this  equation,  one  finds  that 


J  (W)  =  (q$)  - - E -  - 

p  14  '  (1+aL  )  e 

p 


(B.  7) 


Since  the  inversion  region  is  very  close  to  the  metal  surface, 

Eq.  (B.4)  implies  that  the  excess  carriers  generated  in  this  region 
arise  mainly  from  the  absorption  of  short  wavelength  photons.  The 
pertinent  continuity  equation  is 

2 

„  3  An  i — i  9  An  An 

Dn  7 y  "  l^n  97”  "  ~  +  8(x)  =  °»  0  <  x  <  xi. 

ox  n 

Here  we  have  approximated  the  electric  field  in  the  inversion  layer 
by  the  mean  electric  field  7,  derived  from  Eq.  (A. 4). 

Substituting  Eq.  (B.4)  into  (B.8),  yields 


(B.8) 


32An 

3x2 


I  e  |  1  9 An  -ax 

L~  *7  -  -  <D~)e 
c  n  n 


(B.  9) 


kT 


where  ec  =  »  for  Ln  =  0.05  cm,  =  0.5  eV  at  T  =  300°K. 

For  the  case  of  Au-nSi  Schottky  diode, | e/ec j »  1  (see  Table  A.l). 
The  general  solution  of  Eq.  (B.9)  can  be  written  as: 


An  =  A  cosh  ( 


x^-x 


x .  -x  $L 


-ax 


)  +  B  sinh  (-f-)-(-JL) - _ - 

n  (aL  +|e/e  |) 


(B. 10) 


where  L  -  |— |  Lr  is  the  effective  electron  diffusion  length  in  the 
Inversion  region. 

Constants  A  and  B  in  hq.  (B.10)  are  determined  from  the  boundary 
conditions  that  An  =  0  at  x  =  xi  and  x  =  0.  The  results  yield: 


A  =  (-^) 


-ax, 


D  aL  +|e/e 
n  n  1  c 


(B. 11) 


and, 


B  =  (—*)  ( - ~ — 

n  aL  +| e/e 


-ax . 

1-e  1cosh(x./L) 
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I  sinh(x^/L)  ' 


(B. 12) 


n  1  c 1 

The  electron  diffusion  current  density  directed  towards  the  depletion 

region  may  be  evaluated  at  x  =  x,,  the  result  is  given  by: 

-Otx. 

dAn (x  )  q$aL  -ax  1-e  cosh(x,/L) 

J  (x. )  -  qD  — - -  =  ( - )[e  1  _  _ *  i  /R  ion 
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For  a  Au-nSi  Schottky  barrier  photodiode,  x±  «  L,  Eq.  (B.13)  reduces  to: 


Combining  Eqs.  (B.l)  through  (B.14),  the  total  photocurrent 
density  produced  in  the  diode  is  thus  given  by: 

q$0tL  -Otx  .  -  Ox .  -rvv 
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n  1  c 


~ax,-  i  ax  -ax . 

-)[e  bl  +  te 


-aW 

«-i5) 

P 


and,  the  quantum  yield  is  given  by: 


„  -  J)>  _  ,  aLn  .  ,  axi  .  i  -<«,  -axJ  „-aW 

n  "  7$  ~  C - Z - "He  -(-—-)  (1-e  1)1  +  [e  1  -  ® _ _wR1M 

q  aL  +  e/e  I  ax,  n+m.  v1  (B*16' 

n  c 


(1+aL^ 


Eq.  (B. 16)  will  reduce  to  the  expression  for  n  given  by  Gartner  (see  Eq.  (1» 
if  Xi  '  0  (no  inversion  layer).  Eq.  (B.16)  provides  a  theoretical 
description  of  the  fall-off  of  T)  seen  as  the  wavelength  of  the  incident 
photon  flux  decreases.  Qualitatively,  this  fall-off  occurs  for  the 
following  reason.  The  carriers  generated  by  photons  of  short  wavelength 
appear  mainly  near  the  surface,  in  the  inversion  layer,  as  is  implied  by 
Eq.  (B.4).  Thus  the  current  resulting  from  short  wavelength  excitation 
must  come  mainly  from  carriers  generated  in  the  inversion  layer.  Most  of 
the  holes  and  electrons  generated  there,  however,  ultimately  flow  toward 
the  surface  and  recombine.  The  charge  of  these  holes  and  electrons  thus 
cancel  and  contribute  no  net  current.  The  holes  are  directed  toward 
the  surface  by  the  electric  field  in  the  depletion  and  inversion  regions. 

The  electrons  are  directed  toward  the  surface  by  diffusion,  which 
dominates  the  drift  tendency  over  most  of  the  inversion  layer.  Only 
near  the  plane  x  =  xi  do  drift  and  diffusion  cooperate  to  direct  electrons 
into  the  depletion  region.  Thus  only  a  small  fraction  of  the  carriers 
generated  in  the  inversion  layer  contribute  to  the  photo-current  J 
We  have  termed  the  resulting  contribution  J^)  and  have  displayed^ 
functional  dependence  in  Eq.  (B.13). 


Eq.  (B.16)  has  been  used  to  calculate  the  quantum  yield  for  a 
Au-nSi  Schottky-barrier  photodiode,  for  different  bias  voltages  and 
doping  concentrations.  Note  that  the  present  derivation  does  not 
include  the  reflection  loss  at  the  metal  surface,  and  thus  the  actual 
quantum  yield  should  be  lower  than  the  predicted  values  by  Eq.  (B.16), 

by  a  fraction  of  R  (where  R  is  the  reflection  coefficient  at  the 
metal  surface). 
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D.  DETERMINATION  OF  ELECTRON  AND  HOLE  CAPTURE  RATES  IN  NICKEL-DOPED 
GERMANIUM  USING  PHOTOMAGNETOELECTRIC  AND  PHOTOCONDUCTIVE  METHODS 
(H.  F.  Tseng  and  S.  S.  Li) 


I .  Introduction 

It  is  known  that  a  substitutional  nickel  atom  introduces,  in 
agreement  with  the  tetrahedral  covalent  bond  model,  two  acceptor  levels 
in  the  forbidden  band  of  germanium.  One  df  them  is  located  at  0.30eV 
below  the  conduction  band'  edge,  and  the  other  is  0.22eV  above  the  valence 
band  edge;  both  levels  are  known  to  serve  as  efficient  recombination 
centers.  Nickel  atoms  can  be  introduced  into  germanium  crystal  by 
ordinary  diffusion  process. 

The  properties  of  nickel  doped  germanium  have  received  considerable 

attention  since  the  recombination  property  was  first  investigated  by 

Burton  et.  al.^  They  were  able  to  interpret  the  electron-hole  recombination 

2 

in  terms  of  the  first  acceptor  level.  Tyler  et.  al.  later  found  the 

second  acceptor  level,  but  failed  to  show  that  both  levels  were  associated 

3 

with  the  same  substitutional  nickel  atom.  Subsequently,  Battey  and  Baum, 

in  a  study  of  the  temperature  dependence  of  lifetime,  deduced  the  electron 

capture  rates  in  both  levels.  They  found  that  was  independent  of 

the  temperature,  while  increased  exponentially  with  temperature.  On 

the  other  hand,  Okada^  showed  that  photoconductivity  in  nickel-doped 

germanium  crystals  can  be  interpreted  in  terms  of  the  temperature 

independent  capture  rates  and  Kalashnikov  and  Tissen^  showed  that 

the  photomagnetoelectric  (PME)  lifetime  in  p-type  crystals  can  be 

interpreted  in  terms  of  a  two-level  model  assuming  that  and  C n2  are 

independent  of  temperature  and  with  ratio  of  Q^/Qn^  =  However, 

Wertheim,^*  in  his  study  of  the  bombardment-conductivity  decay-time 

experiments,  concluded  that  C  ^  and  C^2  a1-'2  independent  of  temperature 

but  with  the  ratio  C  ,/C  „  =  1/6,  which  is  contradictory  to  the  result 

ml  nl  r 

of  Kalashnikov  and  Tissen.  This  discrepancy  was  later  explained  by 
Eliseev  and  Kalashnikov7  as  due  to  dislocations  in  their  early  germanium 


crystals,  and  concluded  that  C^2  is  larger  than  Cnl  with  the  ratio 
Cnl/Cn2  =  1/6* 

Note  that  the  ratios  of  electron  capture  rates,  Cnl/Cn2 
f)  7 

by  Wertheim  and  Eliseev  et.al.  are  consistent  at  room  temperature 
(i.e.,  Cnl/Cn2  =  ~).  However,  the  values  of  Cnl  and  Cn2  reported  by 
Eliseev  et.al.  are  a  factor  of  2  smaller  than  that  of  Wertheim  (see 


Table  I).  The  values  cf  hole  capture  rates  reported  by  previous  authors 


are  even  more  widely  spread  and  less  conclusive.  The  main  reason  for 
this  is  due  to  the  fact  that  most  of  these  results  are  determined  from 
the  temperature  dependence  of  lifetimes,  which  is  further  complicated 
and  prone  to  lead  incorrect  results  if  some  of  the  recombination 
parameters  (e.g.,  capture  rates)  become  temperature  dependent. 

In  this  paper,  we  shall  demonstrate  that  by  deluding  the  excess 
carrier  trapping  effect  in  the  theory  of  photomagnetoelectric  (PME)  and 
photoconductive  (PC)  effects  and  by  elaborately  controlling  the 
compensation  ratio  of  nickel  and  shallow  impurity  densities  in  germanium, 
we  are  able  to  determine  carrier  capture  rates  more  accurately  at  room 
temperature  from  PME  and  PC  measurements.  Furthermore,  by  closely 
examining  the  previous  results,  we  are  able  to  interpret  consistently 
the  present  results  with  those  reported  by  Wertheim  and  Eliseev  et.al. 

II.  Theory 

The  energy  levels  of  nickel  atom  in  germanium  and  the  transitions 
among  bands  and  the  nickel  levels  are  shown  schematically  in  Fig.  1. 

The  nickel  atom  in  the  three  possible  charge  states  (.with  zero,  one 
and  two  electrons  in  the  two  acceptor  levels  of  nickel  atom)  are 
designated  by  Nq,  and  N2>  which  are  further  related  by  the  total¬ 
ity  condition 

No  +  N1  +  N2  =  (1) 

where  is  the  total  concentration  of  nickel  atoms.  The  parameters 

nl»  pi»  n2  and  p2  rePresent  the  carrier  densities  when  the  Fermi  level 
coincides  with  acceptor  levels  E^  and  E^ ,  respectively.  They  are 
related  to  the  densities  of  impurity  charge  states  through  the  fol¬ 
lowing  relations^: 

°1  =  Vo/pl 


o 


(3) 


where  we  have  used  the  superscript  "o"  to  denote  the  equilibrium 
densities  and  is  the  Fermi  level.  The  above  equations  also  pro¬ 
vide  the  relationship  of  the  equilibrium  nickel  densities  in  different 
charge  states  to  the  Fermi  level.  From  Eqs.  (1),  (2)  and  (3),  the 
density  of  nickel  atoms  in  each  charge  states  can  be  determined  from 
the  Hall  effect  and  the  conductivity  experiments. 

From  Fig.  1,  the  rate  equations  in  the  bands  and  the  nickel  ac¬ 
ceptor  levels  can  be  written  as 


The  charge  neutrality  condition  expressed  in  terms  of  the  excess 
densities  from  the  equilibrium  values  is  given  by: 


Ap  -  An  -  Aiv^  -  2  AN2  =  0 


Under  the  condition  of  small  injection^  (i.e,,  An«  n  +  n,. 

°  1 


Ap  P0  +  pi  ♦  J  “  2),  and  the  steady  state  case,  the  lifetimes  for 

J  o 

electrons  and  holes  are  defined  respectively  by  : 

2 

T  =  An/  j  R  . 
n  m 


T  -  ip/  l  R 


FroraEqs.  (4),  (5.),  (6)  and  (7),  the  trapping  ratio,  T,  and  the 
lifetimes  can  be  expressed  by^ 

r  =  An  =  in  _  1  "  glg2  +  (1  +g2^Jpl  +  ^1+[V1Jp2 


Ap  Tp  1  -  (3^2  +  (1  +^2^ni  + 


where 


"n'1  '  <po+r'V  'CnlVSo0+Nl°>  /H1  +  Cn2Cp2(N1°+N'20> 

'  <e„+r"V  f  !'cPoi  ("o+ni)  +  T„oi  W1"' 

+  l\o2  (no  +  n2)  +  T„o2  <po  +  p2)rlt  (9) 

H1  =  Cnl  (no  +  V  +  Cpl  (pD  +  Pl}  <10> 

H2  =  Cn2  (no  +  n2}  +  Cp2  (Po  +  P2}  Cll) 

ynl  “  CnlNo  ^Hl*  pn2  =  Cn2Nl  ^H2  (12) 

yPl  =  CplNi°/Hl>  ^p2  =  Cp?-N2°/K2  C*3) 

31  =  No0/  <No  +  Nl°>  (14) 

32  -  N°/  (N°  +  n2°)  (15) 

Tnol  =  1/Cnl  (No°  +  N?  *  Tno2  "  ^2  (Ni°+  N2>  <“> 

Vl  ■  1/Cpl  <N0  +  Kr>  ■  Tpo2  ■  1/Cp2  <N1°+  N2)  (17) 


For  steady  state  case,  =  Rp^  =  R^,  R^2  ■  Rp2  =  R2’  anc*  t^le 
relative  importance  of  the  net  steady  state  recombination  rates  of  the 

g 

two  levels  can  be  obtained  from  the  following  expression  : 

^9 

Fj  =  C„2.  (no  +  PlCpl/Cn.l.)/Cpl(po  +  C„2n2/Cp2)  <18) 

For  n-type  cample,  nQ  »  Pq,  P^/C^  »  Cn2n2/C  2,  therefore 

R  >>  R. ;  the  upper  level  dominates  for  all  temperatures,  and  the  re- 
1  1  XI 

combination  process  reduces  to  the  simple  Shockley-Read  model, 
reduces  to  Eq.  (A7)  of  reference  (11). 


Eq.  (9)  then 


For  p-type  samples, 

VR.  -  — —  — i 

Cnl  Po  (19) 

The  relative  importance  of  the  two  acceptor  levels  depends  on  the 
ratio  of  C^/C^  and  the  resistivity  of  the  sample. 

For  the  case  of  small  injection,  by  taking  into  account  the  trap¬ 
ping  effect,  the  steady  state  PME  open  circuit  voltage  can  be  expressed 
in  terms  of  the  photoconductance  ’  : 


n  I’-1  +  p 


V,  ,  o_ _ 

PME  Mp  (b  +  r1)  nob  +  P0 


o 


1  D  1 

(£) 


n  o 

where  VpME  is  the  PME  open  circuit  voltage  per  unit  length,  0 


(20) 


9P  =  ypHB  are  Hal1  angles  for  electrons  and  holes,  respectively.  AG 
is  the  photoconductance  per  unit  length  to  width  ratio  of  the  sample, 
and  Gq  is  the  dark  conductance.  Other  parameters  not  defined  here 
have  the  conventional  meanings.  By  meanings  of  Eqs.  (8),  (9)  and 
(20) ,  we  can  compute  the  capture  rates  and  lifetimes  from  the  PME 


and  PC  data. 

HI.  Experimental  Details 


N-type  germanium  single  crystals  were  cut  into  rectangular  bars 
with  typical  dimensions  of  0.8  X  5  X  12mm3.  The  shallow  donor  dens¬ 
ities  were  determined  from  the  Hall  effect  and  the  conductivity  mea¬ 
surements.  These  samples  were  cleaned  by  chemical  process  and  coated 
with  nickel  by  vacuum  evaporation.  The  samples  were  then  placed  in 
the  furnace  for  diffusion.  At  the  end  of  the  diffusion  period,  the 
samples  were  quenched  to  room  temperature.  The  nickel  concentrations 
were  controlled  by  the  diffusion  temperature.  The  nickel-doped  sample 
was  mounted  in  a  AC-310  liquid  helium  refrigeration  system  for 
measurements.  The  Hall  effect  and  the  conductivity  measurements  were 
taken  by  standard  d.c.  technique.  The  PME  and  PC  measurements  were 
taken  by  using  a  chopped  light  source  and  an  a-c  detecting  system, 
so  that  the  temperature  gradient  effect  can  be  eliminated. 

The  samples  for  investigation  were  divided  into  two  groups: 
one  group  has  the  nickel  densities  in  the  range  that  <  2N.. 

In  this  way,  the  nickel  doped  germanium  samples  remained  as  n-type, 
but  with  high  resistivity.  The  second  group  was  prepared  with 
N4  >  ND’  s0  that  samPles  became  p-type  after  nickel  diffusion.  The 
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nickel  densities  can  be  determined  accurately  from  Hall  effect  and  the 
conductivity  measurements.  We  shall  next  discuss  the  results  of  our 
measurements  on  these  two  groups  of  samples  separately. 

N-type  samples 

For  this  group  of  samples,  the  first  acceptor  level  (i.e., 
level)  of  the  nickel  atom  is  fully  occupied  by  electrons,  and  the  se¬ 
cond  acceptor  level  (i.e.,  E2)  is  partially  occupied.  The  occupation 
probability  of  the  second  level  and  the  densities  of  nickel  in  each 
charge  states  can  be  determined  from  Eq.  (3).  In  order  to  show  that 
nickel  is  the  only  dominant  impurity  in  these  samples,  we  measured 
the  Hall  coefficient,  electrical  resistivity  and  electron  densities  as 
a  function  of  temperature  before  the  PME  and  PC  measurements. 

Fig.  2  shows  the  temperature  dependence  of  Hall  mobility  for  tem- 

15  -3 

perature  between  180°k  and  300°k,  for  samples  Ge  -  N2  (N^31  3.5  X  10  cm  , 

N-  =  2  X  1015cm-3).  The  Hall  mobility  is  found  nearly  independent  of 

'L 

temperature  in  this  temperature  range.  The  electron  drift  mobility 

13 

u  can  be  deduced  from  the  relation  that 
n 


Here  y^H  and  are  the  low  field  Hall  mobility  and  Hall  coefficient 
respectively,  and  R  is  the  Hall  coefficient  at  very  high  field.  For 
electron,  this  ratio  is  nearly  equal  to  unity  in  the  temperature  range 
of  interest^13*  1^.  Fig.  3  and  Fig.  4  illustrate  the  temperature 
dependence  of  resistivity  and  electron  concentration,  respectively, 
for  sample  Ge  -  N2.  Since  the  electron  mobility  is  nearly  independent 
of  temperature,  the  activation  energies  deduced  from  the  slope  of  these 
two  curves  are  identical  and  equal  to  0.30  eV.  This  represents  the 
activation  energy  of  the  second  acceptor  level  of  nickel.  (i.e., 

E2  =  0.30  eV  below  the  conduction  band  edge.) 

The  PME  and  PC  measurements  were  taken  at  room  temperature,  and 
the  result  is  displayed  in  Fig.  5,  for  B  =  3  kG. 

The  linear  relationship  between  Vp^  and  AG/Go,  as  predicted  by 
Eq.  (20),  is  observed  in  Fig.  5. 

To  calculate  the  recombination  parameters  from  Eqs.  (8),  (9)  and 
(2)  ,  we  examined  over  ten  samples  with  different  donor  and  nickel 
concentrations.  The  parameters  n^,  p^,  n2  and  P2  used  in  computations 
have  the  following  values  at  300°k.: 
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o  Ge-N2 


A  Ge-Pl 


«r  SLOPE  =  0.22eV 


10  /T  (  K) 


Figure  4 

Carrier  densities  versus  temperature  for  samples  Ge-N2  and  Ge-Pl.  The 
activation  energy  for  Ge-N2  is  found  equal  0.30  eV  which  is  the  second 
acceptor  level  of  nickel  atom  in  Ge.  The  activation  energy  for  Ge-Pl  is 
found  equal  0.22  eV  which  is  the  first  acceptor  level  of  nickel  atom  in  Ge. 


mm 


nj  =  3.0  X  1011  cm"3 
P-L  =  1-2  X  1015  cm"3 
n2  =  9.7  X  1013  cm"3 
P2  =  3.8  X  1012  cm"3 


The  electron  diffusion  coefficient  Dn  is  obtained  from  ^  through  the 

use  of  Einstein’s  relation  that  D  =  rl m  The  ratio  of  electron 

n  J  Mn‘ 

and  hole  mobilities  is  taken  from  the  results  reported  by  Prince14. 

For  n-type  samples,  the  upper  level  is  dominant  and  the  lower 
level  has  little  effect  on  the  trapping  ratio  T  (since  p  is  very 
large),  Eqs.  (8)  and  (9)  can  be  simplified  into  the  forms  that: 


and 


r  = 


no  +  n2  +  y2(po  +  P2  +  N2> 

no  +  n2  H'  Y2(Pq  +  ?2)  +  Nx 


(21) 


T  _1  =  (Po  +  r\>'[po2 


n 


Cno  +  n2>  +  Wo  +  P2 D  <22> 


where 


Y2  Cp2/Cn2 ’  and  \,o2  and  Tpo2  are  given  b?  E9S-  (16  —  1  7). 


From  Eqs.  (20)  and  (22),  it  is  noted  that  VpME  and  x^1  depend  greatly 
on  the  trapping  ratio  T  which  in  turn  depends  on  the  ratio  y  . 

From  the  present  experimental  data,  we  can  calculate  y As  well 
33  ^n2  and  %2*  result  yields; 


Cn2  "  ®  X  ^  cm3/sec 
Ep2  =  4^  ^  cm3/sec 
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For  comparison,  we  summarize  our  results  together  with  the  results 
reported  by  other  investigators  in  Table  1. 


For  this  group  of  samples,  the  first  acceptor  level  (i.e.,  E^)  are 
partially  occupied,  and  the  second  acceptor  level  are  completely  empty. 
The  results  are  shown  in  Fig.  2  to  Fig.  5  for  sample  Ge  -  PI  (with 
=  1.4  X  lO^cm  =  7.5  X  lO^cm  ^) .  Fig.  2  shows  the  hole 

Hall  mobility  versus  temperature.  The  temperature  dependence  of  the 

hole  Hall  mobility  is  found  to  be  T  ~  which  is  in  good  agree- 

—1  6  1  A 

ment  with  the  T  *  dependence  of  hole  mobility  reported  by  Prince1  , 

The  temperature  dependence  of  the  resistivity  shown  in  Fig.  3  yields 
a  slope  of  0.20  eV  which  is  less  than  the  activation  energy,  =  0.22eV, 
of  the  first  acceptor  level.  This  is  due  to  the  fact  that  hole  mobility 
increases  with  decreasing  temperature.  Fig.  4  shows  the  hole  concen¬ 
tration  versus  temperature,  and  the  slope  of  this  plot  yields  the 
activation  energy  of  the  first  level,  which  is  found  to  be  0.22eV 
above  the  valence  band  edge.  The  hole  drift'  mobility  is  deduced  from 
the  Hall  mobility  and  the  ratio  R^/  which  ranges  from  1.4  to  1.8 

**  1  <3  i  / 

depending  on  the  temperature  and  impurity  concentration.  -  The 
linear  relationship  between  Vnx_,  and  AG/G  is  also  observed  for  p- 
type  samples,  as  is  ; hown  in  Fig.  5  and  predicted  by  Eq.(20). 

For  p-type  samples,  pQ  »  nQ,  also  p^  »  n^,  from  Eqs.  (9)  and 
(20),  it  is  noted  that  the  trapping  ratio  T  has  little  effect  on  the 


values  of  Vp^  and  T^,  also  can  be  ignored  in  the  expression  of 

T  .  Thus,  for  p-type  samples,  we  were  unable  to  determine  C  ,  and  y, 
n  Pi  1 

(where  y^  =  Cp^/C^)  accurately.  The  relative  importance  of  the 

recombination  in  both  levels  for  this  case  depends  not  only  on  the 


resistivity  of  the  samples,  but  also  on  the  ratio  of  C  „/C  , .  By 

n2  nl 

varying  the  resistivity  of  the  samples,  which  is  achieved  by  changing 


the  nickel  doping  densities,  we  can  obtain  the  ratio  of  C  _/C 

nz  nl 

From  this  result,  we  conclude  that  C  „  is  larger  than  C  .  with  the 

n<£  nl 

ratio  of  cn2/cni  =  fi. 

For  samples  with  pQ  <  p^,  the  second  acceptor  level  (i.e.,  E2) 

dominates  the  recombination  process,  and  t  is  nearly  equal  to  T 

n  ,  no2  , 

In  this  case,  the  value  of  obtained  from  the  p-type  sample  agrees 
well  with  that  obtained  from  n-type  samples . 
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,-y.  rj-  .  , 

For  samples  with  PQ  <  p  ,  the  effect  of  the  first  acceptor  level 
(i.e.,  E1)  becomes  important.  The  value  of  C  x  can  be  determined 
from  the  PME  and  PC  measurements .  The  result  yields 

;  _  =  1.4  x  10  ^  cm^-sec  \ 

O' 

For  comparison,  the  above  result  is  also  included  in  Table  1. 

IV.  Discussions  and  Comparison  with  Previous  Work 
,r,  In  computing  carrier  capture  rates  from  lifetime  measurements  for 
the  nickel-doped  germanium  crystal,  the  reasons  for  the  wide  spred  in 
the  previous  results  may  be  attributed  to  the  following  facts:  (1) 
experimental  difficulties  in  determining  accurately  the  density  of 
electrically  active  nickel  atoms;  (2)  failure  to  include  the  excess 
carrier  trapping  effect  in  the  PME  theory;  (3)  improper  assumptions  on 
the  temperature  dependence  of  hole  capture  rates,  C^;  and  (4)  neglecting 
term  due  to  hole  capture  rates,  Cp2  in  the  lifetime  expression  for  n-type 
samples,  which  results  in  underestimating  the  electron  capture  rate,  C 
at  room  temperature. 

In  an  attempt  to  avoid  all  the  possible  errors  cited  above,  we  have 
employed  a  unique  approach  in  determining  the  carrier  capture  rates  for 
nickel-doped  germanium  crystals  at  room  temperature.  First,  we  have 
elaborately  controlled  the  compensation  ratio  of  nickel  atoms  and  the 
shallow  impurity  densities  (see  Section  III)  such  that  the  density  of 
nickel  atoms  can  be  measured  accurately.  Secondly,  instead  of  computing 
capture  rates  from  temperature  dependence  of  lifetime  data,  we  have 
determined  the  electron  and  hole  capture  rates  at  one  single  temperature 
(i.e.,  room  temperature)  in  both  n-  and  p-type  samples;  this  has  the 
advantage  of  eliminating  possible  errors  arising  from  (3)  discussed 
above.  Thirdly,  we  have  included  the  excess  carrier  trapping  effect  in 
our  PME  and  PC  theory;  our  result  has  shown  that  the  effect  of  trapping 
is  more  prominent  in  n-type  crystals  than  in  the  p-type  samples. 
Implications  of  the  above  criteria  to  our  present  work  as  well  as  previous 
work  are  discussed  as  follows. 

The  ratio  of  electron  capture  rates,  C^/C^,  obtained  by  us  is  1/6 
at  room  temperature,  which  is  in  excellent  agreement  with  that  reported 
by  Wertheim  and  Eliseev7  et.  al.,  respectively.  However,  the  values  of 
Cnl  and  Cn2  obtained  her®  are  somewhat  higher  than  that  of  Wertheim  and 
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Eliseev  s,  while  value  of  Cp2  obtained  here  is  much  lower  than  that  of 

Wertheim  but  agree  well  with  Burton  et.  al.1  This  discrepancy  can  be 

explained  by  the  following  facts.  Wertheim  deduced  the  ratio  of  C  /C 
-  nl  n2 

from  the  two  lifetime  plateaus,  one  at  low  temperature  (i. e.  C  )  and 

nl 

the  other  near  room  temperature  (i.e.,  C^).  Although  C  is  independent 

of  temperature  due  to  its  neutral-center  nature,  it  has  been  reported  by 

Belgaev  and  Miselguk  that  C^2  decreases  with  decreasing  temperature 

for  140°K  <  T  <  280°K  due  to  Coulomb  repulsion,  and  is  essentially 

independent  of  temperature  for  280°K  <  T  <  350°  due  to  the  tunnelling 

penetration  of  the. Coulomb  barrier  by  electrons.  The  room  temperature 

value  of  Cp2  reported  by  Wertheim  was  taken  from  his  low  temperature 

lifetime  data  by  assuming  that  Cp2  is  independent  of  temperature.  There 

has  been  no  temperature  dependent  data  on  Cp2  reported  as  yet.  However, 

it  is  expected  that  value  of  C  _  will  decrease  with  increasing  temperature 

i'a 

due  to  its  attractive-center  nature.  Thus  the  value  of  C  2  reported 

by  Wertheim  at  room  temperature,  was  overestimated,  and  C  _  should  be  less 
-9  3  _i  -  p2 

than  100  x  10  cm  -sec  at  T  =  300°K.  In  computing  C  „  at  room 

nz 

temperature  from  lifetime  data  for  n-type  crystals,  Wertheim  employed 

Eq.  (1)  of  his  paper  and  neglected  term  due  to  Cp2.  This  results  in 

underestimating  the  value  of  C^,  because  Cp2  is  not  negligible  at  room 

temperature.  In  fact,  we  have  substituted  values  of  C  _  and  C  obtained 

n2  p2 

in  this  paper,  into  Eq.  (1)  of  his  paper,  and  found  that  lifetimes  computed 
from  this  result  are  in  excellent  agreement  with  his  measured  lifetimes 
for  samples  532  at  room  temperature.  Therefore,  the  value  of  C  2  reported 
by  Wertheim  is  somewhat  underestimated,  while  the  value  of  C  „  deduced 

p2 

from  his  low  temperature  data  is  overestimated  at  room  temperature. 

Our  present  results  support  such  conclusions. 


Next,  we  examine  the  values  of  Cnl  and  C^2  reported  by  Eliseev  and 

Kalashnikov.  Their  analysis  of  the  lifetime  data  is  essentially  similar 

to  that  of  Wertheim,  however,  they  obtained  values  of  C  _  and  C  „  which 

nl  n2 

are  a  factor  of  two  lower  than  that  of  Wertheim.  This  is  due  to  the  fact 
that  they  determined  C^  and  C^2  from  the  PME  method,  which  they  failed 
to  include  the  trapping  effect  in  their  PME  theory.  They  pointed  out 
that  the  lifetime  data  observed  from  photoconductivity-decay  experiment 
was  only  about  50  ~  60%  of  their  PME  lifetimes.  If  we  use  the  lifetime 
data  they  obtained  from  photo-decay  experiment  to  analyze  Eliseev's  data, 
the  results  for  and  C^2  will  be  close  to  the  values  of  ours.  Thus, 
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the  values  of  C  ^  and  reported  by  Eliseev  et.  al.  are  underestimated, 
which  is  due  to  the  fact  that  they  neglected  the  excess  carrier  trapping 
effect  in  the  PME  theory. 

Finally,  the  results  of  capture  rates  reported  by  Burton  et.  al.3 

were  obtained  from  lifetime  data,  which  they  determined  from  diffusion 

length  measured  over  a  wide  range  of  conductivity  at  room  temperature. 

Although  the  value  of  hole  capture  rate  reported  by  them  agrees  quite 

well  with  ours,  it  is  noted  that  this  data  was  available  before  the 

existence  of  the  upper  level  had  been  recognized.  The  present  result 

confirms  that  the  hole  capture  rate  due  to  upper  level  C  „  is  equal  to 
-9  3-1 

40  x  10  cm  -sec  at  room  temperature. 

V.  Conclusions 

From  the  present  study  of  the  capture  rates  for  nickel-doped 

germanium  crystal  determined  by  PME  and  PC  methods  at  T  =  300°K,  we 

conclude  that:  (1)  The  ratio  of  C  ^/C  2  equal  to  1/6,  which  is  in 

good  agreement  with  that  of  Wertheim  and  Eliseev  et.  al.;  (2)  Our  result 

yields  C^  =  1.4  x  10  ^  cm3-sec  3,  and  C J  2  =  8  x  10  ^  cm3-sec  3;  this 

result  is  consistent  with  Wertheim *s  result  if  C  „  was  included  in  his 

p2 

Eq.  (1)  in  computing  C^  from  room  temperature  lifetime  data,  and  agrees 
well  with  Eliseev  et.  al.'s  result  if  the  trapping  effect  was  included 
in  their  PME  theory.  (3)  Our  result  for  C^2  (40  x  10-9  cm3-sec_1)  is 
in  excellent  agreement  with  Burton  et.  al.'s  value,  although  their  study 
was  completed  before  the  discovery  of  the  upper  level;  the  present  result 
is  also  consistent  with  Wertheim 's  if  his  low  temperature  data  on  C^2 
was  corrected  to  the  room  temperature  value  without  approximation. 

In  short,  we  have  demonstrated  that  by  including  the  trapping  effect  in 
the  PME  and  PC  theory  and  by  elaborately  controlling  the  compensation 
ratios  of  nickel  atoms  and  shallow  impurity  densities,  we  are  able  to 
determine  more  accurately  than  previous  results  three  of  the  four 
capture  rates  for  nickel  atoms  in  germanium  at  room  temperature  from 
PME  and  PC  measurements. 
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E.  AMORPHOUS  SEMICONDUCTOR  FILMS 


(P.  K.  Chaudhari,  E.  R.  Chenette,  and  A.  van  der  Ziel) 


This  report  is  concerned  with  the  results  of  an  intensive  study 
of  the  electrical  and  optical  properties  of  amorphous  films  of  a 
typical  chal cogen ide  semiconductor.  The  report  is  divided  into  two 
almost  independent  parts  for  convenience. 


PART  I:  OPTICAL  PROPERTIES  OF  As2Se3 * Sb2$e3  FILMS. 


A.  INTRODUCTION 

The  band  structure  and  the  transport  properties  of  chalcogenide 

semiconductor  materials  have  attracted  increasing  interest  in  tecent 

12 

years,  especially,  since  new  electronic  devices  ’  have  been  reported. 

The  use  of  these  materials  as  photodetectors  is  of  great  practical 
interest.  The  understanding  of  the  carrier  transport  behavior  in  these 
material  is  still  in  its  infant  state.  Many  of  the  properties 
connected  with  the  electric  conductivity  seem  to  vary  widely  from 
sample  to  sample.  The  reasons  are  not  fully  understood.  Investigation 
of  optical  absorption  and  photoconductivity  has  provided  valuable  insight 
in  the  physics  of  crystalline  semiconductors.  The  main  contribution  of 
work  reported  here  is  to  apply  these  methods  to  vitreous  semiconductors 
to  obtain  information  about  electronic  structure.  The  emphasis  is  on 
the  optical  properties  deduced  from  optical  transmittance  data  on 
evaporated  films  of  As2  Se3  *Sb2  Se3  at  300°K.  These  data  are  analyzed 


to  determine  optical  constants  (refractive  index,  extinctinction 
coefficient  k  or  absorption  coefficient  a  =  %k/A) .  An  analysis  of  the 
absorption  edge  has  also  been  carried  out  to  obtain  the  optical  band-gap, 
and  the  nature  of  the  transitions  involved  (direct  or  indirect).  An 
interpretation  and  discussion  of  the  results  that  are  related  to  the 
energy  band  diagram  and  electron  transport  behaviour  is  presented 
in  the  discussion  of  electrical  properties  (Part  II)  which  follows. 
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B.  EXPERIMENTAL  TECHNIQUE 


The  compounds  As^e^and  Sb2Se3  used  were  of  high  purity 
(semiconductor  grade  99.999+  %) .  Each  compound  was  weighed  separately 
with  precision  and  then  mixed  together  and  ground  to  a  fine  powder  in  a 
mortar.  The  resulting  material  was  mixed  thoroughly  in  a  steel 
cylinder  with  a  steel  ball  by  vibrating  for  three  minutes  to  obtain 
a  homogeneous  mixture.  Disks  of  I  cm  diameter  and  about  1mm  thick 
were  formed  in  a  press  at  7000  lb/ in. 2  at  room  temperature  under  vacuum. 
These  alike  were  then  used  as  sources  for  evaporation.  The  evaporation 
was  carried  out  in  a  beel  jar  system  at  a  pressure  of  about  10-^ 
torr  from  a  quartz  crucible.  The  substrates  were  at  room  temperature. 

A  shutter  arrangement  permitted  deposition  to  occur  only  after  the 
desired  constant  evaporation  rate  had  been  established.  Ten  separate 
runs,  under  identical  conditions,  were  made.  The  films  were  deposited 
on  glass,  and  on  NaC£  crys.als.  ihe  glass  substrates  were  degreased 
thoroughly  and  washed  with  deionized  water.  They  were  then  placed 
in  an  ultrasonic  cleaner  containing  hot  trichloroethylene  (TCE)  for 
15  minutes,  and  then  rinsed  with  TCE  followed  by  acetone  and  stored  in 
methyl  alcohol.  The  NaC£  crystals  were  polished  with  running  deionized 
water  and  rinsed  with  alcohol  and  also  stored  in  methyl  alcohol. 

The  film  thickness  was  determined  by  a  Sloan  Angstrommeter  with  an  accuracy 
of  +  100  A.  For  films  greater  than  3  y  the  interferrometr ic  thickness 
was  used  to  calculate  the  optical  constants  of  the  amorphous  films. 

X  -  ray  diffraction  on  two  films  indicated  a  vitreous  structure. 


Transmittance  was  obtained  using  10  different  samples  at  room  temperature, 
from  2 . 5“ 1 6  y  using  a  Beckman  model  IR-10  spectrophotometer.  All 
the  data  were  taken  under  identical  conditions  with  a  reproducibility 
better  than  \/Q.  Optical  absorption  was  measured  at  room  tempeature 
and  at  liquid  nitrogen  (LN2)  temperature  from  0.7-2  y  using  a  Carry 
model  l*i  spectrophotometer.  For  LN2  measurements  a  specially  constructed 
cryostat  having  a  brass  finger  for  sample  mounting  and  with  appropriate 
quartz  windows  for  light  transmittance,  was  placed  in  the  sample 
compartment  of  the  Carry  U.  Absorbance  curves  were  measured  for  6 


d iff rent  samples  at  LN^  temperature  and  for  10  different  samples  at  room 


temperature,  one  after  the  other,  under  identical  conditions  of 
Carry  14  settings.  Because  there  might  be  a  difference  in  composition 
from  run  to  run  due  to  the  Variations  in  substrate  temperature  and  evaporation 
rate,  which  were  not  controlled  precisely  during  this  study,  the  absorption 
coefficient  wa*  calculated  for  individual  samples. 


C.  THEORIES 


Analysis  of  the  optical  data  can  best  be  approached  by  first 
summarizing  the  various  theoretical  results  and  their  relationships  to 
the  physical  phenomenon  of  absorption.  These  results  will  be  employed 
wherever  appl icable. 

If  no  interference  fringes  are  observed,  the  optical  density  D  *  of  a  f 
is  given  by^ 


D  =  ad/2.3  (1) 

Here  a  is  the  absorption  constant  and  d  is  the  film  thickness.  This 
relation  assumes  double  reflection  at  the  film-air  and  substrate-film 
interfaces  and  an  exponential  law  for  the  absorption  process.  The 
absorption  coefficient  a  is  related  to  the  extinction  coefficient 
k  by 

a  =  4itk/X  *  1.01  X  105  kE  cm"1  (2) 


where  A  and  E  -  hu  are  the  free  space  wavelength  and  photon  energy  in 
electron  volts,  respectively.  Thus,  k  can  be  found  simply  from  the 
measured  value  of  D 
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If  a  is  small  enough  to  allow  the  light  wave  to  traverse  the  sample 
several  times,  interference  fringes  result  in  the  transmittance  response 
■  The  values  of  the  wavelengths  at  which  T  shows  a  maxima  and  minima 
can  be  used  advantageously  to  determine  the  refractive  index  n.  The 
condition  for  the  occurrence  of  extremum  values  for  T  is4 


4nd  =  mA  1  - 


[(4af  *.  4fs)  /raJ 


(3) 


(maxima  for  even  m;  minima  for  odd  m).  Here,  m  is  an  interger. 

The  parameters  4>f  and  «fj  are  phase  angles  due  to  reflections  at 
the  air-film  (af)  and  film-substrate  (fs)  interfaces  and  are  given  by 


'af 


=  tan 


1  £  2k/ ^n2  -  1  +  k^yj 


C4a) 


’fs 


=  tan 


1 


nl  +  *2 


)] 


(4b) 


where  n,  is  the  refractive  index  of  the  nonabsorbing  (k  =  0)  substrate. 

If  k  is  small,  the  S's  may  be  neglected.  In  our  films,  for  n  =  2  26 
k  i  10  J  and 
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nl  -  1.50  which  are  typical  numbers;  the  values  of  6^  and  6^  are  negligibly 
small  (<  10  ^)  .  Thus  the  simplified  form  of  Eq.  (3), 


4nd  =  mX 


(5) 


may  safely  be  used  to  determine  n,  provided  m  is  known.  If  the 
value  of  n  does  not  change  very  rapidly  with  the  wavelength  as 
was  found  to  be  true  in  the  present  study,  then  by  assuming  ^nd 

-  (m  +  l)Xj  and  ^nd  =  (m  +  1)^  for  two  successive  interference  maxima, 

we  find  that 

"  "  V(h  -  *2)  (6) 

Hence  Equation  (5)  gives  n  at  the  various  wavelengths. 

The  classical  theory  oT  absorption  and  dispersion  of  a  free  carrier 
gas  in  the  amporphous  state  leads  to  the  following  expressions  for 
the  optical  constants:”* 


2  ,  2 

n  -  k  =  e 


1  -  T2/(l  *  w2t2V| 


(7) 


and 


2nk  = 


t/ 


(8) 
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Here  q  is  the  electronic  charge,  em  is  the  dielectric  constant  arising 
from  the  combined  effects  of  the  bound  charges,  T  is  the  carrier 
density,  and  m*  is  the  effective  mass  of  carriers.  These  equations 
may  be  simplified  for  the  limiting  cases  of  high  and  low  frequencies. 
In  the  first  case  (to  T  »  1),  and  n  »k  ,  Eq.  (7)  reduces  to 


n  =  e 


«  e 


m 

^Kq2/m*e^  ^1/u2) 


(9a) 


(9b) 


2  2 

Here  top  =  Nq  /mAe^is  the  plasma  frequency. 


2  _2 


For  low  frequencies  (u  r  «  1)  and  Eq.  (7)  reduces  to 


(10) 


2 

Hence  n  does  not  depend  any  more  on  X  (if  T  is  independent  of  the  energy). 
T  is  related  to  (he  free  carrier  mobility  by  the  expression 
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When  the  energy  of  the  incident  photon  exceeds  the  energy  E 

(optical  band  gap)  an  electron  can  be  directly  excited  from 

the  valence  band  into  the  conduction  band.  In  the  case  where 

the  valence  band  maximum  and  the  conduction  band  minimum  appear  at 

the  same  point  in  the  Brilloin  zone  (namely  at  =  0)  the  transitions 

from  the  oneband  to  the  other  are  referred  to  as  direct  transitions 

between  parabolic  energy  bands.  The  energy  dependence  of  the  absorption 

coefficient  a  near  the  edge  of  this  fundamental  absorption  is  then  given 

by 


T  *  m*y/q 


(11) 


This  relation,  when  plotted  as  a  versus  hv,  will  yield  a  straight  line 
giving  a  value  of  E  t  by  extrapolation  to  a  =  0. 
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0.  RESULTS  AND  DISCUSSION 


A  typical  transmission  curve  in  a  region  extending  from 
2.5-16;i  is  shown  in  Fig.  1.  The  transmission  of  the 
film  continues  to  show  interference  fringes  up  to  16  y, 
indicating  that  plasma  resonance  is  absent  in  amorphous 
films.  This  is  expected  because,  due  to  the  highly 
disordered  lattice,  the  time  between  collisions  is 
considerably  lowered  so  that  uit  «  1  should 
characterize  such  a  damped  system.  As  we  shall  see 
later,  dispersion  behavior  of  the  refractive  index  at 
long  wavelengths  supports  this  conclusion. 

Interference  maxima  in  the  transmittance  curves  were 

used  to  obtain  the  optical  constants.  The  dispersion 

of  the  refractive  index  is  illustrated  in  Fig.  2.  The 

curve  shows  a  relatively  sharp  increase  towards  the 

short  wavelength  region  and  approaches  a  constant  value 

of  2.26  at  iong  wavelengths,  so  that  the  relative 

2 

dielectric  constant  e  =  n‘  =  5,11.  This  value  is 
typical  of  the  chalcogenide  glass.  The  constancy  of 
the  refractive  index  at  long  wavelengths  is 
characteristic  of  a  highly  damped  system  (u>T  «  1). 

Wavelength  dependence  of  n  may  also  be  described  by  the  limiting  cases 

of  the  classical  Eq.  (7)  for  the  short  wave  side  and  the  long  wave 

side  as  shown  in  Fig.  3.  It  is  clear  that  n  is  independent  of  X 

in  the  long  wave  region,  in  good  agreement  with  Equations  (9)  and 

(10)  respectively.  The  transistion  between  these  two  regions  occurs  at 

1  k 

Xx  =  3.08  y  or  0)T  =  6.1*»  X  10  .  According  to  the  classical  theory 

this  frequenct  should  be  equal  to  the  reciprocal  mean  free  time  1/t 
of  the  free  carriers.  Therefore  t  =  1.63  X  lO’1^  sec.  The  intercept 
at  X  =  0  also  gives  em  =12.  The  magnitude  of  the  mobility  is 
determined  by  using  Eq.  (11)  and  is  given  by 
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WAVELENGTH  X  (microns) 


Fig.  2  Refractive  index  vs  X  for  sample  P 


p 


qx/m* 


2.S6  cm^/Vsec 


(13) 


where  m*  =  m. 

Figure  shows  complete  curves  of  the  dependence  of  a  on 

the  photon  energy  obtained  from  the  absorbance  data  below  the 

interference  effect  regime,  taken  at  room  temperature  and  at  Ll^ 

temperature.  (Note  that  the  energy  range  covered  here  is  the  same 

as  that  explored  by  Platakis  et  al,  ^  on  this  composition  in  bulk 

form.)  All  the  samples  essentially  follow  the  same  law.  Note  that  there 

is  no  observable  exponential  tail.  This  agrees  with  the  data  reported 

by  Platakrs  et  al.,  on  the  same  composition  in  bulk  samples,  and 

by  Donovan  and  Spicer  ^  on  amorphous  germanium  films.  The  observed 
2 

a  ^  hv  dependence  indicates  allowed  direct  (momentum  conserved) 
transitions  between  the  conduction  and  valence  bands.  When  extrapolated 
to  a  *  0,  the  optical  -^ap,  averaged  over  10  samples,  is 
EQpt  =  ].kk2  eV  at  at  room  temperature  and  Eopt  =  1.558  eV  averaged 

over  6  samples  at  LN_  temperature.  Assuming  that  the  absorption  edge 

^  8 

depends  linearly  on  temperature  (such  a  linearity  has  been  assumed 

for  As0  Se_),  the  thermal  coefficient  of  the  absorption  edge  is 

3  =-  5-3  X  10  eV/°K,  which  agrees  with  the  value  found  by 

g 

Kolomiets  and  Manontova  for  As2Te^  glass. 

In  order  to  understand  the  observed  optical  behavior  it  is 
necessary  to  have  a  suitable  band  model.  A  plausible  band  structure 
consistent  with  both  the  optical  and  electrical  data  is  proposed  in 
Part  II.  For  the  present  purpose  it  suffices  to  briefly  describe  the 
model  consistent  with  the  optical  data.  Unfortunately,  there  is  no 
quantitative  calculation  available  regarding  the  band  structure  of  amorphous 
materials.  Band  models  proposed  so  far  are  based  on  speculative  assumptions 
that  the  band  structure  of  the  energy  spectrum  is  essentially  retained^’^ 
and  is  not  determined  by  the  long  range  order,  buy  is  dependent  upon 
the  short  range  order.  The  avsence  of  long  range  order  influences, 
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however,  the  detailed  shape  of  the  band  structure  in  that  the  band 

edges  are  not  well  defined,  and  there  are  present  in  the  gap  a  large 

number  of  localized  levels  due  to  the  random  potential  fluctuations  at 
12 

the  atomic  sites.  Localization  presents  low  mobility  because  of  trapping 
effects.  Consequently,  the  region  between  the  limits  of  the  localized 


levels  in  the  conduction  and  valence  bands  acts  as  a  pseudogap, 

1 3 

usually  defined  as  the  mobility  gap.  Thus  by  assuming  a  localization 


range  equal  to  1.50  eV  (obtained  from  the  electrical  data  in  Part  II), 
2 

a  ^  hv  dependence,  that  yields  E  ^  =  l . 442  eV,  can  be  explained  as 
follows. 


The  optical  absorption  coefficient  a  is  given 


o 


NC(E)  N'V(E  +  hv)  dE 


04) 


Here  P.f  is  the  transition  probability,  and  N  (E)  and  N  (E  +  hv) 
are  the  densities  of  states  in  the  conduction  and  valence  bands  respectively. 
For  transitions  beyond  the  mobility  gap  i.e.,  involving  the  nonlocal ized 
states  only  and  based  on  the  assumption  that  the  range  of  the  localized 
states  is  small,  Eq.  (l4)  is  used  to  interpret  the  optical  absorption  data 
on  amorphous  germanium1**  and  selenium  films1"’’1^.  Optical  transitions 
between  two  localized  states  are  not  likely  to  contribute  to  absorption 
because  they  a'e  spatially  separated.  Transitions  between  occupied 
localized  states  in  the  valence  band  edge  and  nonlocal ized  empty  states 
in  the  conduction  band  across  the  mobility  gap  edge  could  be  important 
however,  and  these  have  not  been  considered  by  any  of  the  above  authors. 

In  fact,  the  observed  optical  data,  on  amorphous  As2Se2‘Sb2Se^ 
films  investigated  here,  can  indeed  be  interpreted  in  terms  of  these 
transitions.  Assume  that  these  transitions  contribute  to  the 
absorption.  It  follows  that  a)  Nv(E)  ^  E1^,  and  that  b)  Nc ( E ) 
is  large  near  the  mobility  edge  and  can  be  taken  as  constant. 

Furthermore,  the  transition  probability  P.^,  as  observed  from  the  data 
(see  Fig.  k.)  is  constant. 


Substituting  these  results  in  Eq.  (14)  shows  that 


In  the  case  of  d i rect-a i lowed-trans i t ions ,  the  same  dependence 
is  found  by  Smith17  for  crystalline  semiconductors. 

Equation  (15)  represents  the  data  observed  in  Fig. 4.  This 
in  good  agreement  with  the  proposed  band  diagram  in  Part  II. 
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SUMMARY  AND  CONCLUSION 


The  results  of  the  optical  investigations  of 
sb2Se3  evaPorate<^  amorphous  films  may  be  summarized  as  follows: 


As2Se^  + 


5. 


The  real  part  of  the  complex  dielectric  constant  shows  a  few 
features  that  have  been  predicted  by  the  classical  theory  of  absorption 


and  dispersion  by  free  carriers;  the  square  of  the  refractive 


index  n  is  proportional  to  X  in  the  short  wavelength  region. 


2.  The  transistion  between  these  two  regions  occurs  at  the 


frequency  a>x  which  is  equal  to  the  reciprocal  mean  free  time 

~15 

sec;  was 


1/x  of  the  free  carriers.  The  value  of  x  =  1.63  X  10 
determined 


3.  The  mobility  calculated  from  the  value  of  x  is  <  5  cm2/Vsec. 


The  absence  of  plasma  resonance  observed  in  the  long  wavelength 


region  enabled  us  to  calculate  the  value  of  the  staic  dielectric 
2 

constant,  e  =  n  =5,11. 


The  optical  absorption  was  interpreted  as  being  related  to 
direct  transitions  involving  the  localized  states  in  the 
conduction  band.  The  optical  energy  gap  extrapolated  to  a  =  0  is 


E  .  =  1.^2  eV. 
opt 


6.  The  thermal  coefficient  of  energy  was  found  to  be  dE  /dT  = 
-L  opt 

-  5.3  X  10  4  eV/°K. 
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PART  II:  ELECTRICAL  PROPERTIES  OF  As^ * Sb^  FILMS 


A.  INTRODUCTION 


During  recent  years,  a  physical  model  of  the  electrical 
structure  of  amorphous  materials  has  evolved  largely  through  the 
work  of  Mott,^  ^  Cohen,  Fritzsche,  and  Ovshinsky^  and  Cohen. 

Although  the  general  features  appear  to  be  in  agreement  with  experimental 
evidence,  there  still  is  a  serious  lack  of  experimental  information 
related  to  the  electronic  transport  properties.  This  part  of  the  report 
is  concerned  with  the  results  of  measurements  of  electrical  conductivity 
as  a  function  of  temperature,  applied  field,  and  illumination.  The 
results  of  these  studies,  together  with  those  of  Part  1,  are  interpreted 
in  the  light  of  a  plausible  energy-band  diagram.  Conductivity 
measurements  and  the  spectral  response  allowed  us  to  make  a  number 
of  deductions  about  the  electronic  distribution  and  the  transport 
mechanism. 
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B*  EXPERIMENTAL  nrTfl 1 1  c 


a  ,The  9'"eral  exPerime"“'  details  for  the  preparation  of 
“?  3AS2SV2SV°3  fi'"s  »«•  already  been  described 

-  JC,e  r;;;  pr nin9  *? the  — - 

deposited  h  1  ms  or  electrical  measurements  were 

and  I  s  ,  eTat'°"  U,traS°"i“"^  -  Corning  02„  g.ass 
I  e  I  "  rh  3  9r°M  S'°2  '  *  thick,  insulating 

st  :  2  j  2  SUb;trateS  haJ  3  P-deposited  aluminum  or  gold 

to  the  W.  6  ”'"S  '0n9>  Wh'Ch  ““d  f°"  «  -tact 

d|  t  °r  °p  e  ectr0des'  a'ofidum  or  gold  dots  of  varying 

thus VrS-Were  evap°rated  0n  the  n,ms  deposited  on  Corning  0211  glass 
thus  forming  a  sandwich  structure  Three-mn  a-  9 

ucture.  I nree-mi 1  diameter  copper  leads 

Du;:r;;2TI|;:  the  "ith  —««■»  electronic  grade 

r  4922  S"Ver  paste-  Electrical  conductivity  measurements  in 
thw  temperature  range  of  230  to  ^7n°k’  r 

Model  t  enfi  . \  m  K  Were  performed  using  a  Statham 

Model  S06  oven  w,  th  an  accuracy  of  Thirty  to  forty  minutes 

rmre  a„  d  to  stabilize  the  temperature  before  any  readings  were 

taken.  Current  was  measured  using  a  Hewlett  Packard  Model  t,3-A  micro. 

micro-ammeter  connected  in  series  with  the  sample.  The  majo  t  e 

Z  ™S  made  "ith  a  *°''-<"°P  o'  '  V  across  the  lample 

pill  ,  "*  0bta'ned  USi"3  the  —  -tap  except  that  a  higher 

uZ  lit-  ^  Ph°t0  5peCtral  ~P~.  -s  measured 

i nco-Grat ing  monochromator  Model  #4-8400  with  a  650  W 

tungsten  movie  projector  lamp.  The  light  was  focused  on  the  back  of  the 

sample  (sOT, transparent  substrate,  mounted  inside  the  monochromator 

If  lie  data0.  ”'"ate  ^  ^  most 
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RtSULTS 


1.  Conductivity 

Figure  1  shows  a  typical  curve  of  conductivity  versus  1/T. 

Most  of  the  data  were  on  films  with  a  bulk  conductivity  of  about 
“9 

10  mhos  at  room  temperature.  The  conductance  of  these  films  was 
proport ional  to  the  area  and  inversely  proportional  to  the  thickness  as 
would  be  expected  for  a  well-defined  conduction  mechanism.  Anomalous 
effects  occurred  in  films  of  lower  apparent  conductivity.  Some  of  these 
films  showed  an  erratic  "swi thching"  behavior.  The  conductance  then 
did  not  show  a  well  defined  dependence  on  film  geometry.  The  conductivity 
follows  an  exp  -  [E/KgT]  relation  with  activation  energy  E.  The  value 
of  the  actibation  energy,  obtained  from  the  neasurements  averaged  over 
ten  samples,  E,  is  0.82  +  0.02  eV. 

2.  Field  Dependence  of  Conductivity 

Conductibity  for  relatively  high  applied  fields  (5  X  10"*  V/cm  maximum) 

1  /2 

and  at  different  temperatures  was  measured.  Semi  log  plots  of  I  vs  V  for 

1.35  Mm-thick  3As2Se^"2Sb2Se2  film  at  296.5,  288,  and 

282°K  are  shown  in  Fig.  2.  It  is  clear  that  at  296. 5°K, 

1  /2 

the  current  varies  according  to  exp[8FV  J  for  fields  greater 
4  ' 

than  about  7  X  10  V/cm.  A  similar  variation  is  observed  at  other 
temperatures.  The  intersection  of  the  three  curves,  when  extrapolated 
to  the  high  voltage  region,  yields  an  estimate  of  the  breakdown  field 
strength.  The  significance  is  discussed  later. 

3.  Photoconductivity 

The  photo  spectral  response  for  a  typical  sample  is  shown  in 
Fig.  3.  The  photocurrent  is  normalized  to  the  per  unit  incident 
power.  It  is  seen  that  thin  films  of  2Sb2Se^  are  sensitive 

photoconductors  over  the  range  of  visible  light. 
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Photoconductivity  response  per  unit  incident  power  vs 
wavelength  for  sample  p  . 


DISCUSSION 


In  order  to  interpret  the  data  on  amorphous  2As2Se  ‘2Sb  Se 
films,  it  is  necessary  to  have  an  energy-band  diagram  for  the  3 
material.  The  usual  derivation  of  an  energy-band  assumes  an  ideal  periodic 
crystalline  structure  .  However,  this  assumption  of  long-range 
order  is  not  required.  The  band  structure  of  one  and  three-dimensional 
disordered  lattices  has  been  studied  theoretically  by  Gubanov8, 

and  others  ’  .  The  detailed8  analysis  has  shown  that  the 

essential  features  of  the  crystalline  band  structure  are  retained. 

However,  the  absence  of  long  range  order  produces  tailing  (also  called 
"smearing"  or  "fuzziness")  of  the  band  edges  with  a  large  number  of 
localized  states  *  .  Consequently,  the  interval  of  localized  states 

acts  as  a  pseudogap  and  is  usually  defined  as  the  "mobility  gap"4. 


1.  Proposed  Energy-Band  Diagram 


The  proposed  energy-band  diagram  is  shown  schematically  in  Fig'.  4. 
An  appropriate  dens i ty-of-states  diagram  shown  schematically  in  Fig.  5. 

It  will  be  shown  that  this  diagram  is  consistent  with  the  experimental 
resul ts. 


2.  Low  Fields  Conduction 

According  to  the  proposed  band  model,  there  can  be  two  processes 
leading  to  conduction  in  amorphous  3As?Se3'2Sb2Se^  films:  (a) 

An  intrinsic  excitation  of  nonlocal ized  states;  and  (b)  a  hopping  process 
in  the  localized  states.  If  it  is  assumed  that  the  Fermi  level  lies 

at  the  middle  of  the  mobility  gap,  the  dc  conductivity  in  the  first  case 
vvi  1 1  vary  as 


- GAP  WIDTH  (FOR  STRICTLY 

PERIODIC  SYSTEM) 

////////  LOCALIZED  STATES 

NONLOCALIZED  STATES 


Fig.  4  Proposed  energy  band  diagram  of  amorphous 
3As2Se^*2Sb2Se2  films. 


442ev  =  E,„, 
opt 


Eopt  ABSORPTION  EDGE 
Eg(T)  MOBILITY  GAP 

VM  LOCALIZED  REGION 


NONLOCALIZED  REGION 


i 

Fig.  5  Proposed  density  of  states  for 
amorphous  3As„Se„ *2Sb_Se„  films. 


where  oq  is  a  slowly  varying  function  of  temperature,  h  is  the 

Boltzmann  constant,  T  is  the  temperature  in  degrees  Kelvin,  and 

14  15 

Eg (T)  is  the  energy  gap  which  generally  depends  on  temperature  * 

If  therefore,  we  write  Eg  (T )  =  Eg (0)  +  3T  ,  where  3  is  the  temperature 

coefficient  of  energy,  then  Eq.  (1)  becomes 

o-o  exp  -[$/2k„]  exp  - [E  (0) /2k  T]  .  (2) 

o  B  g 

-4 

Hence,  with  E^  (0)  =  2E  *  1 .64  +_  9*92  eV,  and  with  S  =  -  5.3  X  10  eV.°K 
(from  optical  data  in  Part  l),  the  value  of  the  energy  gap,  at  room 
temperature,  E^  ,  is  1.50  eV.  In  the  second  process,  the  conductivity 
caused  by  hopping  would  show  a  temperature  dependence  that  could  be 
expressed  as 

^  exp  -  [EH/k  T]  ,  (3) 

where  Eu  is  the  temperature  dependent  activation  energy  required  for 

H  i  _  i  g 

the  hopping  process.  Experimental  results  indicate  that  E^  is 

generally  very  small  (<  0.2  eV)  and  is  not  unique  in  value  since  it 

depends  on  the  hopping  lingth.  Because  our  results  show  a  unique  and  high 

activation  energy,  we  conclude  that  phonon  assisted  hopping  is  not 

significant  in  amorphous  3As0Se_ ‘2Sb,Se- ,  down  to  230°K.  Many 

*  j  j  ]  g 

amorphous  semiconductors  show  a  behavior  of  type  O  ^  C  exp  -  [E (0 ) /k  T] 

This  is  similar  in  form  to  Eq.  (2),  if  C  =  aQ  exp  -[3/k]  .  For  values 
of  energy  just  on  the  nonlocal  ized  side  of  E^.  Mott's  theoretical 

estimates  show2^  values  of  o  from  10.10^  ,  Hence  C  in  the  range  10 

- 1  - 1  °  . 
ft  cm  is  expected.  The  experimental  value  of  C,  found  <n  the  present 

work,  is  2.20  x  102  ft  '  cm  '  (see  Fig.  1).  Mott  and  Davis  2'  show  that 

if  this  mechanism  (a)  is  predominent  at  room  temperature,  the  range  AE 

of  localized  states  must  be  less  than  0.2  eV.  The  present  study  yields 

AE  =  E  -  E  =  0.06  eV  (see  Fig.  5).  Furthermore,  for  such  a  transport 
g  opt  22 

behavior  the  upper  limit  to  the  mobility  is  expressed  as 

* 

p  <_  q  a2  vgj/kT  ,  (4) 

Here  a  is  the  interatomic  distance;  vgj  is  the  electronic  frequency  defined 
by  =  Eg/h  »  Eg  is  the  energy  of  the  electron;  and  h  is  Planck's  constant. 


.-.■ATfti -At 
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Typical  values  are  v  ^  —  10  9  per  sec  and  a  =  3  The  value  of  the 
upper  limit  to  this  mobility  is  y  <_  5  cm2/Vsec.  The  deduced  value  of 
the  mobility  obtained  from  the  optical  data,  in  Part  I,  agrees 
reasonably  well  with  this  magnitude.  It  is  clear  from  the  above  results 
that  the  exponential  decrease  in  conductivity  is  due  to  intrinsic  conduction 
across  the  mobility  gap.  The  dc  conductivity  is  also  largely  insensitive 
to  the  addition  of  impurities23,  2**  a  id  this  is  further  evidence  for 
intrinsic  conduction. 

3.  High  Field  Conduction 

2 c  _  27 

It  is  a  well  established  experimental  fact  that  at  fields  in 

excess  of  some  10^  V/cm,  many  dielectric  films  as  well  as  amorphous 
semiconductor  films'8,  28,  29  exhibit  current  voltage  characteristics 
of  the  form 

I  ~  exp  [sFV1/fj 


with 


where  I  is  the  current  in  amperes,  V  is  the  voltage  drop  across  the  sample, 

e  is  the  relative  dielectric  constant,  eQ  is  the  dielectric  constant 

of  vacuum,  and  d  is  the  thickness  of  the  film.  This  is  evidence  of  either 

30 

the  Schottky  or  the  Poole- Frenkel  (PF)  effect. 


It  is  clear  that  the  curves  of  Fig.  2  obey  Eq.  (5).  The  slope  of 
the  straight  lines  increases  with  decreasing  temperature.  This  change 
is  consistent  with  the  PF  theory.  According  to  this  theory  the  slope 
is  determined  from  Eq.  (6).  For  a  typical  temperature  of  296. 5°K,  the 
value  of  this  slope  from  the  plot  of  Fig.  2  is  0.96.  Using  this  value 
in  Eq.  (6),  the  dielectric  constant  is  e  -  6.65,  which  is  a  reasonable 
value  compared  to  5.11  obtained  from  the  optical  data,  the  lines 
extrapolated  in  the  high  field  region  meet  at  V  =  Vg  -  256  V. 


This,  then,  yields  the  field  strength  for  the  electric  breakdown, 

F  =  1.96  x  106  V/ cm.  The  curves  also  indicate  that  the  slopes  of  Hn 
I  versus  1/T  have  a  tendency  of  increasing  from  negative  towards  zero 
Hence,  the  correct  form  of  Eq.  (5)  should  be 
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(7) 


This  evidence  of  a  field  enhanced  thermal  emission  from  a  localized 

level  or  trap  level  at  Vt.  The  depth  of  the  trap  level  is  determined 

from  Eq.  (7)  using  V  =  Vg,  to  be  Vfc  -  0.446  eV. 

The  trap  depth  may  correspond  to  a  maximum  in  the  statistical 

distribution  of  localized  levels.  These  levels  may  be  due  to  structural  ^  ^ 

defects  and  short  range  order  fluctuations31.  Results  of  Kolomiets,  et  al  . 
for  T£2SeAs2Te3  and  As2Se3  glasses  also  show  evidence  of  trap  levels.  The 
physical  model,  in  agreement  with  the  observed  high  field  transport  behvior, 
for  3As2Se3‘  2Sb2Se3  films  is  schematically  displayed  in  Fig.  6. 

4.  Photoconductivity 

The  spectral  response  of  the  photoconductivity  is  displayed 
in  Fig.  3.  There  is  no  wavelength  at  which  the  response  vanishes  completely. 

In  such  a  case  the  band  gap  is  defined3**  as  the  energy  at  which  the 

response  has  fallen  to  half  value.  This  yields  Eph  =  1.70  +0.02  eV 

as  the  photoconductivity  gap,  and  implies  an  intrinsic  gap  because 

the  absorption  due  to  localized  states  would  not  give  rise  to  photoconductivity. 
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High  field  emission  model. 


W 


6 


V. 


SUMMARY  AND  CONCLUSIONS 


The  results  of  the  electrical  investigations  of 

2Sb2Se^  evaporated  amorphous  films  may  be  summarized  as  follows: 

1.  The  observed  dc  conduct ibity  was  found  to  be 
consistent  with  intrinsic  semiconducting  behavior. 

2.  A  thermal  activation  energy  (mobility  gap)  at  room 
temperature  Eg  =  1 . 50  eV  was  found. 

3.  The  high  field  conductivity  was  found  to  obey  the 

field  enhanced  thermal  emission  (PF  law)  ;from  a  deep  trap  at 
0.446  eV. 

4.  The  dielectric  constant  estimated  from  the  PF-mechanism 
agreed  reasonably  well  with  that  determined  from  the  optical 
data. 

5-  The  photoconductivity  was  observed  to  be  intrinsic, 

yielding  the  band  gap.  Eph  »  1.70  (+0.02)  eV  at  X]/r 
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j .  ten... 


The  important  optical  and  electrical  parameters  of 
amorphous  3As2Se2*2Sb2Sej  films  obtained  during  this 
investigation  are  compiled  in  Table  I. 


Table  I.  Summary  of  the  measured  electrical  and  optical 
constants  of  amorphous  BAs^e^ *2Sb2Sej  films. 


Resistivity  (ohm-cm)  at  300°K 

=  109 

Mobility  (cm2  V  *  sec  ')  at  300°K 

<  5 

Optical  band  gap  (eV)  at  300°K 

1 .442 

Optical  band  gap  (eV)  at  77°K 

Temperature  coefficient  of  optical 

1.568 

band  gap  (eV/°K) 

-5.3  x  10-4 

Electrical  band  gap  (eV)  at  300°K 

1.50 

Breakdown  strength  (V/cm) 

=  2  x  I06 

Photoconductivity  gap  (eV)  at  300°K 

1.7  (+0.02) 

Refractive  index  (  at  5  M  and  300°K) 

2.26 

Dielectric  constant  (static) 

5.11 

Relaxation  time  (sec) 

1.63  x  10~11 
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guinier  x-ray  cameras 

REACTION  KINETICS 


By 


D.  E.  Clark,  G.  J. 


Scott  and  L.  L.  Hench 


Introduction 

Many  investigators  have  examined  the  kinetics  of  solid 
state  reactions.  The  methods  employed  usually  consist  of 
either  an  optical  microscopic  determination  of  reaction  layer 
thickness  in  diffusion  couples'1'2'  „r  quantitative  x-ray  anal¬ 
ysis  using  standard  x-ray  diffractometer  methods.”-6'  Both  of 
these  methods  limit  the  analysis  of  the  reacting  materials  to 
e  later  stage  of  reaction  where  the  fraction  of  product  formed 
S  greater  than  about  3  weight  percent.  However,  in  many  areas 
o  current  research  interests,  the  identification  of  less  than 
percent  of  a  phase  is  crucial  to  the  characterisation  of  the 
microstructure,  properties  and  behavior  of  ceramic  materials. 

Several  examples  of  the  areas  where  the  examination  of  the 
products  of  early  stages  of  reactions  are  important  include: 

evitrification  of  glasses,  metastable  phases  and  nuclei  in 
glass-ceramics,  interfacial  compounds  in  composite  materials, 
oxidation  products,  grain  boundary  compounds,  sintering  addi¬ 
tives,  precipitates  from  solid  solution,  fission  products  and 
inorganic  compound  nucleation  in  biological  systems 


For  the  purposes  of  this  paper,  early— stage  reactions  are 
defined  as  those  that  involve  nucleation,  surface  diffusion 
and  interfacial  rearrangements  and  which  are  not  lattice  diffu¬ 
sion  controlled;  consequently  the  fraction  of  reaction  product 
involved  is  quite  small.  Later— stage  reactions  are  derined  as 
those  that  involve  lattice  diffusion  processes  and  therefore 
the  fraction  of  reaction  product  is  usually  appreciable. 

The  objective  of  this  paper  is  to  establish  a  method  for 
the  quantitative  analysis  of  solid  state  reaction  kinetics 
through  the  use  of  Guinie r-de Wo 1 f f  and  Guinier-Lennd  x-ray  cam¬ 
eras.  Analysis  of  early-stage  solid  state  reactions  requires 
the  use  of  the  Guinier-deWolf f  x-ray  camera.  A  method  for  em¬ 
ploying  a  high  temperature  Guinier-Lennd  x-ray  camera  in  the 
analysis  of  later— stage  solid  state  reactions  is  also  presented. 
And  finally,  the  range  of  applicability  of  the  Guinier  cameras 
in  comparison  with  other  x— ray  techniques  is  experimentally 
established . 

Solid  State  Reaction  Kinetics 

Numerous  diffusion  models  have  been  proposed  for  describ¬ 
ing  the  later  stages  of  solid  state  reaction  kinetics  involving 

two  reactants.  It  has  been  experimentally  established  by  Bran- 
(  3  ) 

son  that  the  most  realistic  of  the  various  models  was  devel¬ 
oped  by  Valensi^  and  Carter.  ^  The  model  is  unique  in  that 
it  allows  for  the  volume  difference  between  reactants  and  prod¬ 
ucts.  The  following  equation  was  used  by  Valensi  for  describ¬ 
ing  the  oxidation  kinetics  of  nickel  spheres. 

K'c  ■  srr  -  u-“>2/3  -  ITT  [1  +  <4-l)a]2/3  (1) 

where  K'  =  the  reaction  rate  constant  for  a  given  temperature, 
t  =  the  reaction  time, 

A  =  the  volume  of  product  formed  per  unit  volume  of 
reactant  consumed,  and 
a  =  the  reacted  fraction  of  the  reactants. 


(3) 

Branson  has  shown  that  equation  (1)  is  valid  for  describing 
the  later-stage  kinetics  of  the  ZnO-Al^  system.  His  data, 
when  plotted  on  a  graph  of  K't  versus  t,  yielded  a  curve  with 
constant  slope  that  extrapolated  through  the  origin  for  t  =  0. 
The  significance  of  these  results  is  that  K'  is  indeed  a  con¬ 
stant  for  a  given  temperature  and,  furthermore,  the  initial 
conditions  are  satisfied  because  initially  no  product  is  indi¬ 
cated  to  be  present.  Other  models  fail  to  satisfy  both  of 
these  requirements.  The  Valensi  model  has  also  been  used  suc¬ 
cessfully  by  Duncan  and  Stewart(10)  for  describing  the  kinetics 
of  zinc  ferrite  formation.  For  these  reasons  the  Valensi-Carter 
model  will  be  used  throughout  the  development  of  the  later-stage 
reaction  techniques  in  this  paper. 

Early-stage  solid  state  reactions  have  not  been  examined 
as  extensively  as  later-stage  reactions.  This  is  primarily 
because  adequate  equipment  has  not  been  available  for  quantita¬ 
tive  analysis  of  very  small  fractions  of  the  product  phase, 
i.e.,  <3  percent.  In  fact,  little  early-stage  reaction  kinetics 
data  have  been  collected  for  any  system  to  the  authors'  knowl¬ 
edge,  even^ though  several  early-stage  models  have  been  pro¬ 
posed.  ’  For  this  reason,  the  authors  suggest  that  selec¬ 
tion  of  an  appropriate  kinetics  model  is  an  integral  part  of 
early-stage  kinetics  evaluation.  Thus,  in  later  discussions 
"early-stage  model"  is  used  in  its  general  sense. 


Description  of  the  Two  Guinier  Cameras 

X-ray  identification  of  crystals  produced  in  the  early 
stages  of  a  reaction  sequence  can  be  obtained  through  the  use 
of  a  vacuum  Guinier-deWolf f  camera.  The  camera  is  an  evacuated 
powder  camera  with  a  monochromatic  x-ray  beam  incident  on  the 
specimen.  The  focused  beam,  specimen  and  film  are  arranged  so 
as  to  produce  a  focusing  circle  (Figure  l).^11^  The  film  is 
on  the  circumference  of  this  circle,  and  a  7-inch  filmstrip 
subtends  the  arc  between  approximately  -1°  and  45°  0.  The 
focusing  circle  has  an  effective  radius  of  114.6  mm,  so  that 
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4  mm  on  the  film  corresponds  to  1°  of  Bragg  angle.  Important 
characteristics  of  the  camera  are:  (1)  high  resolution  due  to 
eliminating  the  white  radiation  background;  (2)  high  intensity 
from  a  monochromatic  beam  focused  by  a  bent  quartz  single  crys¬ 
tal;  (3)  reduced  air  scattering  background;  (4)  a  larger,  more 
representative  sample  volume;  and  (5)  four  powder  samples  are 
exposed  simultaneously  to  x-ray  irradiation. 

The  wavelength  of  radiation  on  all  four  sample  holders  is 
the  same,  but  an  intensity  decrease  on  the  two  outer  samples 
results  in  a  lesser  photographic  density  from  them.  The  photo¬ 
graphic  effects  can  be  empirically  normalized  to  those  of  the 
center  pair  by  testing  a  film  of  four  quantitative  samples  of 
the  same  substance. 


The  Guinier-Lenne  camera,  in  contrast,  holds  only  one  sam 
Pie,  but  employs  a  furnace  and  a  moving  film  cassette.  The  x- 
ray  optics  are  the  same  as  for  the  Guinier-de Wol f f  camera,  but 
the  x-rays  pass  largely  through  the  ambient  atmosphere  as  only 
approximately  3  cm  on  either  side  of  the  sample  is  evacuated. 
Both  incident  and  diffracted  beams  pass  through  beryllium  win¬ 
dows  that  transmit  83  percent  of  incident  CuKct  x-rays.  The 
furnace  is  pla t inum- rhodium  resistance  heated,  and  its  shell 
is  water  cooled.  The  film  cassette  moves  normal  to  the  focus¬ 
ing  plane  of  the  x-rays,  which  permits  reaction  effects  to  be 
recorded  on  film  as  they  occur.  Consequently,  one  of  the  most 
important  features  of  the  Lenne  camera  is  that  it  permits  the 
study  of  high  temperature  reactions  in-situ.  For  example,  the 
Lenne  camera  is  especially  useful  for  analyzing:  (1)  thermal 
expansion  vs.  temperature  curves,  (2)  temperatures  of  phase 
changes,  (3)  crystallization  of  amorphous  substances,  (4)  de¬ 
composition  reactions,  and  (5)  solid  state  reactions. 


The  temperature  of  the  Guini e r-Lennd  camera  may  be  set 
up  to  1200°C  for  isothermal  runs  for  an  indefinite  time  period, 
or  motor  driven  variacs  may  be  engaged  to  cause  the  furnace  to 
heat  from  a  preselected  low  temperature  to  any  higher  one  at  a 
controlled  rate, or  for  a  similarly  controlled  cooling  cycle. 


The  sample  holder  temperature  is  read  us-ing  a  built  in  Pt-Pt- 
lORh  thermocouple,  which  is  fed  through  the  vacuum  chamber  wall 

The  elapsed  times  between  initiation  and  completion  tem¬ 
peratures  can  be  selected  to  be  3,  4,  6,  12,  30,  40,  60  or  120 
hours.  Pates  of  film  advance  available  are  0.1,  1.0,  1.5,  2.5, 
10,  15  and  20  mm/hour,  and  slit  widths  of  1,  3  and  5  mm  can  be 
used.  Consequently,  a  large  variety  of  exposure  times  are  ob¬ 
tained  by  varying  the  slit  widths  and  the  rate  of  film  advance. 

Several  potential  complications  in  the  use  of  both  Guinier 
cameras  include:  (1)  varying  or  uncontrolled  particle  orienta¬ 
tion  can  yield  misleading  results,  especially  when  examining 
small  fractions,  and  (2)  the  presence  of  A/2  radiation  can 
cause  both  extra  lines  on  the  film  and  can  excite  characteris- 
tic  fluorescence  of  the  sample,  thus  limiting  the  exposure  time. 
However,  it  is  well-known  that  these  problems  are  not  unique  to 
Guinier  cameras  and  occur  with  other  x-ray  methods  as  well. 

The  choice  of  film  greatly  affects  the  sensitivity  and/or 
convenience  of  both  Guinier  systems.  As  can  be  noted  from 
Table  I,  very  long  exposure  times  may  he  required  to  obtain 
quantitative  data  when  very  small  percentages  of  compounds  are 
involved,  even  when  an  extremely  fast  x-ray  film  is  used.  A 
trade-off  must  be  made  between  the  available  contrast  for  a 
low- to -mode rate  photographic  density  of  the  analysis  line  over 
background,  and  the  time  reauired  to  create  its  image. 

Kodak  No-Screen  offers  a  satisfactory  compromise  for  most  appli¬ 
cations.  The  film  type  giving  the  next  significant  step  up  in 

contrast  from  Kodak  No-Screen,  Kodak  AA ,  requires  six  times 

(12) 

longer  exposures.  An  additional  advantage  of  the  smaller 

grain  size  of  the  slower,  higher  contrast  film  types  is  that 
it  gives  a  smaller  noise  amplitude  and  period,  which  assists 
interpretation  of  mi cro den  si t oie t e r  graphs. 

Although  the  films  of  cither  Guinier  camera  can  be  read 
with  standard  optical  x-ray  readers,  a  Joyce-Loebl  recording 
microdensitometer  provides  maximum  accuracy  in  analyzing  the 
exposed  film.  Such  accuracy  is  especially  important  in  early— 
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stage  kinetic  studies.  This  t»o-bea.,  continuous  w-dge  instru¬ 
ct  features  accurately  controllable  specimen  beam  intensity, 
and  electronically  monitored  scanning  speed  for  optimum  pen 
response.  A.  set  of  levers  produce  adjustable  graph-distance 
to  film-distance  ratios.  For  instance,  the  breadth  of  a  peak 
produced  on  the  densitometer  output  graph  will  be  100  times  the 
breadth  of  the  same  peak  on  the  film  when  the  100  lever  is  used. 
Lever  ratios  of  1,000:1,  200:1,  100:1,  50:1,  20:1,  10:1,  5:1, 

2:1  and  1:1  are  available.  Examples  of  a  1:1  lever  ratio  trace 
and  a  200:1  trace  from  a  film  produced  in  a  Guinier  deWolff 
era  are  shown, in  Figure  2. 


Quantitative  Sample  Preparation 

.  The  diffracted  intensity  equation  for  Guinier  cameras  given 

(13)  . 

by  Sas  and  deWolff  is 

2  (2) 

I  =  LABS  G  TkpF  P 


where 


L  = 


,  the  Lorentz  factor  for 

5  J 


(3) 


sin20  cos0  ’  Guinier  cameras,  and 
A  is  the  absorption  correction  factor  which  is  expressed  as 

1  cos'!)  •  cosX  rexp  (-yd/cosij))  -  exp  (p d/cosX)  ]  (4) 

A  yd  cos^-cosX 


for  X  =  I  20  - 


where  26  =  Bragg  diffraction  angle, 

^  -  30°  for  commercial  Guinier  cameras, 

„  .  the  linear  absorption  coefficient  for  x-rays,  and 

d  =  the  sample's  thickness. 

In  equation  (2)  B,  the  absorption  factor  for  the  film’s  pape 
wrapper,  is 

(  6 

B  =  exp [-y pdp / tosX ] . 
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The  absorption  by  the  Be  dif f racted-beam  window  in  the  Guinier- 
Lennd  camera  furnace  body  is  also  calculated  from  formula  (6), 
as  should  be  the  atmospheric  absorption  between  the  window  and 
the  film. 

In  equation  (2)  the  second  emulsion  factor  for  integrated 
intensities,  S,  is 

[1  -  exp{-(p1d1/cosX) }]  [1  +  exp{-(y1d1  p2 d2 ) / cosX  }  ] 

S  [1  -  exp(-p1d1)]  X  [1  +  exp {- (y1d1+y2d2) }]  ’  ^ 

giving  the  effective  intensity  exposing  double  emulsion  films. 
The  hkl  beam  intensity  reaching  the  rear  emulsion  has  been  re¬ 
duced  by  passing  through  the  front  emulsion  and  the  base.  In 
equation  (7)  the  thickness  of  one  layer  of  emulsion  and  of  base 
are  d^  and  d2  respectively,  and  and  p2  are  the  correspond¬ 

ing  absorption  coefficients. 

The  geometric  factor  accounting  for  the  oblique  incidence 
of  diffracted  beams  onto  the  film  is  given  by 


G  =  1/cosX 


(8) 


2 

T,  k,  p,  F  and  P  used  in  equation  (2)  are  the  standard  powder 

camera  temperature,  scale,  multiplicity,  structure  and  polari- 

(13) 

zation  factors,  respectively. 

For  optimum  quantitative  use  of  the  Guinier  cameras,  it 
is  essential  to  use  high  density  samples  that  are  consistent 
in  both  weight  and  texture.  The  major  factor  involved  in 
achieving  optimal  samples  is  the  absorption  correction  factor, 
which  contains  much  of  the  sample’s  effects  on  diffracted  in¬ 
tensity,  since  it  alone  contains  the  sample's  x-ray  opacity 
("pd"  product),  see  equation  (4). 

We  specify  the  problem  by  defining  the  sample's  thickness 
as  being  the  total  powder  weight  per  unit  area  of  sample  holder, 
divided  by  the  resulting  actual  bulk  density,  i.e., 


d 


W/A 

s 

P' 


(9) 
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Next,  by  summing  each  phase’s  absorptivity  (y"),  adjusted  for 
the  weight  fraction  present,  one  obtains  a  total  absorptivity 


V 


Q 

E 

n  =  l 


E 

(  E 
i  =  l 


n 

w  . 

y,  — ) 

i  w 

n 


(10) 


In  equation  (10)  Q  and  E  are  the  total  numbers  of  phases  and 
elements,  y.  is  the  absorptivity  of  the  i^  element,  is  the 
actual  weight  of  the  i1"  element  in  the  n  phase,  and  w^  is 
the  total  weight  of  the  nC^  phase. 

Furthermore,  the  structure,  multiplicity  and  Lorentz  fac¬ 
tors  can  be  combined  into  a  scattering  coefficient  per  unit 
weight  since  these  terms  contribute  to  the  diffracted  intensity 
from  any  specimen.  The  number  of  effective  electrons  in  the 
phase,  which  do  the  work  of  diffracting  x-rays,  is  proportional 
to  the  weight  of  the  phase  in  question.  The  structure  factor 
counts  these  electrons  and  accounts  for  the  interference  the 
hkl  beam  encounters  from  the  x-ray  scattering  of  all  the  atoms 
in  the  unit  cell.  The  mulciplicity  factor  gives  the  number  of 
equivalent  planes  which  can  cause  a  given  reflection.  Thus  a 
100  beam  from  a  cubic  crystal  can  be  excited  from  [100],  [100]  , 

[010],  [OlO] ,  [001]  and  [001]  planes,  so  the  multiplicity  factor 

is  6,  The  Lorentz  ("velocity")  factor's  value  is  controlled  by 
the  rate  that  the  hkl  node  of  the  reciprocal  lattice  penetrates 
the  Ewald  sphere.  Reference  14  gives  a  complete  discussion. 

Introducing  the  Bragg-angle-dependent  structure  and  Lorentz 
factors,  of  course,  limits  our  discussion  to  the  diffracted  in¬ 
tensity  of  a  single  line.  We  thus  restrict  0  to  the  value  given 
by  arcsin (X/2dhkl) ,  and  we  restrict  the  line  with  indices  "hkl" 
to  be  the  one  who.-e  variation  of  photographic  density  with 
weight  of  product  phase  will  be  the  analytical  basis  of  kinetics 
studies . 

Thus  the  diffraction  intensity  is  proportional  to 
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for  a  mass  scattering  coefficient  om  and  a  weight  o.  material 
Whkl  P^Perly  oriented  to  diffract  the  hkl  line.  Note  that  the 
diffracted  intensity  now  depends  on  the  weight  of  diffracting 
material  in  a  complex  way  (equation  (9)).  For  a  single  phase 
sample,  the  intensity  value  is  dominated  by  the  yd  product  of 
the  denominator,  i.e.,  the  difference  of  exponentials  rises 
slowly  for  typical  Bragg  angles,  while  the  graph  of  the  whole 
function  falls  monotonically  and  approaches  the  opacity  axis 

for  large  yd.  Thus,  there  is  no  optimum  sample  weight  for 
single  phase  samples. 

For  multiphase  samples,  the  direct  proportionality  between 
intensity  and  in  equation  (11)  indicates  that  a  sample  of 

maximum  total  weight  should  be  used.  The  limits  on  the  total 
sample  weight  that  may  be  used  are  placed  by: 

1)  tne  fluffiness  of  the  starting  material  or  product. 

2)  the  presence  of  highly  absorbent  elements  (e.g.,  Pb , 
or  Fe  or  Mn  when  using  a  Cu-target  x-ray  source)  which 
would  make  a  well-packed  sample  frame  excessively 
absorbent. 

3)  A  slightly  overfull  sample  holder  requires  much  more 
mixture  handling  than  one  slightly  underfull.  The 
amount  of  extra  pushing  of  the  powder  may  be  quite 
variable,  leading  possibly  to  a  different  preferred 
orientation,  especially  if  powder  particles  are  flat 
or  otherwise  anisotropically  shaped.  In  terms  of 

equation  (11),  may  change  as  the  amount  of 

handling  increases. 

Thus  the  optimum  quantitative  sample  Is  one  with  the  maxi¬ 
mum  weight  that  can  be  easily  and  consistently  handled,  and 
whose  x-ray  opacity  is  unity  or  less.  Note  that  for  „d  -  1,  a 
single  phase  sample's  absorption  correction  factor  is  approxi¬ 
mately  1/3.  Consequently,  the  optimal  experimental  sample  will 
produce  an  intensity  which  is  at  the  leas t  one- th i rd  that  of  the 
reference  sample  which  has  aero  absorption.  If  a  sample  exceeds 
the  "„d  -  1"  criterion  it  will  'nevitably  require  substantially 


longer  exposure  times.  Such  exposure  times  may  not  be  possible 
with  the  Guinier-Lennd  camera.  For  the  deWolff  camera,  much 
longer  exposure  times  are  possible.  The  ultimate  limit  is 
determined  only  by  the  time  at  which  the  background  becomes 
darker  than  the  hkl  line  image.  Fluorescent  radiation  from  the 
sample  or  from  the  partial  vacuum  (a  typical  vacuum  is  10-50 
raillitorr)  can  contribute  to  the  background,  as  can  inelastic 
scattering  from  the  sample. 

In  order  to  maintain  optimal  quantitative  information  from 
the  x-ray  films  it  is  essential  that  a  consistent  method 
of  storage  and  developing  be  maintained. 

Method  of  Analysis 

This  section  presents  a  general  method  for  analyzing  data 
obtained  from  Guinier-deWolf f  films.  A  schematic  of  this  method 
is  shown  in  Figure  3. 

Figure  3a  is  a  plot  of  film  line  density  as  a  function  of 
sample  reaction  time.  Samples  are  fired  at  various  temperatures 
and  Guinier-deWolf f  films  are  made  from  these  specimens.  The 
exposure  time  for  these  films  are  the  same  as  for  the  standard 
films  in  order  to  eliminate  exposure  time  as  a  variable.  The 
densitometer  is  used  to  obtain  line  density  as  a  function  of 
reaction  time.  Figure  3b  shows  a  graph  of  the  film's  hkl  line 
density  as  a  function  of  weight  fraction  of  product  phase  pres¬ 
ent  in  a  material.  This  graph  is  obtained  by  making  standards 
with  known  weight  fractions  of  second  phase  present.  Guinier- 
deWolf  f  films  are  made  by  irradiating  the  standards  for  a  fixed 
period  of  time  at  constant  milliamperage  and  voltage.  The 
microdensitometer  is  employed  to  determine  the  relative  line 
densities  of  the  films  for  the  various  weight  fractions. 

The  graph  in  Figure  3c  is  obtained  from  one  of  the  kinetic 
reaction  models  presented  in  the  literature,  as  discussed  above. 
The  equations  describing  these  models  usually  relate  K't  to  a, 
where  K'  is  the  reaction  rate  constant  and  a  is  the  percent  of 
reaction  that  has  been  completed.  The  tables  relating  K't,  a 
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and  weight  fraction  values  for  various  models  are  easily  gener¬ 
ated  by  computer. 

The  graph  in  Figure  3d  is  obtained  from  graphs  a,  b  and  c. 
The  times  for  Figure  3d  are  taken  directly  from  Figure  3a.  The 
experimental  photographic  densities  for  the  product  phase  are 
compared  with  those  of  the  "standards'  curve,  Figure  3b.  The 
weight  fractions  inferred  are  converted  to  K't  values  in  Fig¬ 
ure  3c.  These  K't  values  are  used  as  the  ordinate  values  in 
Figure  3d.  For  this  procedure,  we  have  been  considering  samples 
run  at  one  temperature,  for  the  times  which  are  on  the  abscissa 
of  Figures  3a  and  3d.  The  additional  lines  in  Figure  3d  cor¬ 
respond  to  data  obtained  at  the  first  and  two  additional  tem¬ 
peratures.  The  slope  of  each  line  gives  the  reaction  rate 
constant,  K',  for  the  particular  temperature. 

Figures  3e  and  3f  are  graphs  of  Jin  K'  versus  reciprocal 
temperature.  Figure  3e  will  be  applicable  for  the  early  stage 
reaction  kinetics  in  which  nucleation  and  surface  diffusion  may 
be  rate  controlling.  Figure  3f  is  characteristic  of  analyses 
of  later-stage  reaction  kinetics  where  diffusion  is  rate  con¬ 
trolling.  The  respective  activation  energies  are  obtained  from 
the  slopes  of  these  graphs. 

The  method  outlined  in  Figure  3  can  be  used  for  analyzing 
the  Guinie r— Lennd  films  once  the  density  values  for  the  stand¬ 
ards  in  this  camera  are  obtained.  The  values  of  the  line  den¬ 
sities  are  plotted  against  the  weight  fraction  of  the  second 
phase  to  obtain  a  standard  curve.  The  densities  of  the  hkl 
line  for  various  reaction  times  from  high  temperature  isother¬ 
mal  runs  on  unknowns  are  obtained  in  the  same  manner.  The  em¬ 
pirical  correction  for  hkl  line  density  as  a  function  of  sample 
temperature  is  then  made  as  shown  in  the  next  section. 

The  zero  time  of  reaction  is  found  by  tracing  across  a 
reactant  line  at  various  points  near  the  beginning  of  the  reac¬ 
tion.  The  point  at  which  this  line  density  becomes  zero  over 
background  is  the  zero  of  reaction  time.  The  Vernier  scale  on 
the  densitometer  allows  the  distance  corresponding  to  a  given 
reaction  time  to  be  accurately  determined.  By  knowing  the  film 
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speed  during  the  time  of  exposure  and  the  distance  .he  film  has 
moved,  the  reaction  time  for  a  given  density  can  be  determined. 

These  methods  have  been  used  successfully  to  determine 
powder  compact  reaction  kinetics  of  equimolar  NiO-Al  0  mix¬ 
tures.  The  results  will  be  published  in  a  future  paper. 

High  Temperature  Modification  of  the 
Valensi-Car ter  Equation 

A  high  temperature  modification  of  the  mass  density  of 
materials  must  be  employed  when  the  Valensi-Carter  equation  is 
used  to  express  high  temperature  data  collected  in  the  Guinier- 
Lennd  camera.  Recall  that  the  Valensi-Carter  model  pictures  a 
sphere  of  one  of  the  reactants  being  consumed  by  the  product. 
Thus,  the  volume  change  on  formation  of  product  (the  A  factor) 
is  the  one  which  exists  at  the  reaction  temperature.  Using 
room  temperature  values  for  density  of  reactants  and  products 
to  calculate  A  is  therefore  inaccurate  when  using  data  from  the 
Gumier-Lenn<5  camera.  This  factor  may  be  accommodated  in  the 
Valensi-Carter  equation  by  considering  that  A  is  a  function  of 
temperature  such  that 


Art-[1  +  3“B<T)(T°K  -  300 °K) ] 

Aht  [1  +  3ctA(T)(T°K  -  300°K)  ]  (14> 

where  Aht  =  fractional  volume  change  between  reactant  A  and 
product  B  at  the  reaction  temperature, 

Art  ""  fractional  volume  change  between  reactant  A  and 
product  B  at  room  temperature, 

aA  -  mean  coefficient  of  thermal  expansion  for  reactant 
A  over  the  range  300°K  to  T°K,  and 
aB  =  same  for  the  product  phase. 

Thus  for  an  in-situ  high  temperature  analysis  equation  (15)  must 
be  used. 


K' t 


(1-a) 


2/3 


2/3 


[l+Uht-l)«] 


(15) 


The  potential  error  associated  with  the  A  effect  can  be 
important  when  K'  is  evaluated  at  small  percentages  of  reaction 
Table  I  shows  thedependence  of  the  AhJ.  effect  on  percent  com¬ 
pletion  and  magnitude  of  A^/A^.  The  calculation  of  Table  I 
was  based  on  the  formation  of  nickel  aluminate  from  the  com¬ 
ponent  oxides  where  A^  =  1.75.  This  table  applies,  of  course, 
for  deviations  from  Ahfc  for  any  reason,  including  use  of  impre¬ 
cise  molecular  weights  and/or  densities. 

At  1 2 00 °C ,  for  the  above  reaction,  Ah  is  0.4%  greater 
than  Art,  so  the  first  column  in  the  table  applies.  It  is  ap¬ 
parent  that  only  a  small  effect  is  observed  for  the  nickel  alu¬ 
minate  system.  For  other  reactions  with  approximately  this 
Art  value>  thermal  expansion  may  cause  Ah(_  to  be  as  much  as  8% 
greater  than  A^.  Several  ranges  of  deviations,  expressed  as 
±100[(Aht/Art)-l] ,  are  the  column  heads  of  Table  I  ,  and  five 
percentages  of  reaction  completion  are  the  row  headings.  The 
body  of  the  table  shows,  in  the  upper  left  of  each  block,  the 
percentage  error  from  the  true  value  of  K't,  and  the  corres¬ 
ponding  percen tage-of-f ull-scale  error  in  the  lower  right.  The 
tabulation  shows  that  when  the  Ahfc  deviations  become  large,  the 
error  becomes  quite  noticeable  at  large  percentages  of  reaction 
completion . 

Comparison  of  Standard  X  —  ray  Techniques 
With  the  Guinier  Camera  ™ 

A  comparison  of  the  sensitivity  of  conventional  quantita¬ 
tive  x-ray  techniques  was  performed  using  Ni0-Al203~NiAl  0 
powder  mixtures.  The  purpose  of  this  experiment  was  to  compare 
the  limits  of  detection  of  various  x-ray  systems  with  those  of 
the  Guinier  cameras.  Standard  mixes  of  NiAl^  in  equimolar 
Ni0-A1203  powders  were  used  for  these  studies.  The  samples 
examined  with  an  x-ray  goniometer  were  prepared  in  the  follow¬ 
ing  manner:  0.2  gram  standard  mixes  were  spre.ad  evenly  over 
the  flat  surface  of  a  1.5  cm  diameter  silicon  wafer.  The  pur¬ 
pose  of  the  wafer  was  to  eliminate  amorphous  scattering  from 
the  sample  holder  that  would  add  to  background  intensity  and 
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thus  make  peaks  from  small  weight  percentages  of  Ni.A.1  0  diffi¬ 
cult  to  detect.  The  400  NiAl^  diffraction  beam  washed  for 
the  analysis  of  all  diffraction  studies.  The  scattering  ability 
of  this  plane  is  approximately  three-fourths  that  of  the  highest 
intensity  scattering  plane  in  NiAl^  which  could  not  be  used 
because  it  is  covered  by  a  NiO  line.  The  second  phase  was  con¬ 
sidered  detected  when  a  peak  0.3  inch  higher  than  the  highest 
peak  of  the  background  appeared  at  the  Bragg  angle  for  the  peak. 

For  quantitative  analysis  an  intensity  approximately  twice  this 
value  is  desired. 

A  standard  Philips  diffractometer  with  a  copper  x-ray  tube 
was  the  first  technique  to  be  investigated.  This  arrangement 
featured  a  graphite  diffracted  beam  monochromator  mounted  in 
front  of  a  proportional  counter.  The  monochromator  reduces  the 
background  and  sample  fluorescent  intensity  that  reaches  the 
detector.  Using  unfiltered  white  radiation,  the  smallest  quan¬ 
tity  of  NiA1204  that  could  be  detected  was  1.0  weight  percent, 
as  shown  in  Table  II  .Using  this  same  arrangement  with  a  nickel 
filter  mounted  in  front  of  the  receiving  slits,  the  detection 
limit  was  improved  to  0.7  weight  percent  NiAl^.  Quantitative 
analysis  was  possible  at  5-7  weight  percent  NiAl^. 

Another  technique  using  a  LiF  main  beam  monochromator  crys¬ 
tal  to  focus  the  copper  radiation  onto  the  sample  was  employed. 
This  method  also  involved  a  goniometer  arrangement.  The  fo¬ 
cused  x-ray  beam  should  eliminate  white  radiation  background  and 
enhance  the  x-ray  intensity  reaching  the  sample  and  the  argon 
proportional  counter.  As  shown  in  Table  II  the  limit  of  detec¬ 
tion  for  this  technique  is  no  improvement  over  the  other  dif¬ 
fractometer  methods. 

The  Debye-Scherrer  method  using  a  57.3  mm  (small)  camera 
was  investigated,  also.  Using  a  copper  x-ray  tube  and  nickel 
filter,  the  400  diffraction  line  of  NiAl^  from  a  3.0  weight 
percent  mixture  was  visible  on  the  film  after  .a  45-minute  ex¬ 
posure.  The  same  line  for  a  1.0  weight  percent  mixture  was 
not  visible  after  a  1.5-hour  exposure.  Longer  exposure  times 
were  not  feasible  because  background  radiation  from  air 
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scattering  darkened  the  film  and  covered  the  400  diffraction 
line  image. 

Table  II  also  shows  the  detection  limits  determined  for 
the  Guinier  cameras.  The  criterion  for  "detection  only"  was 
visual  detection  of  a  faint  line  at  the  position  where  the  400 
line  appears.  Illumination  was  with  a  soft,  somewhat  dimmed 
fluorescent  lamp  from  behind  the  film.  Line  position  was  con¬ 
firmed  with  the  direct  reader  supplied  with  the  camera.  A 

2 

sample  of  47;  mg/cm  ,  for  a  total  x-ray  opacity  of  0.46,  was 
used  for  all  Guinier  camera  samples.  The  fluffiness  of  the 
Linde  A  A1203  powder  prevented  using  more  material.  This  sam¬ 
ple  density  corresponds  to  20  mg  for  the  Lenne  camera  sample 

frame  and  30  mg  for  the  deWolff  camera.  Cu  K  x-radiation  with 

a 

Kodak  No-Screen  was  used  with  operating  conditions  of  40  KV  and 
20  ma . 

For  quantitative  analysis  of  second  phases  it  is  necessary 
to  take  a  microdensitometer  reading  of  the  diffraction  line. 

The  sensitivity  of  the  instrumental  reading  is  less  than  that 
of  the  human  eye  and  therefore  the  detection  limits  for  semi- 
quantitative  and  quantitative  analyses  are  a  larger  weight 
percent  of  second  phase,  as  shown  in  Table  III.  The  microdensi¬ 
tometer  gave  signal-to-noise  amplitude  ratios  of  less  than  ap¬ 
proximately  1.5/1  for  the  Detect  Only  figures.  The  criterion 
for  "Semiquantitative  Analysis"  was  signal-to-noise  ratios  be¬ 
tween  2-5.  Films  showing  ratios  of  5  or  greater  were  judged 
suitable  for  quantitative  analysis. 

The  results  shown  in  Table  Ulshow  that  the  Guinier- 
deWolff  camera  offers  only  a  factor  of  about  six  improvement 
in  detection  level  for  quantitative  analysis  over  the  best  dif¬ 
fractometer  methods.  However,  over  a  factor  of  30  increase  in 
sensitivity  for  second  phase  detection  without  the  use  of  the 
microdensitometer  is  obtained.  Thus,  the  major  attraction  of 
the  Guinier-deWol f f  camera  is  for  qualitative  detection  of 
second  phases.  It  is  suitable  for  quantitative  early-stage 
kinetic  studies  only  with  sufficient  care  in  operation. 
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The  Guinier-Lennd  camera  also  exhibits  a  factoi  of  two 
increase  in  sensitivity  in  second  phase  detection,  in  addition 
to  its.  important  advantage  of  in-situ  reaction  analysis. 
Analytical  criteria  for  Guinier-Lennd  camera  films  are  the 
same  as  for  Guinier-deWolf f  films. 

Summary 

It  is  demonstrated  that  conventional  methods  of  analyzing 
kinetics  data  are  inadequate  for  early-stage  reactions.  A 
method  for  graphical  determination  of  reaction  rate  constants 
is  proposed,  from  which  experimental  activation  energies  for 
both  early-stage  and  later-stage  reactions  can  be  obtained. 

With  the  aid  of  a  Guinier— deWolf f  camera,  one  can  study  the 
early-stage  reactions  (<3%)  involving  nucleation,  surface  dif¬ 
fusion  and  other  phase  boundary  reactions. 

The  Guinier-Lennd  camera  is  useful  for  performing  in-situ 
analyses  of  solid  state  reactions  when  the  high  temperature 
modification  of  the  Valensi-Car ter  equation  is  employed. 

Thus,  a  complete  and  precise  kinetics  analysis  is  possi¬ 
ble  with  the  joint  use  of  the  Guin ie r-deWol f f  and  Guinier-Lennd 
cameras . 
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B.  A  CORRELATION  OF  THE  SHORT  RANGE  ORDER  IN  NigFe  WITH  CERTAIN  PHYSICAL 
PROPERTIES  OF  THIS  ALLOY  (W.  L.  Wilson  and  R.  W.  Gould) 


Introduction 


The  existence  of  atomic  order  in  alloys  near  the  composi¬ 
tion  Ni^Fe  was  first  proposed  by  Dahl  (1936)  ,  in  a  Study  of  the 
effects  of  cold  work  and  heat  treatment  on  Ni-Fe  alloys.  Later 
x-ray  diffraction  studies  by  Leech  and  Sykes  (1939)  confirmed 
the  formation  of  an  ordered  structure  of  the  Cu^Au  type  upon  an¬ 
nealing  below  500°C.  Several  more  recent  investigations  have 
shown  a  strong  correlation  between  the  amount  of  order  present 
in  the  system  and  various  physical  properties.  Examples  of  these 
studies  are  the  calorimetric  work  by  Iida  (1952,  1954,  1955),  the 
mechanical  property  studies  by  Vidoz,  Lazaveric,  and  Cahn  (1963), 
and  the  magnetic  studies  by  Bozorth  (1953).  In  each  of  these 
studies,  a  strong  correlation  between  atomic  order  and  physical 
properties  has  been  shown. 

While  x-ray  scattering  measurements  provide  the  most  com¬ 
plete  information  concerning  atomic  order  in  an  alloy,  the  x-ray 
data  on  this  system  have  been  quite  sparse,  especially  during  the 
early  stages  of  the  order  reaction.  This  lack  of  data  is  mainly 
attributable  to  the  extremely  low  levels  of  order-modulated  in¬ 
tensity  obtained  from  alloys  of  this  system.  The  reason  for  these 
low  intensity  levels  is  found  in  the  equation  given  by  Cowley 
(1950),  which  may  be  written  for  the  cubic  case  as: 

[1] 

2 

I(h^,h2>h3)  =  Z  E  E  NX^Xg(f^-fg)  ct£mn  cos2irh^£  cos27Th2m  cos2TTh3m 
order  *  m  n 

In  this  equation,  XA  and  Xg  are  the  atomic  fractions  and  fA  and 
fB  are  the  atomic  scattering  factors  of  A  and  B,  respectively. 
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k  -  2tt  (S - S0 )  / A ,  and  £,  m,  and  n  are  related  to  the  interatomic 
vector  R£mn  =  £a^/2  +  +  na^/2  (a^  =  normal  crystal  axes). 

The  quantities  h^,  1^  ,  and  h^  are  continuous  variables  in  the 
reciprocal  lattice.  The  term  a£mn  is  the  short  range  order  par¬ 
ameter,  and  is  defined  as: 

a£mn  "  1  ‘  P£mn/XB 
BA 

where  P£mn  is  the^conditional  probability  of  finding  an  atom  of 
Type  B  at  a  site  R£mn  from  an  origin  A  atom.  From  Eq.  1,  it  can 
be  seen  that  elements  whose  atomic  scattering  factors  are  simi¬ 
lar  will  produce  low  levels  of  order  modulated  intensity  since 
the  term  (f^-fg)  will  be  small.  One  means  of  overcoming  this 
problem  is  the  use  of  the  anomalies  in  the  scattering  factors  of 
Ni  and  Fe  which  occur  when  using  CoKa.  Using  the  values  of  the 
dispersion  corrections  for  CoKa  radiation  for  Ni  and  Fe  given  by 
Cooper  (1963)  ,  the  value  of  this  term  was  found  to  be  approxi¬ 
mately  double  that  found  with  copper  radiation.  An  additional 
gain  in  intensity  was  obtained  by  using  a  pyrolytic  graphite 
primary  beam  monochromator.  Renninger  (1954)  has  shown  that  a 
natural  graphite  monochromator  could  yield  as  much  as  a  600%  in¬ 
crease  in  diffracted  intensity  over  a  LiF  crystal,  and  later  work 
by  Gould,  Bates,  and  Sparks  (1968)  showed  that  a  flat  crystal  of 
pyrolytic  graphite  showed  a  300%  intensity  increase  over  LiF  in 
the  same  wavelength  range.  With  the  combination  of  the  Co  radi¬ 
ation  and  graphite  monochromator,  the  intensity  of  the  diffuse 
scattering  was  brought  to  a  measureable  level. 


Experimental 


The  short  range  order  a  parameters  were  obtained  by  determ¬ 
ining  the  distribution  of  diffuse  x-ray  intensity  from  a  single 
crystal  of  NijFe  and  performing  the  Fourier  transform  of  Eq.  1. 


In  terms  of  the  fee  Ni^Fe  structure  this  becomes 


£mn 


I (h^ ,h2 »hj) 

hr0  h2  =  0  h3  =  0  NW£A-V 


1 

=  Z 


1 

z 


1 

z 


[2] 


cos27Th^£  cos2irh2m  cos2iTh2m 
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The  mapping  of  the  intensity  distribution  throughout  the  entire 
first  cell  of  the  reciprocal  lattice  would  require  enormous 
amounts  of  time.  However,  symmetry  elements  in  the  intensity 
distribution  allow  a  much  smaller  volume  to  be  measured,  with 
the  resultant  intensities  reflected  to  fill  the  first  cell.  The 
minimum  volume  of  measurement  necessary  for  the  fee  case  as  given 
by  Sparks  and  Borie  (1965)  is  shown  in  Fig.  1. 

The  intensity  distribution  within  this  tetrahedral  volume 
was  obtained  by  using  a  General  Electric  Single  Crystal  Orienter 
in  conjunction  with  a  modified  G.E.  XRD-5  diffractometer.  Since 
the  equipment  available  allowed  only  one  drive  (29)  to  be  auto¬ 
mated,  it  was  decided  to  collect  data  on  a  series  of  planar  areas 
lying  within  the  volume  studied.  This  data  was  then  extropolated 
by  digital  computer  methods  onto  a  cubic  subgrid  of  points  within 
this  volume. 

Conversion  of  diffuse  intensity  into  absolute  units  was  done 
using  polystyrene  according  to  Sparks  and  Borie  (1965)  .  In  the 

final  computer  calculation,  the  corrected  diffuse  x-ray  intensi¬ 
ties,  divided  by  the  Laue  monotonic  unit  X^XgCf^-fg)  ,  were  re¬ 
flected  through  symmetry  relationships  to  fill  the  first  cell  of 
the  reciprocal  lattice  and  Eq.  2  was  integrated  using  a  three- 
dimensional  Simpson's  rule  method.  These  calculations  resulted 
in  the  short  range  order  a  parameters  characteristic  of  the  state 

of  order1'*'. 

The  reliability  of  the  three-dimensional  interpolation  and 
intensity  reflection  routines  was  tested  by  calculating  a  series 
of  a  parameters  from  an  assumed  model  of  short  range  order  and 
synthesizing  intensity  data  at  the  points  to  be  measured,  using 
Eq.  1.  This  synthesized  intensity  data  was  then  processed  through 
the  computer  programs  and  the  resulting  a's  compared  to  the  input 
values.  This  test  resulted  in  the  a's  agreeing  within  31,  the 
major  cause  of  error  being  the  small  number  of  input  a  s  and  the 
rapidly  forced  termination  of  the  intensity  series. 

The  effects  of  temperature  diffuse  scattering  were  not  in¬ 
cluded  in  the  calculation  of  the  a  parameters  in  this  study.  The 
value  of  the  diffuse  intensity  around  the  Bragg  nodes,  where  the 


the  sufficient  intensity  data  to  obtain  the 


.  1S  Concentrated,  was  extrapolated  to  zero  at  the  nodes.  The 
size  effect  including  Huang  scattering  were  also  neglected.  The 
similarity  in  atomic  size  and  chemical  behavior  for  Ni  and  Fe 
along  with  the  absence  of  a  peak  shift  at  the  (100)  and  (110)* 
positions,  suggested  that  the  size  effect  was  minimal. 

The  x-ray  measurements  and  diamond  pyramid  microhardness 
measurements  were  made  on  the  same  single  crystal  sample.  This 
sample  consisted  of  a  disc  approximately  1"  diam.  X  1/16"  thick, 
cut  to  a  (211)  orientation.  Analysis  in  our  laboratory  showed  ’ 
the  composition  to  be  74/26  Ni/Fe.  Measurements  of  the  electri¬ 
cal  resistivity  and  magnetic  coercive  force  were  performed  on 
polycrystalline  wires  of  composition  75.151  Ni,  balance  Fe 
Both  of  these  tests  were  performed  by  the  standard  methods 'for 
wire  specimens  at  room  temperature.  For  both  the  single  crystal 
and  polycrystalline  samples,  annealing  was  carried  out  in  vacuo 
at  a  temperature  of  480°C,  approximately  10°C  below  the  critical 
ordering  temperature,  for  varying  amounts  of  time. 

Results  and  Discussion 

Following  the  procedures  outlined  above,  x-ray  and  micro¬ 
hardness  data  were  taken  for  the  single  crystal  sample  as  quenched 
rom  1000  C  and  at  annealing  times  of  20  min.  to  40  hours  at  480°C. 
The  values  of  the  a  parameters  for  the  first  several  co-ordination 
shells  as  a  function  of  annealing  time  at  480°C  are  shown  in  Fig. 

2,  along  with  the  values  for  the  perfectly  ordered  structure 

In  the  present  work,  it  was  found  that  an  error  of  less  than 
one  count  per  second  in  the  measured  background  intensity  could 
change  the  value  of  aQ00  from  greater  than  1.5  to  1.0,  while  the 
other  a ' s  were  scarcely  effected. 

Referring  again  to  Fig.  2,  one  should  note  that  for  anneal¬ 
ing  times  up  to  8  hours  at  480°C,  there  are  several  discrepancies 
m  the  signs  of  some  of  the  a's  as  compared  to  perfect  order. 

This  effect  has  been  found  by  Cowley  (1950)  in  his  work  on  Cu  Au 
and  has  been  discussed  as  showing  the  order  relationships  to  hold 
for  the  most  part  only  over  nearest  neighbor  distances.  The  er- 
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Perfect 


t(  hr ) 

Fig.  2:  The  short  range  order  a  parameters  for  Aran  =  110,  200, 
112,  and  222  as  a  function  of  annealing  time  at  480°C. 
Note  the  "wrong"  signs  for  06^12  a222  ^uri-n£  the 

early  stages  of  the  annealing. 


ror  in  sign  stems  from  the  tendency  for  each  atom  to  surround 
itself  by  unlike  nearest  neighbors  during  the  early  stages  of 
order,  a  pattern  which  is  not  entirely  compatible  with  the  fee 
Cu^Au  type  structure.  At  times  greater  than  12  hours,  the  signs 
of  the  ot's  can  be  seen  to  progress  toward  those  found  in  the  long 
range  ordered  structure,  although  their  magnitudes  show  that  the 
arrangement  is  still  far  from  perfect  order. 

The  results  of  the  microhardness,  electrical  resistivity, 
and  coercive  force  measurements  also  show  a  change  to  take  place 
at  roughly  12  hours  of  annealing.  This  is  shown  in  Fig.  3.  The 
rate  of  change  of  hardness  with  time,  although  the  least  sensi¬ 
tive  of  the  measurements  performed,  shows  a  definite  break  to  oc¬ 
cur  between  10-14  hours  of  annealing.  In  a  similar  manner,  the 
rate  of  change  of  resistivity  with  time  shows  a  major  change  dur¬ 
ing  this  time  period.  The  coercive  force  data  show  that  the  co- 
ercivity  decreases  with  time  upon  annealing  at  480°C  for  times  up 
to  12  hours.  After  reaching  a  minimum  value  at  this  time,  the 
coercivity  then  increases  with  further  annealing. 

The  change  in  the  ordered  structure  from  one  dominated  by 
nearest  neighbor  interactions  to  one  with  longer  range  correla¬ 
tions  shows  an  excellent  correlation  with  the  proposed  sequence 
of  ordering  given  by  Iida  (1955).  Based  upon  his  calorimetric 
data,  he  proposed  that  the  period  of  annealing  from  the  quenched 
state  up  to  approximately  10  hours  produced  a  short  range  ordered 
structure,  with  a  small  amount  of  long  range  order.  The  period 
lasting  from  5-60  hours  of  annealing  was  proposed  to  include  the 
formation  of  the  long  range  ordered  structure.  This  would  seem 
to  agree  quite  well  with  the  explanation  of  the  signs  of  the  a's 
given  earlier. 


Conclusions 


Through  the  use  of  x-ray  diffuse  scattering  methods,  the 
local  atomic  arrangements  during  the  early  stages  of  the  order 
reaction  in  Ni^Fe  have  been  quantitatively  determined.  These 
data  have  shown  that  for  annealing  periods  less  than  12  hours  at 
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480°C,  the  nearer t-neighbor  interactions  dominate  the  distribu¬ 
tion  of  the  atoms  on  the  lattice  sites.  For  longer  times  of  an¬ 
nealing,  this  gradually  changes  to  a  more  cooperative  interaction 
extending  over  longer  distances  and  approaching  the  arrangement 
found  in  the  fully  ordered  structure. 

The  changes  in  the  hardness,  resistivity,  and  coercive  force 
which  occur  at  or  near  12  hours  of  annealing  support  the  conten¬ 
tion  that  a  change  in  the  type  of  atomic  arrangement  is  taking 
place  at  this  time. 
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IV.  Discussion 


The  research  of  this  ninth  semi-annual  period  has  provided  technical 
findings  that  bear  on  problems  of  interest  to  the  Department  of  Defense. 
Moreover,  from  this  research  have  come  foundations  undergirding  our 
future  studies. 

Two  sections  of  this  report  (II-B  and  C)  deal  directly  with  devices 
for  photodetection;  a  third  section  (H-D)  deals  with  fundamental  properties 
exhibited  by  germanium  doped  with  nickel,  which  provides  information  of 
potential  use  in  the  design  of  photodetection  and  other  devices  which  find 
wide  employment  in  military  applications.  To  place  into  perspective  the 
properties  of  our  new  type  of  photodiode,  we  note  that  the  fast  response 
time,  high  responsivity,  flat  spectral  response  between  0.4  and  1.0  microns, 
and  low  noise  performance  reported  here  apply  to  a  device  whose  design  has 
yet  to  be  optimized.  Work  toward  optimization  progresses,  together  with 

studies  seeking  a  fuller  understanding  of  the  details  of  the  observed 
behavior. 

MOS  transistors  will  likely  form  the  heart  of  many  military  electronic 
systems  of  the  future,  and  their  current-voltage  behavior  and  noise  properties 
need  to  be  understood  for  proper  application.  Here  we  have  described  certain 
aspects  of  their  performance  resulting  from  avalanche  multiplication. 

Additional  related  studies  are  now  in  progress.  One  study  seeks  a  more 
detailed  understanding  of  the  noise  behavior  attending  avalanche  multiplication 
in  semiconductor  devices.  Another  will  describe  the  alteration  in  the  noise 
performance  of  junction  field-effect  transistors  after  exposure  to  an 
irradiation  environment. 

Properties  exhibited  by  amorphous  semiconductors  make  them  potentially 
useful  in  a  variety  of  applications  that  may  involve  exposure  to  radiation. 

In  this  report  we  have  described  certain  basic  experiments  that  lend  an 
understanding  of  the  behavior  of  these  materials.  Our  future  effort  will 
concentrate  on  determining  the  radiation  sensitivity  exhibited  by  amorphous 
semiconductors  and  on  assessing  the  ultimate  limitations  and  promise  of 
these  materials  in  device  applications. 
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